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FORSWORD 


This second publication of the Space Shuttle GN&C Design Equation Document 
contedns baseline equations for approximately fifty percent of the GN&C 
computation requirements as specified in the GN&C S/W Functional Requirements 
Document (MSC-O 369 O Rev. B) . This document supercedes the original MSC-04217 
and the subsequently published revision. Additions or corrections to this 
document since its original publication are indicated in the Table of Contents 
by asterisks in the margin. 

It is planned to republish this document in a new revision in approximately 
four months time. At that time it is anticipated that equations will be 
available for virtually all reqiiirements. The new revision will be issued 
with format changes intended to stress 5.nterdependency of related submittals 
and to eliminate duplication to the greatest degree practicable. 

This issue has been modified to reflect the shuttle-structure cmd avionics- 
configuration changes which have occurred subsequent to the first issue. A 
significant change is that orbiter control of the booster has been added as 
a requirement. Decentralization of the ccc^utations and allocation to sub- 
systems is the current trend with the MARK I & MARK II shuttle configurations. 

The confutation requirements for shuttle vehicles and missions may be much less 
than those allowed for in this document. However, since the configurations are 
very fluid at this state in the shuttle development, the approach adopted in 
this document is to include as conflete a set of design equations as possible 
to cover reasonable possibilities. Therefore, subsets of equations may be ex- 
tracted from this document to fora specifications for specific vehicles, com- 
puters and missions. 

The GN&C Design Equations document is the result of the efforts of many people 
frcm NASA and sifport contractors. The list is too long to credit all con- 
tributors; however, contractors which made direct contributions to the document 
are as follows: 

a. TRW Systems Groif, Inc., Houston Operations 

b. Mrr/Charles Stark Draper Laboratory 

c. Lockheed Electronics Co., Inc., Houston Aerospace Systems Division 

d. The Boeing Co., Houston, Texas 

The equations are reviewed by the GN&C Formulation and Iiflementation Panel 
and their comments Included on submitted forms where appropriate. The names 
of equation submitters are included on the submittal sheet in each section. 
Comments on the submltteds should be referred to the individi^ submitter or 
to the responsible NASA engineer. General comments on the document or proposed 
submittals should be referred to the System Analysis Branch, Guldemce and Control 
Division. 
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1. PURPOSE 

The purpose of this document is to specify the equations 
necessary to perform the guidance, navigation and control onboard 
compute tion functiens for the space shuttle orbiter vehicle. This 
equations document will provide as Comprehensive a set ol equations 
as possible from which modules may be chosen to develop Part I 
Specifications for particular vehicles, computers and missions. 
This document is e^cpected to be the source of any equations used 
to develop software for hardware/sof tware feasibility testing, for 
ground-based simulations or flight test demonstrations. 
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ciCOPS 




This document defines a baseline set of equations vhicn 
fulfill the computation requirements for guidance, navigauion 
and control of the space shuttle orbiter vehicle. All shuttle 
mission phases are covered from Prelaunch through Landing/Rollout. 
The spacecraft flight mode and the aircraft flight mode are ad- 
dressed, Equations f ’e included for the Mark I systems and Mark II 
systems through the all-up shuttle configuation. Control of the 
booster during launch is covered. The baseline equations may be 
implemented in a single GN&C cos^uter or may be distributed among 
several subsystem computers, depending upon the outcome of cen- 
tralization/decentralization deliberations currently in progress. 


2-1 



3. APPLICABILITY 


Tills document Is applicable to the guidance) navigation snd 
control (GN&C) computation functions for the space shuttle orbiter 
vehicle. It specifies a set of baseline design equations which 
may be used for the shuttle program softveure specification 6ind 
hardware sizing. It defines the baseline equations for MSC G&CD 
hardwarc/sof tware simulation. 
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9. DESCRIPTIONS OF EQUATIONS 


The detailed equations for the GN&C functloc*^ are defined in this 
section. The organization of this section is tentative and will be 
modified so as to present the equations as they are designed in as clear 
a fashion as possible. As an Introduction to each major subsection 
(usually a mission phase) , the general software functions to be 

implemented will be identified and, where appropriate, a conceptual 
discussion and top level flow of the computations, inputs and outputs 
will be included in order to understand and sumnarize what is to be 
covered. This should be an order of magnitude less detailed them the 
flow diagrams of the equations which come later. 

A Qi&C Equation Submittal sheet will introduce each of the GN&C 
equation submlttetls and summarize the GN&C functions, and identify the 
source and NASA contact for each. 

The detailed data to be presented for each (31&C function within each 
of the major subsections (usually a mission phase) is summarized below. 
Although items 6 through 10 are to be referenced only in the equations 
document, they are required submittals before the equations can be 
approved «nd finalized for flight software development. 

1. Functional Bequiresmnts 

The specific functional requirements (from the QI&C 
Software Functional Requirements Document) which are 
satisfied by the equations should be identified. 

2. Functional Diagram 

A brief functional explanation and description of the 
overall concept and approach. A functional block 
diagram should be used where clarity is enhanced. 

Inputs, outputs, and interfaces will be provided. 

3. Equati<ms and Flows 

Detailed equations and a descriptive text which guides 
the reader through the flows of Section 10 should be 
provided. The minimum frequency ef the computations 
shall be specified and rationale given or referenced. 
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A. Coordinate System 

The coordinate systems used shall be defined. 

5. Constants /Variables Summary 

Constants and variables shall be summarized In tabular 
form with the following Information: 

a. Variables /constants symbols and definitions 

b . Units 

c. Allowable quantization 

d. Range of values 

6. FORTRAN Coding 

The FORTRAN coding of the ftmctlon for verification using 
the Space Shuttle Flight Simulation (SSFS) will be 
referenced . 

7. Simulation 

The SSFS specifications, description and user's guide 
used to verify each GN&C function will be referenced. 

8. Testing 

Test plans and test results will be referenced. 

9. Derivation 

The mathematical derivation of the equations Including 
all mathematical assumptions shall be referenced. 

10. Assumptions 

The following will be referenced: 

a. Avionics baseline system assumed 

b. Reference missions assumed 

c. Vehicle mass properties assumed 

d. Propulsion models assumed 

e. Environment models assumed 

f. Error models assumed 

The major subsections of this section are identified and partially 
expanded in the following. 
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9.12 


DEORBIT AND ENTRY 


The Deorblt and Entry phase of the Shuttle mission begins with 
the Shuttle performing the Deorblt Targeting for the Deorblt (or Deorblt 
phasing. If required) bum and terminates at the beginning of transition 
to low angle of attack. The primary landing point and the required orbit 
maneuvers have been determined from pre-phase planning conq>utatlons . The 
GN&C Software functions to be performed during this phase are: 

1. Deorblt Targeting will be performed In the coast 
period prior to the deorblt bum(s) . Targeting 
outputs (for a specified landing site Input) will 
Include parameters as tlme(s) of Ignition, Initial 
attitude for the bum(s), bum time, total AV 

for each bum, IMU realignment schedule, time of 
reentry, transition time, cross and downrange re- 
quired, jet engines Ignition time, and time of 
touchdown. 

2. Orbital navigation using ground beacon or star/ 
horizon measurements as required to update Shuttle 
State prior to thrusting maneuvers and aerodynamic 
entry augmented by available external measurements. 

3. Powered Fll^t guidance resulting In comnands to 
autopilot (attitude) , thrust vector controls , or 
thrust throttle commands to achieve desired con- 
ditions (state and attitude) for transition and 
limit entry g- loading and heating. 

4. Autopilot computations to compute reaction control 
system commands to achieve the proper attitude. 

t 

I 

Figure 1 displays a functional flow diagram of the Deorblt and Entry 
Software Functions. 
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9.12.1 Debrbit Teu^getlng 


SPACE SHirrTLE 

GK&C SOFTWARE EQUATION SUBMITTAL 


Software Equation Section DeorbAt Tar^^tlng Submittal No. 28 _ 

Function To deter«nlne the deor bit man ttBYer as a function of entry range 

and flight path angle. 

Modulo No. OG-5 Function No. 1.2.A.5 (MSC-03690 Rev. A) 


Subodtted byj Brand & firenzian Co. MIT. N o. 14 

(5SS1 

Date: 24 August 1971 


NASA Contact 


C. Lively 


Approved Panel III 


(chaizmai^ 


Organization EG2 
Date jtl 


Sunaoaiy Description: This routj rw* been designed to step through 

successive solutio^ ^valld opnortupiti «f* t .9 ful fill entry conditions 
withiAnconstraints) . *11 oirin g__t^ crew to make selection based on entry 
cpossrange. tiae-to ignition, reoidred velocity change, l.s. lighting 
conditions, urgency of return, etc. 


Shuttle Configuration: (Vehicle, Aero Data, Sensor, Et Cetera) 


This software is independent of shuttle configuration. Its utility 
requires sigificant cross-range ctoabllity. 


Conments: 


(Design Status) 

(Verification Status) 
Panel Conments: 
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9.12.1 DehPfaiit Targeting (contd^i) 


1. INTRODUCTION 


The large entry crossrange capability of the shuttle permits deorbit to a 
specified landing site to be accomplished with a single maneuver. Since the re- 
quired velocity change is smallest when no plane change is made, the equations 
presented here are designed to target the Powered Flight Guidance Routines (Ref- 
erence 3) for an in-plane maneuver. The ignition time for this maneuver is se- 
lected to satisfy entry interface and landing site constraints with minimum fuel 
expenditure. 

If the shuttle had no crossrange capability, then an in-plane deorbit maneu- 
ver to a specified landing site could only occur when that landing site, which ro- 
tates with the earth, intersects the orbital plane of the vehicle. Assuming the 
landing site latitude is less than the orbital inclination angle, and neglecting the 
effects of precession, the landing site will intersect the orbital plane twice every 
twenty-four hours. However, the time difference between these two intersections 
is in general not twelve hours. In the case when the landing site latitude is equal 
to the orbital inclination there will be only one intersection every twenty-four 
hours. 

Since the shuttle has a high crosarange capability, deorbit does not require 
intersection of the landing site vector and the orbital plane. It is possible when- 
ever the angle between the landing site vector and the orbital plane is less than 
approximately 20 deg. In general, there will be two sets of opportunities 
every twenty-four hours. Within each set, there may be several deorbit oppor- 
tunities occurring on consecutive orbits with varying crossrange requirements. 
When the latitude of the landing site approaches the inclination of the orbit, these 
two sets merge to become one. It should be noted, in addition, that if the landing 
site latitude is greater than the orbital inclination, the landing site may still fall 
within the crossrange capability of the vehicle. With these facts in mind, this 
routine has been designed to continue stepping through successive solutions, allow- 
ing the crew to select a particular deorbit opportunity based upon entry crossrange, 
time-to-ignition, required velocity change, landing site lighting conditions, urgency 
of the return, etc. 

The desired entry range and flight path angle will be considered inputs to this 
routine, since available data relating to footprint size and shape, entry heating at 
various ranges . and optimal entry flight path angle are only preliminary. In future 
revisions, consideration should ti given to computing the optimum values of these 
quantities for the particular situation. 
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9.12.1 Deorblt Targeting (cont'd) 


NOMENCLATURE 


a 


^2 


d 


‘ACR 


‘CR 


DR 


ACS 


Semimajor axis 

Alarm code— failure in Av minimization loop 

Alarm code— failure in Precision Required Ve- 
locity Determination Routine 

Semimajor axis of Fischer Ellipsoid 

Estimated magnitude of the thrust acceleration 

Semiminor axis of Fischer Ellipsoid 

Number of columns of navigation filter weighting 
matrix (set to 0 in this routine since the matrix 
is not required) 

Maximum acceptable crossrange distance of Orbit^r 

Estimated entry crossrange distance 

Entry downrange distance 

Magnitude of the engine thrust 

Magnitude of the attitude control system 
translational thrust 


^OMS 


Magnitude of the orbital maneuvering 
system engine thrust 


h j,j Entry interface altitude (400, 000 ft) 
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9.12,1 Deorblt Targeting (cont'd) 


i Unit vector formed by the cross product of the 

angular momentum and the landing site vectors 

ij,j Unit vector in the direction of r j,j 

First estimate of _i 

i'j,j ^ Z -component of the unit vector (z-axis 

assumed North) 

i ^ Unit vector in the direction of the angular momentum 

i Unit vector in the direction of the landing site 

projection into the orbital plane 

_i Unit normal to the trajectory plane (in the direction 

of the angular momentum at ignition) 

k ^ Sensitivity coefficient used tc compute adjustment 

to time-of-arrival at entry interface 

m Estimated vehicle mass 

n Iteration counter 

n Iteration limit 

max 

n Integral number of complete revolutions to be made in 

the transfer (set to zero in this routine) 

Semilatus rectum of deorbit trajectory 

Secant squared of the desired entry flight path angle 

Secant squared of the offset entry angle used by the 
Powered Flight Guidance Routine 

£q Precision position vector 


PD 


P pFy 
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9.12,1 Deorblt Targeting (cont'd) 


r'^ Entry interface position from Precision Required 

Velocity Determination Routine 


-D 

-El 



- LS 


Position of the impulsive deorbit maneuver 
Entry interface position 
Position vector at ignition 

Estimated landing site position at the time of 
landing 


-PFT Powered flight offset target vector 


s^ „ Engine select switch 

eng 

Switch set to indicate non -convergence of Pre- 
cision Required Velocity Determination Routine 


s pp Switch set equal to one after the first pass through step one 

s Switch indicating which perturbations are to be 

included in the Precision State and Filter Weighting 
Matrix Extrapolation Routine (See Reference 5) 


s Switch set when the target vector must be projected 

proj ” 

into the plane defined by 

t^ Precision state vector time 

t ^ Time of impulsive deorbit maneuver 

1 2 Time-of-arrival at entry interface 

1 2 Estimated time at which in-orbit position vector is 

coincident with the landing site projection into the 
orbital plane 
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9.13.1* Deorblt Targeting (cont'd ) 


ETL 


ig 


LTL 


-0 


-2 


-D 


-El 



-PET 


RD 


Desired earliest time-of-Janding 
Ignition time 

Estimated time-of-landing 
Desired latest time-of-landing 
Precision velocity vector 

Entry interface velocity from Precision Required 
Velocity Determination Routine 

Pre- impulse velocity 

Entry interface velocity 

Ignition velocity vector 

Velocity associated with the pov/ered flight 
offset target vector 

Post -impulse radial component of velocity 


V Required velocity 

-req Required velocity on the coasting trajectory 

V jjp Post-impulse horizontal component of 

velocity 


6 1 A djustment t j At 

Ad^CR Increment added to the acceptable crossrange, 
used in rough cross range check 

Ap^ Difference between the predicted and desired p 
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9.12.1 Deorblt Targeting (cont«d) 


Arproj Out-of-plane target miss due to the projection 
of the target vector 

Atjj^ Transfer time 

At ^2 Transfer time 

At 23 Transfer time (tg-t 2 ) 

Atg Estimated duration of the powered maneuver 

At d£ Time-of-flight difference between (1) the interval 

from deorbit through entry to landing, and (2) the 
time spent in orbit over the same total central 
angle 

At jp Time-of-flight required to transfer through the 

central angle 

Av Required velocity change 

Av p Previous value of j Av j 

AS Increment in in-plane angle Sp 

ASq Initial increment in in-plane angle S^ 

Convergence criterion on Ap 
py . y 

Convergence criterion on angle AS 

y^ Post-impulse flight path angle 

Desired entry flight path angle measured 
from the horizontal 


Xls Landing site longitude 

e In-plane angle between precision state vector 

and entry interface 
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9.12.1 Deorblt Targeting (cont'd) 


='01 


“D 


'IP 




P 





T 


In-plane angle between precision state vector and 
deorbit position 

In-plane angle over which search is made to 
find minimum deorbit Av 

In-plane central angle traversed during entry 

Angle between precision state vector and the 
projection of the landing site into the orbital 
plane 

Previous value of ® 

Gravitational parameter of the earth (product of 
the earth's mass and universal gravitation con- 
stant) 

Landing site latitude 
Orbital period 
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9.12,1 Deorblt Targe tlnf f (cont'd) 


2. FUNCTIONAL FLOW DIAGRAM 


A functional flow diagram presenting the basic approach to the oeorbit target- 
ing problem can be found in Figure 3. In addition to the state vector, the primary 
inputs to the routine are the landing site location (latitude and longitude), the entry 
downrange distance, the entry angle (at 400. 000 ft) and the earliest desired time of 
landing. Since the high crossrange capability may make deorbit possible on two or 
more consecutive orbits, after each solution the crew has the option to recycle 
the program to determine the next possible deorbit opportunity. To give the crew 
the flexibility to evaluate solutions in the future without stepping through all eaurlier 
opportunities, the earliest desired lime-of-landing is included as an input. How- 
ever. the vehicle is aissumed to be in coasting flight until the deorbit maneuver, 
and therefore the effects of any maneuvers prior to deorbit are not accounted for. 

After the vehicle state vector is extrapolated forward to the earliest desired 
time-of-landing, the solution process is initiated. This consists of three major 
steps. During the first step the vehicle state is further advanced until the landing 
site, which rotates with the earth, lies sufficiently near the orbital plane so that it 
is within the crossrange (or out-of-plane) capability of the entry phase. During the 
next step an iterative process is used to select the ignition time for this oeorbit 
opportunity which requires the smallest velocity change, thus minimizing the fuel 
expenditure. Since the first two steps involve several conic approximations to 
minimize the computer time used, the ihird step fine tunes the solution by gener- 
ating a precision trajeciory which satisfies the constraint on the desired entry angle 
while accounting for gravitational perturbations and the non-impulsive nature of the 
deorbit maneuver. After completion of this step the results are displayed to the 
crew. They may then elect to accept t4ie solution, recycle the routine to solve for 
the next deorbit opportunity, or exit. If they accept the solution, a few minor 
computations are required to initialize the Powered Flight Guidance Routines for a 
modified Lambert aimpoint maneuver. 

To aid the reader in understanding the functional flow diagram, each of the 
three major steps in the solution process is discussed in more detail below. 

2. 1 Determination of the Next Deorbit Opportunity (Step 1) 

To determine the next possible deorbit opportunity, it is necessary to calcu- 
late the inertial location of the landing site (which rotates with the earth) at the 
time-of-landing. Then the angle between the orbital plane and the landing site can 
be used to estimate the crossrange required during entry. To accom.plish this, an 
estimate of the time -of-fli|^t difference ^tj^^ between (1) the interval from de- 
orbit through entry to landing, and (2) the time spent in orbit over the same 
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9.12.1 Deorblt Targeting (cont*d) 


central angle is used. Analysis has shown that a constant is probably adequate to 
represent this difference since more precise calculations ip the following step will 
compensate for any error. 


Upon completion of the initialization process, the state vector is extra- 
polated forward to ihe earliest desired time- of -landing. Then the inertial location 
cf the landing site at the present state vector time, biased by the time difference 
AtoE* is computed. This landing site vector is projected into the orbital plane, 
allowing the in-plane central angle 0jp between the vehicle position and the pro- 
jection of the landing site to be determined. 


VEHICLE 

POSITION 



Figure 1. Out-of-plane Geometry 

The conic routines can now be used to determine the time -of- flight ^t|p 
required to coast in orbit through the central angle 8^p . If the state was then 
propagated through this central angle, its position vector would be aligned with the 
previously determined projection of the landing site vector. Unfortunately, the 
landing site will move slightly due to earth rotation while the vehicle transfers 
through the central angle. Therefore, the inertial location of the landing site must 
be recomputed, accounting for the time difference Atj^g explained previously. 
Thus, an iterative process is required to precisely determine the location of the 
landing site at the expected time-of-landing. During the first pass through the de- 
orbit targeting routine, the previously described steps are repeated once to insure 
convergence. However, on subsequent passes no iteration is required, since the 
initial guess achieved by extrapolating the state vector one orbit beyond the previ- 
ous solution guarantees a small value for the time-of-flight correction Atjp. 

Assuming the deorbit maneuver is in-plane, the angle between the orbital 
plane and the landing site location at the estimated time-of-landing can be used to 
measure the crossrange required during the entry phase. If the croesrange is 
within the capability of the vehicle, the solution process continues on to the next 
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9.12.1 Deorblt Targeting (cont'd) 


step. If not. the vehicle state is extrapolated forward one revolution to the next 
potential deorbit opportunity and the process of estimating the crossrange is re> 
peated. 

It should be noted that the process used to determine the crossrange require* 
ment is only approximate, and therefore a small increment is added to the tolerance 
used in the cross range check to allow for this. A small number of cases which 
pass this check will actually lie outside the vehicle crossrange capability, how- 
ever, a more precise check later will screen these out. 

2. 2 Ignition Time Selection (Step 2) 

During this step in the solution process, an ignition time is selected which 
minimizes the impulsive velocity change required. For these computations the 
projection of the landing site into the orbital plane is assumed to be the real landing 
site. Then, based upon the desired entry downrange distance, a target position at 
entry interface which also lies in the orbital plane can be defined. This target 
position is set 400, 000 ft above the Fischer ellipsoid. 



Using this entry interface target, and the desired entry flight path angle, a 
search is made on the central angle Oj) traversed between the deorbit maneuver and 
entry interface to locate the position and time of the minimum A v maneuver. 
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Then the time -of -flight required for the deorbit and entry phases can be accurately 
determined. Using this time -of -flight, an accurate calculation of the inertial loca- 
tion of the landing site at the time -of- landing can be made, and the entry interface 
target can also be updated. Tc preserve the central angle of the deorbit phase, the 
impulsive maneuver time is adjusted. Then the ignition time is biased from the 
impulsive time by half the expected length of the maneuver and the state vector is 
extrap>olated to this time. 

Since the location of the landing site at the time>of-landing is now known 
accurately, the angle between the c bital plane and the landing site is recomputed 
to precisely measure the entry crossrange required. Then a precision check is 
made, and any solution exceeding the crossrange capability is rejected, thus 
returning the routine to step one to search for the next opportunity. 

2. 3 Precision Solution (Step 3) 


During this step a precision integrated trajectory from deorbit to entry inter- 
face is generated which accounts for both the finite length of the thrusting maneuver 
and the effects of gravitational perturbations. Since the time-of-flight from de- 
orbit to entry interface is known, the Precision Required Velocity Determination 
Routine can be used to generate this trajectory. However, the effects of conic ap- 
proximations in the previous steps and the finite length of the maneuver can cause 
significant error in the reentry angle. Therefore, the resulting entry angle is 
checked and if it is in error, a slight modification is made in the time-of-flight 
from the deorbit maneuver to entry interface to adjust the entry angle. Then the 
precision trajectory is recomputed. After satisfying the flight path angle constraint, 
pertinent data relating to the maneuver can be displayed to the crew or transferred 
to the Mission Planning Module. 


9.12-U 




9.12.1 Deorblt TargetlPf^ (cont'd) 


ENTER 



A 


Figure 3a. Functional Flow Diagram 
9 . 12 - 5 ^ 
















9J2.1 Deorblt Targeting Coont«d) 


3. INPUT AND OUTPUT VARIABLES 


Input Variables 

Time associated with the vehicle state vector 
£o* —0 State vector 

t Earliest desired time-of-landing 

t Latest desired time-of-landlng 


m 


®eng 



‘^ACR 



-PFT 

PpFY 



®proj 


Current estimated vehicle mass 
Engine select switch 
Landing site latitude 
Landing site longitude 
Desired inertial entry angle 
Desired entry downrange distance 

Maximum acceptable crossrange distance 
Output Variables 
Ignition time 

Time of arrival at entry interface and time associated 
with offset target 

Offset target vector 

Parameter defining the desired conic entry angle 

Integral number of complete 360° revolutions 
(n^^^ = 0 for deorbit) 

Unit normal to transfer plane in direction of 
angular momentum vector 

Switch indicating whether the initial and target 
vectors are to be projected into the plane defined 
by the unit normal 
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4. DESCRIPTION OF EQUATIONS 

To minimize the size of the Deorbit Targeting Routine, extensive use is 
made of other routines. Therefore, this routine consists primarily of simple equa- 
tions, logical operations, and calls to other routines. Since most of the complicated 
equations requiring detailed explanation are contained in the description of the 
other routines, this section will be limited to a list of items not covered in the text 
describing the functional flow diagram. These items will be listed in their order 
of occurence, and are intended to supplement the detailed flow diagram in sub- 
section 5. 

4. 1 Selection of Perturbing Acceleration during 

Precision State Extrapolation 

During the first step in the solution process, which may require long term 
state vector extrapolation, it is desirable to maximize accuracy by including all 
significant perturbing accelerations in the extrapolation process. Therefore, the 
switch , which controls the selection of perturbing accelerations in the 

Precision State Extrapolation Routine, is set to 2. During the later portion of the 
routine, referred to as step three, the switch is reset to 1 , thus limiting the dis- 
turbing acceleration to the J 2 term, the second harmonic of the earth's gravita- 
tional potential function. Since extrapolation during step three is limited to the 
interval from the deorbit maneuver to entry interface, the effects of smaller 
perturbing accelerations are not significant. In addition, extrapolation over this 
interval lies within an iterative loop, and thus may be repeated several times. 

The simplified model can therefore significantly reduce the running time of this 
step. 

4.2 Selection of 6 jp Quadrant 

During the discussion of the functional flow diagram, it was mentioned that 
successive solutio.is to the deorbit problem (when successive solutions exist) are 
about one revolution apart. To find succeeding solutions to the problem, the state 
vector is extrapolated forward one revolution and then the in-plane central angle 
0 JP between the state vector and the projection of the landing site into the orbital 
plane is computed. Analysis has shown that for some selections of orbital inclina- 
tion and landing site, the correction to the assumption of one revolution may be as 
large as 29°. A lower limit on O^p of -30° was chosen, thus allowing a small 

margin from the empirically determined limit of -29° . The upper limit on 9 

o 

is +330 Large positive values for 9jp only occur in situations where no solu- 
tion existed on the previous revolution. 

To determine Sjp , the following equation is used, 
e,p - cos-' [unit (r„) . sign [(r„xi'j^p) - _i^ ] 
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where 


-LSP 

Ih = 


vehicle position vector 

unit vector in the direction of the 

landing site projection 

unit angular momentum vector 


This places 0jp between -180° and +180° and therefore an additional test, shown 
in Figure 4b, is made to force 0jp between -30° and +330° . 


In order to make the first entry into step one compatible with subsequent 
entries, the state vector is initially extrapolated forward beyond the earliest desired 
time-of-landing tg-pL, by one -twelfth of the orbital period, to the time tg, where 


*3 = *ETL * 

One -twelfth of the period is nearly equivalent to a central angle of 30° for typical 
(near circular) orbits, and hence makes the first entry into step one compatible 
with later entries. 

4. 3 Effect of Approximate Entry and Deprbit Times- 

of-Flight on Entry Crossrange uaicuiatiorT" 

During the first step in the solution process, an estimate of the time of 
landing is necessary to compute the inertial location of the landing site and the 
associated entry crossrange. Since the parameters of the deorbit trajectory have not 
been computed, the deorbit and entry times -of- flight are not known. To estimate 
the landing time, a constant At is used to approximately represent the differ- 
ence between the sum of the deorbit and entry times -of- flight and the time spent 
in orbit over the same total central angle. Preliminary analysis has shown that 
if an average value is selected for this time difference, the maximum error will 
be about 6 minutes. This analysis, described in Reference 7, did not include 
• ariations in entry time-of-flight for the particular entry range, but further 
analysis is ex ected to show this effect is small. 

During the first step in the solution, this error will affect the calculation 
of the inertial landing site vector and subsequent entry crossrange computation. 

This effect on the crossrange estimate will be largest for deorbit from a polar 
orbit, and result in a maximum error of less than 90 n.mi. To insure that potenti- 
ally acceptable solutions are not rejected due to errors in the initial crossrange 
estimate, the rou^ check on crossrange during the first step uses a test criterion 
90 n.mi. larger than the acceptable crossrange input to the routine. In step two. 
after the time-of-landing has been refined, the crossrange is recomputed and a 
precision check is made. Thus a few cases which pass the first test will be re- 
jected later. 
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4.4 Velocity Change Minimization Method 

Step two of the routine includes an iterative search to determine the location 
of the impulsive maneuver which minimizes the velocity change Av. As shown in 
Figure 4d.this iteration uses 9^ , the central angle traversed between the impul> 
sive maneuver and entry inte’'face, as the inoependent variable. A very simple 
halving step iterator is used to search for the minimum. Although this does not 
converge quickly, it is safe and reliable. The more efficient technique of using a 
slope iteration was not selected because analysis has shown that inflection points 
exist in the relationship of Av and 9^. Thes'* Inflection points would greatly 
complicate any iteration designed to determine the minimum by driving the slope 
to zero. 

4. 5 Required Velocity Equations 

The equations used in the previously described iterative loop to determine 
the required velocity can be found in Reference 2. These equations, shown in 
Figure 4d of the detailed flow diagram, use the initial vehicle position r ^ , the 
entry interface position ££j< and the desired entry angle as follows. First 
the tangent of liie initial (post- impulse) flight path angle is computed by 

tanv^ = (1 - rj^/rj,j) cot (9p/2) - rp/rgjtan(yj,j) 


where 9^ is the central angle between r ^ and £j,j and also the independent 
variable in the search. The semilatus rectum p ^ of the deorbit trajectory can 
then be determined from 

(ro/rEi)^ Py -a + tanyj2) 


The parameter p y . the secant squared of the desired entry angle, is computed 
once during initialization of the routine. 


The horizontal and radial components of the required velocity are then 
obtained from 


''hd ■ 

'' RD “ HD ^1 


The required velocity is then formed and differenced with the premaneuver velocity 
to obtain the impulsive Av. 


V 

- req 
Av * 


-req " -D 
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4. 6 Entry Time-of-Flight Computation (TBD) 

In Figfure 4e of the detailed flow diagram, the time -of- flight At 23 from 
entry interface to landing is shown as a function of entry velocity, flight path angle, 
and range. Functionalization of this time-of-flight will be included later when entry 
guidance analysis is complete. 

4. 7 In-PIane Effect of Approximate Deorbit 

and Entry Times -of -Flight 

As discus^ 3 ed in subsection 4. 3, the first estimate of the inertial location of 
the landing site is dependent upon an estimate of the time-of-landing. A constant 
time difference At used to estimate the landing time, may be in error by as 
much as 6 minutes. This led to a significant error in the crossrange estimate for 
a high inclination orbit. For orbits of lower inclination, where the movement of 
the landing site can be nearly parallel to the orbital plane, this same error can 
affect the definition of the entry interface location used in the Av minimization 
iteration. 

The entry interface location, computed early in step two, is based upon the 
projection of the landing site vector into the orbital plane and the desired entry 
range. After the minimization process is complete, the deorbit and entry times- 
of-flight can be accurately calculated. As shown in Figure 4e, another calculation 
of the inertial landing site position is made, thus removing the error due to the 
At approximation. To maintain the desired entry range input to the routine, 
the entry interface position is recalculated. This new position will be, at most, 

1. 5° (equivalent to 6 minutes of earth rotation) from the entry interface usc d in the 
A V minimization. To maintain the geometry of the deorbit phase, the lime of the de- 
orbit maneuver is adjusted accordingly so that the central angle from deorbit to 
entry interface is preserved. This adjustment in deorbit time fit^^ is computed 
from the following equation 



where i'^j is a unit vector in the direction of the entry interface position used 
during minimization, i is the new value, is a unit angular momentum vec- 
tor, and t/2it is the inverse of the mean orbital rate. The cross product of the 
unit vectors is nearly equivalent to the angle between them, and the dot product 
gives the proper sign. The mean orbital rate is used to calculate the deorbit time 
adjustment from the angular adjustment. Following this adjustment to the impulsive 
deorbit time, the ignition time for the maneuver is biased from the impulsive time 
by one -half the expected length of the maneuver, thus centering the finite thrust 
maneuver about the impulsive maneuver. 
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4. 8 Co mpensation for Oblatenesa and Finite 

Ma neuver Length 

Step three of the solution process contains calculations which account for 
the finite length of the thrusting maneuver on the required velocity change, and 
compensate for the effects of the J g gravitational perturbation on the deorbit 
trajectory. The Precision Required Velocity Deb*rmination Routine is used to 
accomplish these objectives, and the reader should refer to Reference 1 for a 
description of the technique. That routine, however, is designed to maintain the 
terminal (entry interface) time-of-arrival, and this can cause changes in the entry 
angle. Preliminary analysis, described in Reference 7, has shown that the nomi- 
nal entry flight path angle error resulting from the oblateness and finite maneuver 
length is about 0. 2^, but can be as large as 0. 6° in extreme cases. Therefore, 
to preserve the desired entry angle, the time-of-arrival at entry interface is ad- 
justed slightly. Delaying tne time-of-arrival tends to lofi the trajectory and thus 
increase the entry angle. An earlier time-of-arrival will depress the trajectory 
and result in a shallower flight path angle. 

To determine the time-of-arrival adjustment, the approximate senr ivity 
of changes in timo-of-flight to changes in entry angle is used. Analysis has shown 
that this sensitivity varies by a factor of about 13. depending on the characteristics 
of the pre-maneuver trajectory. However, the sensitivity uivided by the Ueorbit 
time-of-flight varies by a factor of less than 3. This variation is sufficiently small 
such that a constant can be used as the sensitivity coefficient for all cases. 

To reduce the computations required to constrain entry angle, both here 
and in the Powered Flight Guidance Routines* , the secant squared of tne entry angle 
p is used rather than the actual angle. In p>articular, no inverse trigonometric 
function evaluations are required. 

The sequence of calculations designed to reduce the entry angle error are 
shown in Figures 4f and 4g. First the error Ap,^ in the secant squared of the entiy 
flight path angle is computed from the following equation; 

, 1 

V ” r T 9 ~ Py 

1 - j^unit ( r 2 ) • unit ( v'^ ) p 

where r g and v'g are the terminal position and velocity determined by the 
Precision Required Velocity Determination Routine and py is the desired value. 

If the error is too large, the entry interface time-of-arrival 1 2 is adjusted as 
follows: 

‘2 ■ ‘2-S^'l2 

The Powered Flight Guidance Routines, aescribed in Reference 3, 
use the same basic technique described here to maintain entry 
angle in the event of off-nominal thrusting conditions. 
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where is the sensitivity coefficient described earlier and is the time -of- 

flight from deorbit to entry interface. After adjusting the time -of -arrival, the 
Precision Required Velocity Determination Routine is recalled with the adjusted 
time-of-arrival and the results are checked. 

4. 9 Offset Entry Angle 

In the process of computing a required velocity, the Precision Required 
Velocity Determination Routine computes an offset target for use during the 
powered fli^t. For the deorbit maneuver, the powered flight guidance also 
requires an offset entry angle. This offset entry angle, actually the secant squared 
of the angle, is computed from the following equaticxi 


Pppy 



unit ( r p p Mj ) 


1 


where r pp..j. is the offset target for the powered flight guidance and is 

the associated velocity. 
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5. DETAILED FLOW DIAGRAM 


This section contains detailed flow diagrams of the Deorbit Targeting Rou- 
tine. 

Each input and output variable in the routine and subroutine call statements 
can be followed by a symbol in brackets. This symbol identifies the notation for 
the corresponding variable in th<> detailed description and flow diagrams of the 
called routine. U^en identical notation is used, the bracketed symbol is omitted. 
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Figure 4a. Detailed Flow Diagram 
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Figure 4c. Detailed Flow Diagram 
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STEP 


Figure 4e. Detailed Flow Diagram 
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Figure 4f. Detailed Flow Diagram 
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Figure 4g. Detailed Flow Diagram 
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6. SUPPLEMENTARY INFORMATION 

As mentioned in Subsection 4. 6, the equation necessary to determine entry 
time-of-flight is not specified. This must necessarily be postponed until further 
entry guidance design and analysis makes functionalization of this parameter pos- 
sible. 

The effect of different azimuths on the entry phase has not be considered in 
this design. Since the crossrange capability of the vehicle is dependent on azimuth, 
and thus may be larger in one direction than the other, this si ould be considered in the 
acceptable crossrange criterion. In addition, varying entry azimuths also effect 
the required inertial flight path ^ngle, and hence it may be more desirable to con- 
strain the relative flight pam angle at entry interface. 

The deorbit targetl.:g technique presented here requires long term state 
vector extrapolation to select a deorbit opportunity several revolutions later. The 
accuracy of th.3 ext opolation is dependent upon accurate knowledge of the current 
state vector. Any errors due to imperfect navigation will tend to be magnified by 
the long term extrapolation. Fortunately, the out-of-plane component of error 
tends to oscillate, and thus can be expected to remain below about 1 n.mi. Con- 
sequently, the prediction of crossrange for a deorbit opportunity several revolutions 
later is not significantly affected by the expected state vector error. Mission 
planning functions, which are sensitive to crossrange requirements, can be carried 
out without being significantly affected by the navigation error. 

The in-plane component of state vector error does increase during long term 
extrapolation, at a rate of about 1 n.mi. per revolution. Therefore the targeting 
should be repeated during the revolution prior to the deorbit maneuver, taking ad- 
vantage of more precise knowledge of the vehicle state vector. 

The deorbit targeting equations presented here are designed to target a 
single minin;um fuel maneuver which satisfies entry interface and landing site con- 
straints. Since deorbit is the final major' maneuver, it may be logical to use 
all the remaining fuel to either effect a faster return or to reduce the entry cross - 
range requirement. This could be accomplished by (1) thrusting out-of-plane or 
(2) using multiple maneuvers to place the vehicle in a phasing orbit prior to landing. 

The first alternative could easily be added to the design presented here. The second 
alternative, using phasing orbits to adjust the timing and hence inertial landing site loca- 
tion. would require an additional logical structure to select the proper phasing 
orbit; however, targeting for the final deorbit maneu"er could still be accomplished 
with the equations presented here. 
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1. INTRODUCTION 

The fast ‘time integration (FTI) guidance approach was 
originally developed for the MSC straight-wing orbiter reentry 
vehicle as a candidate guidance logic. The background study 
leading to the development of the original FTI subroutine was 
documented in reference 3. The details of the guidance 
algorithm as it applied to the single -control variable 
straight-wing orbiter appear in references 4 and 5. Exten- 
sions and significant modifications to the original routine 
have followed the decision to apply the method to the high- 
crossrange vehicle designs. These modifications include the 
addition of a second control variable with the resulting 
capability of controlling not only range but also any other 
parameter which depends only on the trajectory state varia- 
bles. This report is intended to provide current documenta- 
tion of the prediction, the sensitivity, and the guidance 
equations. Also included are a functional flow chart of the 
logic and a listing of the current FTI subroutine. Prelimi- 
nary simulation results appear in reference 7 which is a 
collection of viewgraphs used for an oral presentation of 
the method. Additional simulation results will appear in 
future memoranda. A discussion of the relationship of the 
FTI guidance approach to onboard optimization and to 
parametric guidance methods is provided in reference 6. 

The prediction and sensitivity equations ure presented in 
Section 2. Section 3 discusses the control equations. 

The Appendix providep a fanctional flow chart of the guidance cent- 
putational procedure pnd a listing of the current FTI subroutines. 

(The listing is supplied for definition of factors used, but not 
otherwise defined in the text.) 
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SYMBOL DEFINITION OR DESCRIPTION 

a Break fi'-;:iuency in damping filter. 

b Control variable used to shift the angle of 

:;ttack profile. 

c Weighting factor used in minimum effort control 

logic to scale relative magnitudes of Aa and 
A(^ . 

Drag coefficient. 


d 

D 

DET( ) 

So 

g*(nT) 


h 

i 

K 


K 


0 


Lift coefficient. 

clT 

Constant in damping logic equal to e 
Aerodynamic drag. 

Determinant of ( ) . 

Surface gravitational constant used in 
predictions . 

Impulse response of the damping- loop transfer 
function. 

Altitude . 

Instantaneous orbit inclination. 

Constant used to define the range (heating 
rate) error when minimum effort heating rate 
(range) control is desired. See Section 3. 

Gain in altitude rate damping logic. 


I Great circle crossrange distance. 

L/D Lift to drag ratio, 

m Vehicle mass. 
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SYMBOL 





Q 


R 


r 


E 


DEFINITION OR DESCRIPTION 

Guidance sensitivity matrix defined by equation 3.1. 

Element from row i and column j of the P 
matrix. 

I 

Reference convective heating rate at the stag- 
nation point of a unit sphere; 

Predicted downrange distance at the cutoff 
condition; i.e., R ■ 

Earth radius. 


S 


Reference area. 


s 

T 

u 


V 


V 


s 


Laplace transform variable. 

Guidance cycle time. 

Nondimens ional horizontal velocity; 


Total vehicle velocity relative to a rotating 
earth and atmosphere. 

Circular satellite velocity at 400,000 feet 
altitude. 


X Great circle downrange distance. 

a Angle of attack. 

Specified angle of attack profile used in 
guidance predictions and sensitivity calculations. 
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SYMBOL 

B 


Y 

yH 


6 (p 


<t> 


P 


0 


y 



n 


DEFINITION OR DESCRIPTION 

Reciprocal of atmospheric scale height in 
exponential atmosphere used in guidance 
predictions . 

Flight path angle. 

Sampled flight path angle. 

Damping correction to bank angle. 

Bank angle; vehicle rotation about the velocity 
vector. 

Assumed sea-level density in exponential 
atmosphere model used in the guidance 
predictions . 

Instantaneous latitude. 

Heading with respect to the original downrange 
direction. 

Eartii spin rate. 

Nondimensional ( ) . 
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2. PREDICTION AND SENSITIVITY EQUATIONS 

The equations of motion which form the basis of the 

guidance predictions and the sensitivity calculations are the 

following. 






- (r^ ♦ h)a)g cos^ i 


( 2 . 1 ) 



+ 2o)_ cos i 
£ 


h * V sin Y 

X = V cos Y cos *1' 

i ^ pQ® ^^'vSCpCL/D sin 

i V sin ip cos Y 


V “ (Tg f h) 

( 2 . 2 ) 

(2.3) 

(2.4) 

2b)g sin i sin y (2.5) 

^ ( 2 . 6 ) 
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Note that these equations include rotating earth effects but 
are based on an exponential atmosphere. Computationally the 
atmosphere consists of two separate exponential fits with a 
conti lUu 5s junction at an altitude of 140,000 feet. 

T}(e equations actually integrated by the guidance are exactly 
equivalent to those above but have been transformed sub- 
stantially to permit rapid and efficient integration. The 
algebra of these transformations is tedious but straight- 
forward and therefore is not included here. However, the 
transformation procedure is itemized as follows: 

1. Convert the dimensional equations to non-dimensional equa- 
tions using the vehicle mass, the circular satellite 
velocity at 400,000 feet altitude and the circular orbit 
radius at 400,000 feet altitude as non-dimensionalizing 
factors . 


2. Change the independent variable from time to non- 

dimensional horizontal velocity u by dividing equatiojis 
(2.1) through (2.6) by equation (2.1) and multiplying 


the resulting 

equations 

by 


A 

dv 



1 

du 

cos yT: 

L - 

0 tan Y 

i . e . , 

ILi . 

du 

Zft ^ J dv 
dv/dt du 


(2.7) 


( 2 . 8 ) 
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3. Factor the resulting equations extensively so that quan- 
tities appearing frequently in the equations are evaluated 
only once for each derivative evaluation. Although this 
factorization renders the equations almost unrecognizable 
in the guidance routine, it is perhaps the most important 
coding detail to permit rapid integration. 

In the guidance subroutine, the equations of motion appear 


as follows: 

D(l) = (T56) [(T27) (LOD) - (T5) (T26) + 2(T2)] (2.9) 

D(2) = (T56)((VN)(T10)] (2.10) 

D(3) - (T56) [(T25) (T5) (T99)] (2.11) 

D(4) - (T56) [((T27)/(T7)) SIN(PPHI) (LOD) 

+ 2 (OMEGA) SIN (ANGLEI) SIN (XMU)] (2.12) 

D(5) - (T56)[(VN)(T8)(T5)] . (2.13) 


The factored "T" quantities are defined by the code itself 
which appea.' in l^he Appendix of this report. 

The perturbation equations integrated by the guidance result 
from taking partial derivatives of equations (2.1) to 
(2.4) with respect to the quantities (L/D cos 4) and . 
These quantities are used instead of a rnd directly t ) 
avoid algebraic sign reversals in the sensitivities as the 
bank angle goes through zero degrees. The sensitivities to 
(L/D cos 4) and .o are related to the sensitivities to 

a and by simple chain rule differentiation; i.e., 
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(2.15) 


In the guidance subroutine, the perturbation equations 
for (L/D cos <t>) take the form: 


D(6) = 

(A11)(Y(6)) 

♦ (A12)(Y(7)) 

+ PFLODl 

(2.16) 

D(7) 

(A21)(Y(6)) 

+ (A22)(Y(7)) 

+ PFLOD2 

(2.17) 

D(8) = 

(A31)(Y(6)) 

♦ (A32KYC7)) 

+ PFLOD3 . 

(2.18) 


The perturbation equations for appear as 

D(9) = (A11)(Y(9)) + (A12)(Y(10)) + PFCDl (2.19) 

D(10) - (A21)(Y(9)) + (A22)(Y(10)) ♦ PFCD2 (2.20) 

Dfll) = (A31)(Y(9)) + (A32)(Y(10)) + PFCD3 . (2,21) 

The values of the aerodynamic coefficients C and C and 

L D 

their derivatives with respect to a (C and C ) are 

L D 

a a 

included in tae prediction and sensitivity integrations using 
cubic spline approximations (reference 2). Interpolation using 
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splines is very nearly as efficient as using linear inter- 
polation and provides accurate derivative information at no 
cost. It was found that very simple polynomial approximations 
for angle of attack dependence were too inaccurate to produce 
a net gain in guidance efficiency. Mach number effects are 
included using linear interpolation. 

The crossrange equations (2.12) and (2.13) are integrated by 
the guidance with the bank angle profile multiplied by -1 to 
predict the crossrange capability that remains if the bank 
angle were reversed at the current time. The crossrange 
logic commands reversals when the reversed-bank crossrange 
equations predict termination inside a shrinking interval 
about the target. As a result of using in equations 

(2.12) and (2.18), the value of from the crossrange inte- 
gration is not fed into the downrange equation. Instead, 
the value of ip used in equation (2.11) is obtained as 
follows : 


i (L/D sin (2.22) 

V » 'D/m (2.23) 

Assuming that L/D sin is constant over each increment of 
velocity corresponding tc a single integration step 

♦i+i ■ 
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This simple integration provides a value of sufficient 
accuracy to use in the downrange equation (2.11). 

The value of the latitude used in equation (2.5) is calcu- 
lated as in Figure 1. 

A predictor-corrector integration method is used in the 
currant FTI guidance routine to perform the numerical inte- 
gration. Four steps of a fourth-order Runge-Kutta method 
are used to start the integration. The predictor is the 
fourth-order Adams Bashforth formula and the corrector is 
the fourth-order Adams -Moulton formula. Stepsizes are fixed 
at 50 feet per second. The integration code appears as an 
integral part of the guidance routine to avoid calls to 
external and overgeneralized integration subroutines. 

To date, no attempt has been made to strip down either the 
equations of motion or the integration. The synthesis 
philosophy has been to build from formulae that are unques- 
tionably accurate enough and to do so with as much efficiency 
as possible. Complete trajectory integrations are currently 
being performed in approximately 2 seconds of UNIVAC 1108 
execution time. It is probable that an investigation of the 
following items would lead to a substantial reduction in 
computation time: 

1. Simplification of the prediction and/or the sensitivity 
equations based on an order of magnitude examination of 
the terms. 

2. Increase of the integration stepsize. 
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3. Alteration of the order of the integration formulae. 

For example, with current stepsizes, fourth-order formulae 
are probably not necessary. Elimination of the corrector 
alone would almost halve the integration time. Other 
investigators have recently recommended the use of 
higher-order integrators with very large stepsizes to 
improve efficiency. Ihese suggestions should be investi- 
gated. It is anticipated that by using state-of-the-art 
integration methods the integration efficiency could be 
improved by at least a factor of two. 
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sin A . “"C90°-y 

sin 90° sin(90°-Ug) 


CCS a 

o 

cos y 


sin A 


cos c ■ cos X cos 1 


sin A 


sin I 
sin c 


A” - (180°-A) ♦ A* 


cos (90°-y) 


sin y ■ cos c cos(9C®-y^) 

sin c sin(90°“y )cos A" 

o 

cos c sin y -f sin c cos y 
o o 


Figure 1. - Computation of latitude at point : 

(2 of 


( 
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3. GUIDANCE EQUATIONS 

As currently coded, the FTI subroutine issues angle of 
attack and bank angle commands based on the reference unit 
sphere convective heating rate at a future point and the 
downrange distance traveled when the vehicle reaches a speci- 
fied cutoff altitude. It is apparent f. hi reference 1 that 
approximate temperature control can be maintained by speci- 
fying an appropriate reference sphere heating rate profile 
for the nose panel of the orbiter. The guidance algorithm 
attempts to make the predicted reference heating rate at a 
future point on the trajectory equal to a desired value 
obtained from a specified profile. The prediction point is 
determined by subtracting a fixed velocity decrement from 
the current velocity. For example, suppose the specified 
velocity decrement is 400 feet per second. At a guidance 
computation which begins at a velocity of 23,000 feet per 
second, the heating rate at 22,600 feet per second is 
predicted. Also calculated are the sensitivities of the 
heating rate at 22,600 feet per second to the current con- 
trol values. Similarly, when the vehicle is traveling at 
22,400 feet per second the predicted heating rate at 
22,000 feet per second is being controlled. Simultaneously, 
the predicted range at the cutoff altitude is being calcu- 
lated along with its sensitivities to the current control 
values. 

The guidance computes the elements of the sensitivity matrix 
P , where 
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P 



dR 

^1 



(3.1) 


The sensitivities of R are evaluated at the cutoff condi- 
tion and the sensitivities of Q are evalucted at the 
prediction time appropriate for the current guidance call. 
Linear analysis leads to 



LAQJ 



(3.2) 


where AR and AQ are determined by differencing the predicted 
and desired values of R and Q , respectively. The con- 
trol commands are obtained by solving this simple linear 
system simultaneously; i.e.. 


Aa 


ar’ 


■ P 

• 





(3.3) 


To first order, the above control law finds the control 
pair which eliminates the range error and the heating rate 
error simultaneously. 


It is easily seen from equations (2.14) and (2. 15) that P 
is singular when ^ - 0 or 4 ■ 90® . In the subroutine, 
the magnitude of the determinant of P is calculated and 
if necessary limited to a value greater than a preset 
minimum. In addition, if the determinant of P gets small 
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causing large control changes to result from 3.3, these 
changes are limited so that |Aa| <3® and |A4)| £ 10® . 
This allows the controller to make the next prediction with 
a value of 9 for which P is non-singular. 


Since it is necessary to provide for upper and lower limits 
on both angle of attack and bank angle, special control 
saturation logic is used. Each forward integration is per- 
formed with control limiting logic. For example, if the 
controller has shifted the angle of attack profile so that 
in some segment of the trajectory a exceeds its maximum 
allowable value, the trajectory prediction and the sensi- 
tivity calculation are made using the maximum allowable 
value in that segment of the trajectory. Thus, all pre- 
dictions and sensitivity calculations are based on perturbing 
the control only in regions where perturbations can be 
tolerated. This represents another major advantage of 
onboard sensitivity calculation. Other perturbation guidance 
schemes, such as the adjoint variable methods, cannot easily 
adapt to control saturation and therefore provide inaccurate 
sensitivities in those regions. If a control is saturated 
and if the change given by (3.3) does not cause the control 
to leave the control boundary, the FTI control law recomputes 
the other control variable change based on the proper 
single-control variable sensitivity. An example ,.ould help 
clarify this simple saturation procedure. Suppose 

a = a and Aa computed from (3.3) is positive. The corn- 
max * 

manded change in is then recomputed as either 


A^ 



A0 





or 


AR 
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depending on the control segment in \»hich the guidance is 
operating. Similarly, suppose 4> is saturated and A(^ 
computed from (3.3) does not bring off its limit. The 
Aa command issued by the guidance is then recomputed as 


Aa 


. . 

3Q/8a 



or 


Aa 


AR AR 

(3.5) 


Again the equation chosen depends on the current control 
segment. Note that the recomputation of the commands during 
regions of control saturation requires only an additional 
division since the single-control gains are elements cJ the 
sensitivity matrix P which is already computed. 

The control segments referred to in the previous paragraph 
are defined as follows: 

Segment 1 — From guidance startup at .05 g's deceleration to 
a fixed velocity point. The controls in this 
segment are calculated to eliminate the ^ 
error with minimum control effort. 

Segment 2 - From the end of Segment 1 to the point on the 
trajectory where the predicted maximum heating 
rate first becomes less than a threshold value. 

The controls in this segment are calculated to 
eliminate the Q error and the R error simul- 
taneously. If a control saturates in this seg- 
ment, range control is dropped temporarily and 
the ^ error is used to compute the other control. 
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Segment 3 - From the and of Segment 2 to cutoff altitude. 

The controls in this segment ere calculated to 
eliminate the R error with minimum control 
effort . 

The minimum effort criterion is used in regions where the 
control of one parameter takes high priority over the other. 
In Segment 1, the objective of getting onto a desirable 
temperature profile takes priority over range considerations. 
Therefore, the quantity AR is computed as a scalar constant 
times the required AQ and the constant is determined to 
minimize the control change. For example, let the measure 
of control change be given by 

J = Aa^ + cA(f>^ (3.6) 

where c is an arbitrary weighting constant. Assuming that 
range control is to be temporarily ignored, let 

AR - KAQ (3.7) 

Then 

Tao 

u. 

or 

“ crrTt'P2i 


Jkaq 
AQ . 


(3.8) 
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where DET P is the determinant of P . 


Substituting these expressions for La and L<^ into J and 
setting 3J/3K » ^ results in 


K 


P P 
12 22 


cP P 


12 


cP 


11 


(3,11) 


It is easily shown that this value of K represents a rela- 
tive minimum of J . A similar result follows if AQ is 
specified as a constant times AR. 

It has been found in simulation results that the minimum 
effort criterion to specify the "uncontrolled" parameter 
results in excellent targeting in Segment 3 where heating 
problems are no longer of consequence. Although extensive 
simulation verification has not yet resolved the proposition, 
it appears that Segment 1 may not be necessary. The advan- 
tage of Segment 2 cannot be overemphasized. In that segment, 
a desired heating rate profile (and therefore an approximate 
panel #1 temperature profile) can be tracked while the 
guidance is simultaneously solving the range problem. 

Typically, Segment 2 covers from 5 to 10 minutes of trajectory 
time during which heating considerations are important. 

Neglect of direct range control for this length of time and 
at the high-energy end of the trajectory causes a con- 
siderable loss in targeting ability. Simultaneous control 
of range and heating rate in Segment 2 is overridden only in 
regions where the control saturates. These regions can be 
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minimized or eliminated by choosing realistic control pro- 
files for the prediction integrations. The result is that 
Segment 2 provides both path- type and terminal -type control 
and appears to increase substantially the guidance footprint 
areas over that obtained with a temperature-only control 
segment followed by a range-only control segment. 

One of the principal advantages of integrated predictions 
lies in the fact that completely arbitrary control profiles 
can be assumed in the prediction. Currently the program 
uses a constant bank angle profile and a variable angle of 
attack profile of the form 

a ■ b + “p(u) (3.12) 


where 


a (u) = 50° 

u ^ 0 . 8 


P 



= 50° - 

20° sin ^ "2 (0.8 - u)^ 

oo 

• 

o 

V| 

V| 

• 

o 

» 30° 

u < 0.2 

(3.13) 


The guidance adjusts the constant b to increase or decrease 
the entire profile by the amount necessary to null the gui- 
dance errors. The a profile choice was made to allow the 
guidance to stabilize and still permit the pitchover maneuver 
necessary to leave the vehicle at a small angle of attack at 
guidance cutoff. The best choice of prediction profiles has 
not been studied to date but it is anticipated that such a 
study could produce improvements such as enlarged target 
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areas, more nearly optimal trajectories, and reduced sensi- 
tivities to navigation errors. The possibility of using 
several parameters to specify the prediction control profiles 
and of using the guidance to adjust these parameters in flight 
is discussed briefly in references 3 and 6. 

In the FTI guidance routine, both the prediction bank angle 
and the commanded bank angle are modified, if necessary, to 
provide trajectory damping. It has been found that at lower 
angles of attack the delta-wing configurations have a tendency 
to show altitude oscillations that are undesirable. These 
oscillations can be overcome by using a simple flight path 
angle feedback. Since low frequency flight path angle 
changes mi’st be permitted, a simple high-pass filter in the 
feedback loop is suggested. The bank angle command is used 
for damping. The feedback loop has the structure 
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The transfer function from the sampled flight piith angle 
Y* to the bank angle correction 6^ must include the dynamics 
of the hold and is therefore given by 


G(Z) 



(3.14) 



(3.15) 


where d » e*^ ; T is the sample period, and Z represents 
the Z- transform. The inverse of G(Z) is 


g(nT) 


Kq for n ■ 0 

d"”^(d - 1) for n > 1 

VO — 


(3.16) 


Using the convolution summation, the correcticn 6^(nT) is 
obtained as 


f 


{♦(nT) • ^ Y*((i" - n)T]g*(nT) (3.17) 

n»0 

The choice of an appropriate value for the filter break 
frequency, as based on the frequency of the undesirable alti- 
tude dynamics, results in a rapidly converging series. The 
guidance uses four terms of the series with approximately a 
3% error due to truncation; i.e., 

64(t) - Y*(t)KQ ♦ Y^(t - At)g*(t - 2 At) 

Y*(t - 2 At)g*(t - At) +Y *(t - 3 At)g*(t) 

(3.18) 

> 
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This simple damping logic using three backpoints from pre- 
vious guidance passes has proven very effective in smoothing 
the altitude-velocity profiles. Since the damping is incor- 
porated into the prediction and sensitivity calculations as 
well as the commands, no degradation of the targeting accu- 
racy is produced. 
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IP! VW «6T« . *65 ) t»oO E = i 

TTi VN ,LC. .AMT I»0PE = 2 
CUTOPF s 100008 


TPlLTIIlJ 


BC»* _ 

32.17% • {R£n/Rn)*«2 • RO / V0«*2 
TS»tX • vO / RO 




p uai = Tinex 

PTo%) s MMO * RO / 6076.007 
P(65J s XB • RO / 6076.007 
RJ66) = VN 

F(6TI s ALTUCC • 3.28n8 


FI67I . 07 . ATMOSH) 00 TO 
BETA s KTA8 




30 COHTlwUF _ 

P(M» « ASvAUVAtt) • 3.28003 

FfOOl a ABVM.(OOntOT) • 3,20083 / 32.17% 


a rill • BT.20S795 
FITII a YUl • 07,208705 
PITH a rim • BT-.WstoS 
FC7SI a F(0) • 57.20R795 

PlfoVi'llAS ' ~~ 

F175J a PHTSOH 

IF (t|6<»» ,LY. ~ir;8yrinFT<rT5o" 

IP tT(2)*00 .LT. CUTOFF) 00 TO rOO 
-STPtBT r 0 

C Pooor IS TMC VALUE OF THE Frx. hEaTtHO PATE ON thf PBFVtOIK CALj 
IPF>OO eT ,0t. t*S3D ♦is.) ALPnT?T = ALPHTv ♦ 9./ft7,3 

IF! POnOT .LT. P CSOO^S. ) *LP« ir< = ALPv'IW - 9./sT.^ 

iPl ALPFIN .LT. 20./R7.3) »LP"tH =20^/57.^ 

IF|AtPviN_^T._%0^57._?« AtPNt- =%P./57.3 

POOOTaO. 


C* 

c 




RUN6E KUtTA L^'OP 


SO J s 1 

c thp POTrff5E'onowvCT rs TO >nnn>/mE * pl*g to ppohibit tmp c*lcuL< 

C OF SOME OIIANTITIFS IN T^E OEPIVATtVE SECTION ON ALL EwALU*TIONS 
C EXCEPT TMEFlfiST 
lORVCT 7 0 
^CrTO*BEPlV^ 

•f 00 JKPSlfNEQS 


w 

IFI STP CNT .or. 1) en TO hj „ 

CODOt 7 17600, * SrwRT( T3 *«h6n 6/TPM. lO) ♦ VN**3.1S 
«PITC(6»S000I COOOT 


5000 FOPPAtllH ,13PrUPPCNT OOOTs,E20.6) 

_ FttOI g COPOT 

67 COITITmUE 

SSOM = Tl9/AnS(Tl9l 

»»g5o,/vo*sss*r 

X g X ♦ O.S»H 

70 JKXgltNEoS 

70 TU’XP) g Y(JX«) ♦ 0,5 *M*KHjKW> 


J 7 2 

GO TO DEPIV _ 

66~W 65 jkW7i»NEoT" 

F2 (JXW»g O (vlKX) , 

is V(JKM) 7 Y(JKM) - 0,5*H«Ki ( JKXV ♦ 0»S«H*K2IJKX) 

BDTinsiRfv 

90 00 98 JX«7l»NC0S 

K9<JKP)7 OIJKW) 

05 TJJMI) r Y(JKP» - 0.5*H*K2( jKW) ♦ H*K3(JFW) 

X = X ♦ 0,5*H 

J 7 9 


100 00 Its JKMt«Nc05 


K«|JKW)7 nSJXW) 

10s TtJKWI g YtOKW) - H»X3(JXP) ♦ H < . 16666666* ( 1 C JHW )■»?, aX? ( jKw 
i K3«jKwn«9iJKWn 

VMB X/COSIOANmA) 


V X VN*V0 

C APHXMO IS the HCAOfMo USED TO *»00 |Fy THE DOPNPANOE PPEOtCTIQN 


V9S0 .)/V) 

IFIABSI APPX Hr»S7 , 3 ) . OT . AS . ) AFw XHOs^S . /57 . 3«S6N t APpXHn ) 

7 cost appkhBt 

1 F( y( 2)*oQ .’^T . 350000.) GO To A60 

lf( IFC66) - VN) ,6T, f^TV/Vo) 60 TO llO 

e calculation of sensitivities fop path contpol. 

STNSl X Y(6) 

SENS2 7 YI7) 

SEH53 ~ — — 

SENS97 YI9) 




till 


I 




SNS002: S£NS30*Stf «;«: 

SV«0 = OOOT 

UO COUtlNUC 

P63=pe2 __ 

*^«4=P«l 

PeisG*»*m* _ 

ALT = HN • PO 

IP t ALT .LT, cutoff > GO TO ENOpT 
IF^LT ,6T.“ATW5H) go to Jl5 

MTA = eCTA2 

s RHON02 

lis COHTiNUe _ _ 

TN = TN ♦ H/T19 

PACT s alt 

— stfctt s STWJHT ♦ T 

IF |T19 .»E. PTI9> eo TO 120 


= T19 
ACCVN 8 VN 

"lio OOTT I~1760(1. • 5^ fT 3 *ikHON 07 TPHONO» • V9*«"S, 

IF I OOOT .LE. POOOT) GO TO l25 _ 

WOOT 3 00^ 

OOTVN 3 VN 



IF tSTPC NT .CO, GO Tp 1«, Q 

00 IGO JKWrIfNireS 

_ 9 kpthjk w> = Y(Jfw) _ 

WniV9<jNW>= BkDPV3(Jwh) 


loo BK0«V2<jKm 3 Kl(JKM) 

00 T O SO 

IM 00 155 JKW.'liNF05 

"APTK^F) = YUAWI 
155 BKDPVIIJKW) = kIIJki*) 
»•••*«*••••*>**•«**«*•***•. 


>♦*******•••♦*■1 


PWCOIC TOP 


C«**»»*»«»»*»«**« i «»»« *»*»**«*»» *««*»*.>*»«»»»»» a »*«»»» y * *«»»*«»*» *»*»***» 

iM J * 5 
IDIIVCT = 



9.11 2. '3 lEn; 



190 VI JK»0 = YH*T» JKV ) 

~ 6 

60 TO dCfflV 

too DO 205 jy Wsl» NFOS 

205 vijfcw) = Y 5 *VrjK¥r + w*c?*Ju75.*D( J*r»n ♦ I4?7,*wcnqv0l J«flO - 

1 i^’KOfiVl IJKW} ♦ »fl2,»Bi<DRV2t JKW) - 173,»BKnPV3ljKW> ♦ 2T>» 

2 BKDBV9< JKW) ) 

V«ir X/COS(6AMwA) 

V a * VO 

= APKXhD ♦ L0P«T6*Al06< «V»3Q.)/V1 

3) , 6T . 507 > *^*Hd*5n . /5T , :^*SBN I APpXMT ) 

799 s COS( ABaXHO) 

(YI2)*B0 .6T.~550000, ) fid TO %60 


Ca 9 aa***«««*****»w*«« 9 *a*«****««««***«*«ft*««**«* 4 i««*«««*«*«*«*«********* 

irt (riM» - VN) .er. dltV/Voj «o toTiT 

C ^**-^^^^°** ^ SCWSITIV ITTCS F OR PAt H CONT ROL. _ 

SE NS2 = Y(7 > 

S^NS 3 = Yia> - 
SCMS4S TI9I 


l«S5= TllOi 

X<JJJ - 

'OOOT = 17600, • SfiRT(T3*RH0N0/TRHON0) • 

SCNSO O a -17600 «* < 9ETA /2i I •SORT ( T3 «WH0 NC/rRHQN0 >* VN»*3 .i5» RQ 

SVdo = OOOT 


•■■wk 


SCNSoOaSCNS? 


SHtariT- 
210 CONTINUE 



alt s 

IF CT19 
PT19 = 

ACCVN = 

HN • RO 
.6E* PT19) 60 
T19 
VN 


IF (900 


C, MOOTI fiO TO 230 


OOTVN = VN 
230 CONTINUE 


U fST;: 


POlsOANNA 


C 

C C HECK FOR SIWULA TlpW Tt RNiNATIQN 

IF IALT .LT, CUTOFF) fiO TO ENOFT 
TT C»LT ;ft, r.o TO ?90 

_ ttfcTA = BET_A2 

IWONO = RW)N 62 
290 CONTINUE 

TN = tN ♦ W/T19 


C UPDATE STORAGE 

~t 

DO 250 JR«S1*N^ 
PALT r *LT 



WORVi»( JKWI : AK0RV3<jKiit) 
BH0IW51JKW) 3 
dKMVaiJKW) S HKPRVKJKV.) 
H*0 WV1< JKW> = nUPP- il( jKt- ) 
2S0 WPTijJkKIi 3 YtjKji) 
ST^C’*T 3 STPCNT ♦ 1 
60 TO 160 


imterpcl*tion fob cutoff 


^00 PS Ac 3 SCUTOfF - i ' » / jalt* - palt» 

TN S TN ♦ FBAC*H/ii9 

OO~3l0^.iii.HF05 “ ~ 

S16 YIJKM) = BKPTKJKH) ♦ FBAC * (T( JKWI-SKPTl ( JKM) V 
FWftE«6»4007)TrvOb» VN 

MOT FONMATUH rSHSveo-,E2n.f,,6H*T VN^ffSO.f) 

FlOdI 3 SVQO 

60 TO 6UID£ 


ITvALUATlnN OF DEBIyATTVeS 


350 COHTINUF 

lORVCf 3 lORV'CT -s- i 

jrUBSlApHXHDl.eJ. 05. /ST. 31 PHI » -PHI»SS*^J( APBXMQi 

calculation of sound SPECP* bach HO,, ANn ACPOnvNARiC C''FFFTCIFnTS 
^ vti! » 


SNOm B TI2t • 00 

call SPLN2 < NSND . S A , Sg t »;PL5 nO , Vf WO ) _ 

PMACH 3 WN*V0/VSN0(2) 

!F1 PMACH .LT. 2. I pmaCH r 


OATH r 50./57.3 .20./57. 3*SlN( 3. 1A1S926/1 .2*< 0.0-VN) > 

IF< VW ,6T . 0,0 ) PAOATk « 50. /ST. 3 

~nn W ,Lt. ft.* ) ^OATK s 30./57.3 

PAOATK 3 PAOATK ♦ ^PlA9 

IF (PAOATK .6T# 60./57.3) PaOaTK 3 60./57.3 
IF( PAOATK .LT. 20,/«i7.3) PAOaTK 3 20, 


PmACH. PAOATK) 


ro.Tinmj!LT4: 


JaCCULATION of LATiTUPf 

5iijC3 ACbFt'^0rfriy)*BO>BFr*foS(YC^^^ 1 
APSN s ASIN( SIN|Y(S)«B0/PE)/«IN(SNLC> ) 

3 ^,lAlB-»2(i - proa * 

XMU 3 ASIN(C0S(SMLCUSIN(VmU 0)4. SlNt SYLC)*C05(VMU0)« 
1 eOSiA2P)l 


ILTFT IS A BLA6 USC® IN PmLL-«)P OVtnBlPF t.OOic 


iLdTTs r 

60 TO 365 
360 CONTINUE 
365 CONTINUE 


X) 






c 

T1 = 0.5 * RHONO • CO * SN 

~ ilWLTI i 51^- T(«l 

T2 = ON. SAM • C0S(AKSUFI> 

TZTf 5TRCN 

AH - T20 

TI5 = !»Fr*"*^no 

6M = 32 -17s * «REN/RN»**2 • RO/VO**2 


cost PHI I 


^4 s cost PSD 


no 

T12 

— T5 

TIT 

T2ST“ 

TO 


s SlW( S*WtA ) 

3 VN • VN 

S cost SAMPA ) 
w 6M / VM 


S RN / V« 

■ SINI PMI ) 


0 • T?1 

T 3 3 E«Pt-tT l 5 )«HN) 

Y23 s l, / T22 
♦ T23 


- T22 A t21 - T23 

T27 3 T1 * T3 * VN • T7 

■ Tfi = TIO * T25 

_ T20 3 T25 *T5 _ 

"TIV 3 -T1 A T3 • Tl2 - TiO • T30 

IE (TIP r€C. 0> MW3R0 .OT. 200000,1 SO TO »5C 

H66*te WN»|T2t/ »Ct>*T7) • { Cl.*TT-Cn*T10 > ♦ tT2!^-T17» ) 

— PWPt * » l.t/tLOO*A0S(S I M<Php I ) 

WBOhPs tt.S •57.3)/(.l.00*StNIPHI ) )*(^Avma ♦0.ft5*P01 ♦ 0.22S* 

I PS2 _♦ 0.07ft*P6S ) _ _ __ 

IFt ABStPHlO'iP ) ,GT, 20,/57T^ )PHfD'-Pr2o./S7,'3*'^fi» (PWTDMP) 

IFi ABStPHlOMp) ,6E, ABSv'PhI» )PHICMP=0 .«*abS(«»HT )*<;S^'«PHrnMp> 

— SFf OAOlTk •«T. ■SffVTsT.fT'i^tDMP 3 6; 

77 = cost f^l ♦ 7HI0mP*S#»«PmI> ) 

TlTs ti "•^5 p Wi”T7 “ 

369 CONTIMUF _ _ __ 

IE tIFLAB .EC. X) GO TO 370 




• 


ro 

ItOOMs 

1 

PWPBW 

Q\ 

vD 

WBSM 


t FG2 


CPOSSOANGE PRP^ICTION 


:u« 


SVLOC _s_LOP 

MSAVsPMACH 

ACS A VE = T 19 

vM " = V\' 

TRAV 3 ABStPHl) ♦ Pm 1 SSN*FI 9 > 

SPTCHL = -TRAV/C3./57,3r ♦ ACS*VE • VO / PO 
IFLAG 3 1 

■3TP IFTtvNl-VN) SpTCML) so to 380 

PPNI = Ft9l - PM1 SSN*TRAv/SKtCHL*(vNI'*VN) 
IE IS^TCHL .LE, 0.) PPHIs-PhI 
GO TO 3S0 


C PPMI IS THE ROLL A jSLE ,iSe-0 ?N tHE rRO^S o«N6P o9epTC^T''M 

r 



HPMl s -PM I 

coHjimjf. 

tatp = 1, / T19 

Ttl = T24 / T19 

TIS s -Til • T5 • t3o “ 

Tl* - -Til « «T1 • TS • t 15 • Tl? ^ Tin * r?\) 
T*0 = -Tl/Cn • T12 • T3 

TgM Tt7« 0O- T5>T26»a.»T? _ _ _ 

^l4= l,/( fTfl^t5‘.ir^'pTie*T58 ■) 

1188 TsagTse 
Tli*»39Tl 
T*ls T5 ♦ L0r*T7*Tln 

T#is -T60*(-7lO ♦ L0 D*TT*t5) -T28*VN« < T8*T8 • 

I >2t *71*VW«T 5 

TIT s TlW*Ti3*T6i V 
*11= iT 17 - T2t*T2l - T2S»*T10 
TT| = vNpflC 
T7»s TaSjja • T9 

▼•i = -f2f*Ldn*Tl5 ♦ T5/VN*t2l-T23*T23/VN-T5 
TIS X VMptS 
^il X -t28«Tib 
1^ m •T*S/T20*T5 

•lix TSl*TSl-T«8*T8'5*t»ft " 

...•l*5..TS6*J12-TP5*TS5*Ta7 _ 

AlSX'O • 

*aix _ISi4t*J6,'UT85*T7?*T«i 

aits •m«T72«TA7 

M *r »- 

*Sis tS6*T89 -tsS^TTTPT^l 


asss 1 , 

^PLOpix T56*T27 - Tp5*TS»i,(-Tf,0*Tt0) 
i^L002x -T68*T72*(-T60*Tl0) 

gfLOPSs »Te5»T73» (-Tfea> Tta> 

'▼CPlx ^r86*TC7>iOD - T15*TM* f •T61/CD*Tft 1 > 

SEgKS yT1S!»p a>j«Ti Q/CP9TIi 1 

^eosi -TS9*t73-*|-T6o/CnpTll ) 

*31. 8 A3l*T99 

Mt X AS2*T79 

ass s A 3 s »m 

PFUTOS x>rL00S*T99 

^CDS X Pr CDS«T99 

BTai s T56«T T2 T»lO(' -T9*fi?ri2.*T2 ) 

0 « 2 » X T54*» VNPTlOJ 
0l3»xTSiiT25*T5*T99 “ ~ 

OU) X T*16P(T2’/T7 *SlS(PpHX )»LOO * S.^O^riANPltN^ ANSlCI » 

1 ' ^ 

0(8)« T96«VN«TS*T5 

D(IT 8 aii4>vi6)>ai2*v(Tii ♦ aPtoui 

0«7» X a2^l*r(6)*a22*y(7J ♦ PF1.0D 2 

■ Ufl) X *3i*T>l6l<^A32*vTn i~rrt5B5 

OCtI X ail«f(9l-»*12«y( 10) T PPCOI 

t»tlO)x A21«y<9U*22*y< 10) ♦ Pt*’r02 

0(11)8 AM*T(p)T*32ey( 10) 4 P^CPS 

dO "to (oOtsb.^Ofido, 170,200) , 


I 



9.12-71 


COMPUTE CHIFT 


<L/0)COS(PMI) 


C**«***P«i»«««*«*«***«*«****»*V***. ••*«••**«« •«<l>***«4i**«*«««***********»* 

•00 CONTINUE 

■ (WRMft = Y|31 • ^ / f076,O^ 

IF iONRNS ,LT, 0.) 60 TO 460 _ _ _ 

CRSRN6 = YI5> • / f.076,09t^ 

_ _ VN T ABVAL CVOTI) • 3.g* q83 / VO 

“iJtfii = oes!^ - Yi3) 

OC TttO a OCSQO - SVOn 

C CVALTiSTtdN OF THE PARTlALS OF TmE *E«0 MJTM RF4PECT TO *LPM* 


CAU. APODYN«msFV» 
*WIXCL = ENt<l) 
*P»XCD g EWTt ?) 
PCUTs EHTI 7) “ 
PCD* = ENT( Hi 
PCoRAs ^ iSWrSTcDi" 


F (_7) ) 


PCOISAs iSWI^STcDi DC .A ~ -TAP‘PYCL«»PcnAr»/(ADP»C''**?1 
SCNSP s Y (8) 

fFiaMioltwT ,8t. 3T> t^s , *S8n roCYoo) 

1 OL W lOT , ^ .*ft0 76.0qT/8o>DLTPxB0. ♦ 60 76 . Q«*7/ bO»S 6N LOLII* > 
SCLCSTsi. 

PMli = SNStOl • COB (PHI) ♦ M|_0DA 4> SNSOD2 * PCOA 

Pllfif s“iNS« I i 1 -L00*l t C PH I ) ) 

P Mil « P MUi * S CtfST 

= PiiSTT ♦ S'CLCST 

PAP2J = 5ENSP ♦ COS(PHT) ♦.Pt^^OA ♦_ JfJlD ♦„ PCOA 

PAla* s scmP ♦ (-Lob ♦ siH(pMi)) 

QgT s PAPtl * PAPgg j? PAPtI ♦ PAWlg 

CS*«T s 0,2 

*KONlHPAPT^♦PAP?l*♦*PAP^^♦PAP2?^/(PAP^2♦•2♦CONBT• 

1 PAR11P^2 ) _____ 

IF( VH ,LT. ,^»53 .ANr, VN ,6T. ,272) 60 TO «10 


IFi »H)O E.EO,l) OLTPs PWOPoPP OLTOP 

,EO . 2 ) DtTOOroLTtk* ( pAlll 2 *Pil( 22 «^<^T*PAPl l*pAfl| 

t PAP21*P2»C0WST*PMtl** 2) 

■ ~¥io tOHTiwur 

C check fop SYSTEM_SiN6ULAPtTY 

WrAiStOET) .Li, l.E.15 ) beT3i,F-is 
OAI.PHA s U./OeT) ♦ (PAR 22 ^oLTeD*SCLeST-PAP 12 P 0 LTK) 
OFRI F<1,/DcT) ♦ (.PAP21^.tob*BCLCST ♦ PARU*OLtP> 
C LlMlTINt L08IC ON ALPHA ChAME 

1F< aSsTOAlHIA) ,Lt, ♦, 1/57,3 i 80 TO AlS 


OLTPs PROPOA 
.TOOrOLT**ipJ 


ST*PARllPpAflBU/( 


IF( AaS(OALPHA) ,LE. 4,1/57.S ) 80 TO A15 
PPP s 0,/57,3)/ APS1 D^PhA) 

7>1ti>^AsPPP«(}lL^* 

DPMIsPRPttDPH! 

41*^ CONTIn'JF 

LtMlTlND LOGIC ON PHI CHAnOE 

IF! AftsT6PMlTTLF7 10,1/37.3") 60 to 420 

PPPs llO,^,3)/A8SCn(^I) 

ffALPHAsPPP*OAl>tt» 

OPNIsPRP*OPHI 

-5J0^ CWTInUE 

IF( A0ATK_,8T, 58,/57,3 ,AnO. rALP HA , 6T. 
"TFf'AoItK ,Lf7*LPHIN42./S7 3 .and, OALPHA 


q.)r> o TO 42? 

'.LT, 0.)60 70 422 


IFlAPB(PHl) .•T.85,/37, Apb(B6N(OPHT ) -SOmiPht « ) .1 1,F-0*) 

60 T'^ 425 




ir(At’S(Pl-I ) *oS(cPN(nPHI ).5<iM»PMT » ) 1 .P-fl*) 

I 00 TO <*25 

GO TO «2fl 

C ^HI SATUP*TI0N LOGtC 

• f2 CONTmit 

IF(«*OOP ,P0, \ » CPulz ni_T'lp/P«P^2 

IPIMOOE .eg. 2 > PP* Ir PLlP/ >’AP?2 

IE< ABSIOPHl) .6T, 2f>,/57,3» OP»<lrgO . /ST. (Q PMt ) 

ITT Aes«6PKIT .GT. iO,/S*,3jt>PHl = io77sT,SAS6N«noHl) 

GO TO G2S 

ITIWrifJkTuJfiT-iTN Co^T 

•2«i; CONTL'NUE _ _ _ 

irr*«OOC .EO. 1> D*LPH*= ntTOT'/PAoii 

irlMOP? .^0. 2 » PALPHAs P>j.TB/pARgl 

TP ills <D ALPHA ) .gT. 3, /St, 3) oTfuPHA = 3./S7,3*SfiH<nALPMi) 

G28 CONTINUE 

PRI~T 


PHI ♦ oPhT' 

Aj.Bj| A$ = alb I AS ♦ OAlPHA 

IF] VN 
VN 

AOATKx 
~JW 
IF 


SO.'/ST.S -2n./S7. ^•n^'(3.1AlS92A/l .2*«0.A.vwi 

Gli 0,8) TAnATK s 50 ./57.3 

LT, 5,2) TAnATK = 30,/57,3 

TA OATK ♦ auptaj 

,GT, Ae,/57.3) AoATK = A0./57.A 
^LT. ALPmXn ) aOATk s A LPMI N 


) 


noJftir 

aoatk 

IfTaBSjPHD .OT. 8S./57.')PHlr 8S./57.3*s;p,N(PHU 
SFt AySIA5 ,GT, 60./K7.3 . TAOATK ) ALBIAc;: GO./ST,3-TAOA Tk 
^ ALBIAS .LT, ALPMtN - TaOAT(« ) AL"1AS 5 ALP»*t*’ -TAOATk 




vO 


I 


IF 

_ *3fl_C^ TIHU £ 
F«7) 

c 

c 

c**« 


= AOATK 

7* • V **<*7.*7»t7*** *«****R77«a *.* _ *7At« *ai*»**«* *t *T. 

cpocsranc.e logic f^p *‘'On~2CRo ''C«;tPCP rROsPRANGES 


za 37G , S»V*W'*2 -l05.*VN 426. t 

IF IVN ,GT. FOB) GO To -33 

jtH =^8N<Xpt 
UTS s SGN»Y<5) - 
IFC ABS(UYB>^MISGN 


D^POCP) 

) .LE. l.c-io 


) 00 TO i»«0 

GO TO 
UAC 


•So 


<*<* 


«Su 


IF lABS(T<S) - P SRDCW)«W0/Gi76.l G 7 .OT. 2) 

IF (MS(T(S) -“ETshOCP) .67, SfTmAO) 00 TO 
PWI SGW s -PMJSgN 

^6^^n-;o3 

SWTNAO = ABS(Y(5) - P«RDCR)_ 

IF fSwTvAoiRO/GfitS.OG? .LE, 0,5') SWTMAO=0 , 5460t6.0R7/RO 
00 TO ISO 
^HISG' = -^(Vi 
VO8=VN»,03 
CONTINUE 

IF|APS«RLH0) .67, 8S./57.3) PmISANs .r^ON(RLHn) 

PHt = AB^iPMf)-* FiTISGN 
PHIHLOrPMI 


.LE, 


IF( AB<(PHl) 

F«8) r FhI 
f«761 s 2 

F<77) r vj 3)*'>n/6'>7ft, 007 
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9.12.4 A PRELIMINARY DESIGN FOR A BLENDED ENTRY DFCS 
INTRODUCTION 


This chapter presents a preliminary design for SSV attitude 
control during entry, using both the ACPS and ACS (aerodynamic 
control surface) torques in an integrated or ’’blended" approach. 

In addition to providing a means of minimizing ACPS fuel expen- 
ditures (by blending with the ACS torques), a prime objective of 
the design approach is to mnimize variations in the closed-loop 
attitude dynamics due to variations in the open-loop airframe 
dynamics. 

The material presented here is a summary of the design synthe- 
sis and gain computation documented in references 1 through 4. 
Current ch*anges to the latest documented design (given in reference 
1) liave been implemented in this work. 

This chapter is oiganized into three sectioxis: Section 1 de- 
scribes the trajectory parameters defining the entry operational 
environment of the NR 161C vehicle used for controller design; 
section 2 describes the longitudinal control channel; while section 
3 describes the lateral control channel, 

1, Trajectory Parameters 

The trajectory parameters presented here were generated by the 
Statistical Analysis Section (EG2), using a point mass simulation 
with pitch and bank angle modulation to control heati'.g, loads, and 
targeting. Given in Table I are the initial entry conditions for 
two NR 161C vehicle trajectories: A low cross-raige, lot^ down-range 

entry (CR: 370 mi; DR 4000 mi), and a high cross-range, high down- 
range entry CR 1150 ml; DR; 64p0 ni)# herein designated as trajec- 
tories A and B respectively. (Initial conditions are identical.) 
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9. 12.4 Entry DFCS (cont'd) 


Table I. - ENTRY TRAJECTORY INITIAL 
CONDITIOx\*S 


Paranic-ter 





d® 

K 

f 


Lat 


Long 


Value 

2 U, 1<09 I't/sec 
53 deg 
-0.9 deg 
395,052 ft 
90 deg 
0 deg 
81.53 deg 


Shown in figures 1 through 5 are the pertinent entry trajec- 
tory parameters, for both trajectories. Plotted agai ist Mach 
nimiber M in the f^st three figures are: Trim angle of attack 

dynamic pressure ^ , and flight path angle JT© . Plotted against 
time from entry interface t in the last two figures are: toch 

number M and total relative velocity'V^ . It should be noted that 
in addition to a large parametric variation along a parTicular 
trajectory, there is also a significant variation between trajec- 
tories. 

2. Longitudinal Control 

This section presents the control block diagrams and as- 
sociated gains for longitudinal control, or equivalently, control 
of the vehicle’s angle of attack , The basic design procedure, 
described in reference 2, involves a separation of the ACS and ACPS 
control loop synthesis, ’rith the addition of an interface logic to 
define the operational region in which the two control torque 
systems operate simultaneously. Further, the design of the ACS 
control loop is predicated on closed-loop dynamic insensitivity 
to variations in the airframe characteristics. 
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9.12.4 Entry DFCS (cont'c’) 


Section 2.1 and 2.2 present the block diagrams defining the 
aerodynamic surface (i.e., elevator) control logic and the ACPS 
(i.e., pitch jets) control logic, while Sections 2..3 and 2.4 pre- 
sent the fixed and scheduled gain parameters appropriate to these 
controllers. Finally Section 2,5 discuss the inputs, outputs and 
sample rates of the longitudinal control channel. 

2.1 Elevator Control 

Shown in figure 6 is the block diagram for elevator control, 
adaptedfrom reference 1. Basically, contro? consists of attack 
angle error and rate feedback to appropria’c ly position the elevator 
actuator. Use is made of both fixed and trajectory dependent gains 
for transient pitch control, in parallel with a clamped integr-ator 
to assure a trim capability. 

2.2 Pitch ACPS Control 

Shown in figure 7 is the block diagram for pitch ACPS control, 
taken from reference 1. As with the previous logic, attitude con- 
trol is achieved through the use of pitch rate damping, euid the 
appropriate torque signal is sent to the jet select logic. 

The interface logic between the ACS and ACPS subsystems is 
also shown in figure 7. The approach taken here is to inhibit firing 
of the pitch ACPS when there is "sufficient” pitch control acceler- 
ation available from the elevator. This determination of suffi- 
ciency is obtained from the commanded elevator deflection ) and 
hence, is a closed-loop index of elevator effectiveness. 

2.3 Fixed Longitudinal Gains 

Given below are values for the fixed longitudinal gain (and 
limit) parameters, which, from present simulation efforts, result 
in reasonable closed-loop vehicle response throughout the entry 
flight envelope. 
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Figure 6 - Elevator Control Block Diagram 
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figure 7 Pitch ACP8 Control Block 
Diagram 
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9.12.4 Entry Dt'CS (cont'd) 


Table II; Fixed Longitudinal Gain 


Parameter 

Values 

Dimension 

6U 

0.1 

JfC*' 

A 

0.4 

— 

ffmf 

10 


h 







5 





U 

or 



0.1S 


ix 

0.4 

mmm 


It should he n*^-ced t*iat the above gains are applicable to the NR 
161C vehicle, 

2.4 Schedule Longitudinal Gains 

As shown in figure 6, there are three variable gain parameters 
in the longitudinal control cliannel: > fa , and^. This section 

presents gain, schedules for the first two parameters and a simple 
limiting logic for the third. 

As discussed in reference 4> the gain parameters and f x 
may be plotted against Mach number and angle of attack, respect- 
ively, By coD 5 )aring the values for two different entry traject- 
ories (described in section 1) , a -piecewise lineeu? approximation 
may be made for «ach parameter, so that trajectory independent 
gain scheduling! is possible. The results are shCnm in figures 8 
and 9, with^, a function of Mach number andf^^a function of angle 
of attack. 
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9.12.4 Entry DFCS (eont'd) 


In order to avoid computational overflow problems due to div- 
ision by the dynamic pressui*e ^ , a simple limiting logic is used 
to set a lower bound on this parameter, and is own below: 




0.1 


y 


r 

/t- sfcpe 






Limiting Logic for Dynandc Pressure 

A 

In the above diagram ^^^is the dynamic pressure obtained from the 
air data, computer, whi.cL is lower limited to become the dynamic 
pressure ^ used in the controller logic of figure 6, 


2,5 Longitudinal Channel. Inputs, Outputs, and Sample Rates 

The inputs for the longitudinal controller fall into two cat- 
egories; those pertaining to the inner "control" loop, and those 
pertaining to the outer "gain" loop. These two sets of inputs, 
and their corresponding origins, are tabulated below: 

Table’ll Longitudinal Inputs 



1 COfOntOL LOOP ] 

1 LOOP 






%oc 




1 

rmtejyn 

ADC 

ADC 

ADC 
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9.12.4 Entry DFCS (cont»d) 


The outputs for the longitudinal controller are commanded elevator 
deflection,^*', and comnjand :?d pitch ACPS acceleration, it* • To 
obt'^in ele/on commandf ,the elevator and aileron commands (the latter 
from the lateral xogic) must be differenced in an eleven logic as 
shovm in figure 10. In a similar manner, commanded ACPS acceler- 
ations must make use of an appropriate jet select logic, which will 
be the subject of a future memoreuadum. 



There are two sample ratss presently in use for DFCS simulation 
efforts: a high-frequency control loop rate of 8hz, and a low- 
frequency gain loop rate of 0,41iz. It should be noted that this 
implies gain iqjdating every 20 control cycles, and may thus prove 
to be conservative. 

3.0 Lateral Control 

This section presents the control block diagrams and as- 
sociated gains for latei^ control, or, equivalently, control of 
the vehicles sideslip and bank angles, as with the longitudinal 
controller, the basic design procedure, described in reference 2, 
Involves a separation of ACS and ACPS control, with the addition of 
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9.12.4 Entry DFCS (cont'i) 


an interface logic. Further, the ACS oontroller utilizes an 
analogous procedure of gain scheduling to minimize sensitivity to 
open-loop ''shicle parameter cbaxiges. 

Sections 3.1 and 3.2 present the block diagrams defining the 
aerodynamic surface (i.e., rudder and aileron) control logic and 
the ACPS (i.e., y'aw and roll jets.) control logic, while sections 
3.3 and 3.4 present the fixed and scheduled gain paraiiwters ap- 
propriate to these Controllers. Finally section 3.5 discusses the 
inputs^ outputs and sample rates of the lateral control channels. 

3.1 Rudder and Aileron Control 

Shown in figure 11 is the block diagram for rudder and aileron 
control, adapte ’ from reference 1. Basically, control consists 
of sideslip and bank angle attitude and rate feedback to approp- 
riately position the rudder and aileron actuators, with the gains 
chosen so as to minimize lateral coupling and simultaneously 
provide adequate dynamic response with sufficient damping. The 
synthe sis procedure is described in reference 2; as can b® seen, 
the design incorporates both fixed and scheduled gains for lateral 
control with the ACS, 

3.2 Roll and Yaw ACPS Control 

Shown in figure 12 is the block diagram for roll and yaw ACPS 
control, taken from reference 1, As with the previous logic, at- 
tit’jde control is achieved through the use of sideslip and bank angle 


9.12-97 


j 



5j - 


^■</A’Sf«l77€S <»»yeA) 8K : (6) (W/TV StO^ j 

(7) (/v#nr . 

<&^ tAI#ry ; 

(!>; IMfirY SWfS ; 


J,p. ^ 

/,. • er 

/f2 * 

^3 = 


Figure - Rudder and Ailei on Control Block Diagram 


^.12.4 Entry DFCS (cont'd) 















^r"**g/KPf7- 
- ■^^ T^a^<Yjr 


.’ //d; ^2. ; 

(/O i)G4MA{> : 
Ot/TPtjr : 

^2} DeApeAOp ; 

output 

(fSj Su>pe A ’^ ; 


•^/S* ” 

O, tUa 

-^4 ^ -^^CvLo(t 
o,L\A, 

Ai = it^co^o^r 


I iEietr 
I 


Fipure 12- Latarel ACPS Control Block Diagram 













9.12.4 Entry DFCS (cont'd) 


i.. 

rate dajnping, and the appropriate torque signals are sent to the 
jet select logic. 

■ As discus’^e^ in reference 2, three control modes are provided 
for in the later.- ^C?S logic. Early in the entry, both aileron 
and rudder a:.- ’.nf^ff ective , requiring that both the roll and yaw 
jets be Used. Later in the trajectory, the ailerons become ef- 
fective, requiring that only the yaw jets be used. Finally, when 
the rudder eventually becomes effective, no ACPS torques are neces- 
sary, This modal logic is summarized ly the inhibit multiplications 
of figure 12 and the modal parameter definitions of figure 13. Note 
that, as with the longitudinal controller, surface effectiveness 
is measured by commanded deflection. 



.^0 ® ^ 

0 < /, / 

^20 “ /s' 4^5 


Figure /3 - Lateral Mode Parameters 



3.3 Fixed Lateral Gains 

Given in Table 4 are values for the fixed lateral gain (and 
limit) parameters, which, from present simulation efforts, result 
in reasonable vehicle response throughout the entry flight envelope. 


i 
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9.12,4 Entry DFCS (cont'd) 


Tabic 4 : Fixed l4ateral Gains 



As before, it should be noted that the above gains are applicable 
to the NR 161C vehicle. 

3.4 Scheduled Lateral Gains 

As shown in figure 11, there are sevon variable gain paramebers 
in the lateral control channel: g^ tlirough g]^Q (excluding gg),<T 
and Q . This section presents gain schedules for the first six; the 


9.12-100 



9.12.4 Entry DFCS (cont'd) 


dynamic pressure is limited in the same manner as dlcussed in 
section 2,4. 

As with the longitudinal gains, the lateral gains, obtained 
from reference 4> may be plotted against trim attack angle, and then 
fit in a piecewise linear sense. The trajectory independent gain 
schedules are shown in figures 14 through 19 with all gains funct- 
ions of attack angle. 

3.5 Lateral Channel Inputs. Outputs, and Sample Bates 

As ’vith the longitudinal controller, the inputs for the lateral 
controller fall into the two categories of control jjiputs and gain 
inpats. These inputs and their corresponding origins are tabulated 
below: 

Table 5: Lateral Inputs 



1 COmoL LOOP 1 Mw 

u 

iop 

ii^fpur 



P»r 1 $ 

o(r 



IMO 

1 AX> 

ADC 


The outputs for the lateral controller are commanded rudder eind 
aileron deflections,^ and commanded roll and yaw ACPS 

accelerations andit^ < To obtain elevon commands the aileron 
command mgst be differenced with the elevator command as shown in 
figure 10. As with the longitudinal controller, the roll and yaw 
ACPS acceleration commands must make use of a. ^et select logic. 

The sanqple rates for the lateral controller are the same as fcr the 
longitudinal chai_iel: a sampling frequency of 8 hz for the control 
loop and 0.4 hz for the gain loop. 
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9.13 


TRANSITION 


The transition phase starts at the end of the entry phase juSt 
prior to the pitbhover from a high ab|p.p of attack to the cruise angle of 
attack. This period of pitchover constitutes the transition phase. The 
transition phase ends with the attainment of the nominal cruise angle of 
attack. Ivo distinct functions are performed during ti'ansition: 

1. Transition from entry control which is basically 
RCS attitude control to cruising flight control, 
which is basically full aerodynamic control. 

2. Transition from entjy angle of attack (high a) 
to cruising angle of attack (low a) • 

The transition phase is character . zed by a transition in control modes and 
some form of flight path angle control to effect a smooth transition to a 
new equilibrium altitide. This siaggests that some changes to the guidance 
logic wl‘’l also occurr during this period. 

In this Section it; presented the first cut at a Ifcdfiod Digital Autopilot, 

A basic structure is provided for the purpose of encourging commonized coding 
for all control systems. A specific control law is included for only the 
Transition phase. The approach is currently being extended to other mission 
phases, and the Tremsitlon autopilot is boing upgraded as additional inform 
nation beccanes available. The autopilot structure presented here was in- 
tentionally made general enough to compass eJ.1 anticipated requirements. As 
specific autopilot designs mature for the various mission phases, the re- 
quirements for certain aspects of the coding presented here may not pateriali.’^e , 
For example, the bending state estimation is predicated on a reqoirei’jent for 
active control of several bending modes, a requirement that may not materialize. 
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phase. 

Module No. 00-4 Function No. 2 (MSC 03690 Rev. A) 

Submitted by; R. F. Stengel Co. MIT No. 8-71 
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%n3-l AUtijplIot 


1. INTRODUCTION 


The objective of the Unified Digital Autopilot Program 
is to provide rotational and translational control of the space shuttle 
orbiter in all phases of flight, from launch ascent through orbit to 
entry and touchdown. The program provides a versatile autopilot 
structure while maintaining simplified communications with other 
programs, with sensors, and with control effectors by the use of an 
executive routine /functional subroutine format. The program reads 
all e.xternal variables at a single point, copying them into its 
dedicated storage, and ccncrols its major support subroutines to 
be synchronous with the autopilot cycle. As a result, the autopilot 
program is largely independent of other programs in the guidance 
computer and is equally insensitive to the characteristics of the 
processor configuration (dedicated guidance computer vs. shared 
multi- processor). 

The unified autopilot program makes provision for 
sampling rates which are integer multiples of a basic sampling rate, 
using counters to establish the synchronous cycles. Extended 
computations carried out with a low repetition rate must be pro- 
vide I with pre-determined break points, in older that the low rate - 
and high rate - calculations can be interleaved. This requires pro- 
grammer control, but it has the advantage of precluding skipped or 
lost computation cycles (unless interrupting eirternal programs of 
higher priority monopolize the computer's time). 

The sequence of autopilot subroutines is arranged to 
minimize transportation lag, the time interval between receiving a 
measurement and effecting a control force. While this lag may be 
largely due to equipment external to the guidance computer, the 
time required for control computation can be significant. As a 
consequence, the state estimation computations are separated into 
two subroutines. State measurements are incorporated in the esli- 
ma1or )r filter in the early portion of the autopilot cycle, but the 
remaining state propagation or filter "push-down" does not occur 
Uiitil the autopilot commands have been written in the control 
effector output channels. 
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Autobilol/ 


(^ont'ld) 


In order to discuss the Unified Digital Autopilot Program 
beyond the level of the functional flow diagram, it is necessary to 
make specific reference to the autopilot for a particular mission 
phase. The transition from entry to cruising flight has been chosen 
for this purpose, as this phase makes use of all of the types of 
autopilot subroutines. 

The prototype transition autopilot described here can 
be considered an upper limit on the control system computational 
requirements, as it includes adaptive bending mode stabilization, 
high order state estimators and controllers, optimal state estimator 
gain computations, aerodynamic and reaction jet control, inertial 
and aerodynamic parameter estimation, and sensor /effector failure 
detection. It is almost certain that simplifications can be found; 
however, it is equally likely that unexpected complexities will occur. 
It should be emphasized that this is not a final transition autopilot 
design. 



9.13.1 Autopilot (cont‘d) 


NOMENCLATURE 



- 1 . 2 . 3 


H 


Mode stabilization estimation matrix (diagonal) 
Element of A 

Mode stabilization estimation matrix (diagonal only 
for uncoupled elastic modes) 

Element of B 

Mode stabilization estimation matrix (diagonal) 

Roll moment effectiveness ratio 
Pitch moment effectiveness ratio 
Yaw moment effectiveness ratio 
Element of C 

Mode stabilization estimation matrix (diagonal) 
Element of D 

Mode stabilization control matrix (diagonal only 
for uncoupled elastic modes) 

Rigid -body control vectors 

Altitude 


I Identity matrix 

iSl 2 3 4 estimator gain matrices 

k Ratio of specific heats 

k][ 2 Dynamic pressure estimator gains 

M Rigid-body measurement transformation m.atrix 

(longitudinal axis) 


M Mach number 

N Rigid-body measurement transformation matrix 

(lateral -directional axes) 

P Measurement noise covariance matrix 

(longitudinal axis) 
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9.13.1 Autpllot (cont’d) 


tP 

P 

£ 

q 

£ 

R 

r 

£ 

T 

T (H) 
U 


u 

;v 

V 

-yr 


V 


W 


w 


w 


y 



z 


stagnation pressure minus static pressure 

Roll rate (body axis) 

Measurement noise covariance matrix 
(lateral-directional axes) 

Dynamic pressure 

Covariance matrix of x 

Universal gas constant 

Yaw rate (body axis) 

Covariance matrix of v 

Mode stabilization sampling interval 

Air Temperature 

Velocity component along x-body axis 
(earth-relative, positive forward) 

Mode stabilization measurement vector 

Rigid -body disturbance input covariance matrix 
(lateral -directional axes) 

Velocity component along y-body axis (earth- 
relative, positive right) 

Magnitude of earth-relative velocity vector 

Perturbation state estimate vector 
(lateral -directional axes) 

Rigid-body disturbance input covariance matrix 
( longitudinal axis) 

Velocity component along »-body axi.i (earth- 
relative, positive down) 

Perturbation state measurement vector 
(lateral-directional axes) 

Perturbation state estimate vector 

Mode stabilization estimate vector 

Variance estimate for component of y 

Perturbation state measurement vector 
(longitudinal axis) 
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9.13.1 Autopilot (cant’d) 
a 

e 

6 

c 

e 

P(H) 

0 

n 

03 

Subscr ipts 

a 

Control 

DB 

DLC 

e 

Estimate 

Guidarce 

I 

i 

Measurement 

RCS 

r 

SB 

6 


Angle of attack 

Sideslip angle 

Control defection vector 

Bending mode estimatoi’ damping ratio 

Pitch angle 

Air density 

Roll angle 

Yaw angle 

Transformed bending mode natural frequency 
Bending mode natural frequency 


Aileron 

Autopilot output quantity 

Deadband 

Direct lift control 

Elevator 

Quantity estimated by autopilot 

Input quartity from guidance program 

Bending mode index 

Current sample of subscripted quantity 

Measurement input to autopilot 

Reaction control system 

Rudder 

Speed brake 

Derivative with respect to control deflection 
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9.13.1 Autopilot (cont»d) 


Special Notation 


( ) 
O 
(^) 

( )' 
()T 

(J 
{ ) 


Derivative with respect to time 

Mode stabilization control output 

State estimate before measurement update 

De-tuned bending mode estimation quantity 

Transpose of matrix or vector 

Vector 

Matrix 
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9.13*1 Autopilot (cont'dt 


2. FUNCTIONAL FLOW DIAGRAM 

The sequence for a single autopilot cycle is illustrated 
in the functional flow diagram of Figure 1. The autopilot program is 
called on a periodic basis vith a sampling frequency determined by 
the highest bandwidth control mode, which, in turn, is a function of 
the flight phase. Active bending mode stabilization will require the 
highest sampling rate; where this is not necessary, the sampling 
rate wOl be determined by rigid body control requirements. Ihitiali> 
zation branches r.re asynchronous, occurring only when the flight 
control mode changes or when there is a computer restart. All 
other branches are synchronous with the autopilot sampling rate, 
althou^ their sampling intervals may be integer multiples of the 
basic sampling interval. 

Tt'e basic subroutines of the unified digital autopalot 
are the followir.g: 

a) Sequence and foput / Output Initialization Subroutine 
This subroutine establishes the address of a list of 
subroutine addresses according to the flight mode, 
determines the basic sampling interval and th'i 
integer multiples for medium and slow sampling 
rates, and initializes all indices 

b) Re ad Subroutine 

This subroutine copies all inputs to the autopilot 
into dedicated temporary storage. The read list 
and sampling rate are functions of the flight control 
mode. 

c) Filter and Parameter Initial-ization Subroutine 
This subroutine initializes state and parameter 
estimates either at predetermined values or at the 
appropriate values read by (b) as required by the 
restart or control mode change. 

d) Bending Stabilization Subroutine 

This subroutine performs the minimum computations 
necessary to stabilize the betiding mode(s). R 
incorporates new measurements in the bending co- 
ordinate estimates, computes control commands, 
and, if bending pa*'ameter adaptation is performed, 
computes the avei iges required by the parameter 
adjustment law. 1 arameter adjustment, which may 
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9. ■>3.1 Autopilot (cont'd) 


be carried out at a lower rate, is performed in (h). 

In the present scheme, this subroutine does no rigid- 
body control. 

e) Perturbation State Estimation Subroutine - Pcirt I 
Measurements and computed ccrrtrol outputs are 
incorporated in the first part of the estimation sub- 
routine. Perturbations from the flight profile gener- 
ated by the guidance program are estimated here for 
use in the control law which follows. Only those 
perturbation states which are necessary for control 
or for display to the crew are estimated. This may 
include coordinates of bending or sloshing modes 

V '^ich are not sufficiently decoupled from the rigid- 
body modes for stabilization in (d). 

f) Rigid- body Control Subroutine 

This subroutine uses the outputs of (e) to derive 
control effector commands, which are intended to 
null the error between the actual and desired states. 
For proportional control effectors, e. g. , engine 
gimbals or aerodynamic control surfaces, the control 
law may simply consist of scaling and coordinate 
transformation. For RCS thrusters, phase plane 
s’.vitching logic can be used. 

g) Perturbation State Estimation Subroutine - Part II 
The calculations performed here prepare the estimator 
for the incorporation of measurements and control 
outputs on the next appropriate autopilot cycle. In 

the case of an estimator expressed as a constant- 

coefficient digital filter, this consists of 'pushing- 

down" the filter variables, i. e. . storing the i~ 
st 

value in the (i - 1)^ location , etc. For a time- 
vau'ying filter expressed as a state-space estimator, 
the .«?tate must be propagated to the next sampling 
instant, and revised gains must be computed (it is 
possible that the latter be done at a rate which is 
slower that the propagation sampling rate). 
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9.13.1 Autopilot (cont'd) 

h) Bending Parameter F^timation Subroutine 
Should mcertainty or time variation in the bending 
parameters be excessive, this subroutine will revise 
the parameter estimates, either by parameter 
scheduling or parameter tracking, as required. 

This subroutine can oe executed at a slow rate. 

i) T nertial Paraur.eter Update Subroutine 

Allis subroutine will alter ma«s, moment-of-inertla, 
and RCS specific control moment estimates at a 
slow rate as required for revision of estimator and 
control law constants. 

j) Aerodynamic Parameter Update Subroutine 

This subroutine will alter dynamic pressure, Mach 
number, state transition matrices, and aerodynamic 
specific control moments at a slow rate as required 
for revision of estimator and control law constants. 

The ca'*culations madi. in this subroutine are depend* 
ent on the flight control mode. 

<) Failure Detection Subroutine 

This subroutine examines the variances in the state 
and parameter estimates of the previous subroutines 
for "reasonability", detecting failures in sensors or 
control effectors by compa>nng the computed results 
to thresholds for normal operation. The subroutine 
is an adjunct to the failure discretes issued outside 
the autopilot program; since it depends upon an 
averagiiig process, it can bo executed at a slow rate. 

Although not explicitly shown in Figure 1, the option to branch to 
"Closeout" is available at each branch point. This precludes wasted 
testing when, for example, a bending stabilization computation is the 
only calculation required on a given autopilot cycle. It is also implicit 
that lengthy subroutines contain break points, allowing internal branch* 
ing on intermediate autopilot cycles. For example, an aerodynamic 
parameter update subroutine with sampling interval 50 times longer 
than the bending stabilization sampling interval can be entered at 
the higher rate, with computations proceeding from breakpoint to 
breakpoint on sequential passes. 
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9.13.1 Autopilot (cont‘d) 



Figure la FUNCTIONAL FLOW DIAGRAM 
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Figure lb FUNCTIONAL FLOW DIAGRAM 
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Figure Ic FUNCTIONAL FLOW DIAGRAM 
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9.13.1 Autopilot (cont-d) 

3. PROGRAM INPUT-OUTPUT 

The unified autopilot program will require inputs from 
inertial sensors, air data sensors, automatic guidance programs, 
manual controllers, failure monitor devices or programs, and the 
data management system's executive or control program. Control 
could be improved by inputs from control effectors, e. g. . aero- 
dynamic surface or engine glmbal deflection angles, and bending or 
slosh displacement sensors. Autopilot outputs will Include control 
effector commands, crew display variables, and parameters of 
value to guldamce and service programs. 

The input-output lists are dependent on the Qight 
control mode; the lists which follow are based upon the transition 
autopilot prototype. For this case. It is further assumed that 
smgular and translational state vectors are transformed to earth- 
relative. body-axis components before being read by the autopilot, 
that control effector sensors and elastic displacement sensors are 
not used, and that ths space shuttle vehicle Is equipped with con- 
ventional aircraft control surfaces as well as a reaction control 
system. The subject of what types of sensors are used to measure 
the vehicles state is not addressed. It is assumed that the transition 
maneuver is unpowered (no thrust input) and automatic (no manual 
input). Direct lift and drag controls are accounted for in the estimator 
but are not specifically commanded by the autopilot prototype. Mode 
and failure discretes are not specifically identified in these lists, nor 
are crew display variables. 
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9.13.1 Autopilot (cont'd) 


Input Variables 
H 

guidance 

measurement 

AP 


^guidance 

^measurement 

I* 

guidance 

measurement 

T(H) 

t 

t 

go 

U 


guidai-:e 


measurement 


guidance 


V 


measurement 


W 


W 


guidance 

measurement 


DI.C guidance 
^LC Measurement 
guidance 
guidance 
guidance 

®SB measurement 


8 


e 


guidance 

measurement 


e 


e 


guidance 


measurement 

p(H) 


Desired altitude 

Measured altitude 

Stagnation pressure minus static 
pressure 

Desired roll rate 

Measured roll rate 

Desired yaw rate 

Measured yaw rate 

Air Temperature 

Time 

Tirre-to-go 

Desired velocity along x-axis 

Measured velocity along x-axis 

Desired velocity along y-axis 

Measured velocity along y-axis 

Desired velocity along z-axis 

Measured velocity along z-axis 

Nominal aileron deflection 

Desired direct lift control setting 

Measured direct lift control setting 

Nominal elevator deflection 

Nominal rudder deflection 

Desired speed brake setting 

Measured speed brake setting 

Desired pitch attitude 

Measured pitch attitude 

Desired pitch rate 
Measured pitch rate 

Air density 
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Desired roil attitude 
Measured roll attitude 
Dosired yaw attitude 
Measured yaw attitude 

RCS thruster identification (one 
per commanded thruster) 

RCS thruster on-time - roll axis 

RCS thruster on-time - pitch axis 

RCS thruster on-time - yaw axis 

Aileron deflection 

Elevator deflection 

Rudder deflection 





9.13.1 Autopilot (cont’d) 

4. DESCRIPTION OF EQUATIONS 

The equations of this section are applicable to the 
control of an elastic vehicle in atmospheric flight, with particular 
reference to the transition phase of the shuttle's return from orbit . 
The transition phase actually consists of 2 transitions; an angle - 
of attack transition and a control effector transition. The transition 
from high-to low-angle of attack is necessitated by the need to 
meet heating constraints early in the flight and the requirement 
for inc. eased maneuverability toward the end of the flight. 

Similarly, reaction control thrusters provide attitude control 
during hypersonic flight, while aerodynamic control surfaces are 
used during the terminal portion of the flight. Each type of 
transition causes major changes in the dynamic characteristics 
of the system over a relatively short time interval. 

The transition autopilot prototype includes the following 

elements: 

a) Adaptive stabilization of 3 structural modes based 
upon a "classical" resonance filter. 

b) Linear optimal estimation of the rigid-body perturba- 
tion state. 

c) Gain- scheduled control commands to 3 aerod}mamic 
surfaces. 

d) Simplified phase-plane logic for backup control 
using the reaction control system. 

e) Inertial and aerodynamic parameter update, and 

f) Sensor and control effector failure detection. 

4. 1 Bending Mode Stabilization 

It has not been determined that active bend mg stabili- 
zation will be necessary for the space shuttle; however, should 
this be the case, a high oandwidth mode stabilization control loop 
will be required. Thi ; infers not only that sensors and control 
effectors have high bandwidth capability, but that the digital 
compensation be executed at a high sampling rate. If rigid -body 
control and mod? stabilization are combined in a single estimator- 
control law, the resulting number of computations may be excessive. 
As an alternative, it is suggested that the computations be partitioned, 
with the rigid and elastic control equations executed at different 
rates. A simple mode stabilization law is then executed at a high 
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9*13.1 Autopilot (cont'd) 


rate, while the more complex rigid-body control occurs at a slower 
rate which is consistent with the frequencies of rigid-body motion. 

The mode stabilization technique is divided into 3 parts; 
identification of the bending frequencies, estimation of the bending 
mode components in the state measurement, and control to oppose 
the bending deflection. The present design neglects the transporta- 
tion and dynamic lags of the control loop but could be modified as 
these quantities are identified. It also assumes tiiat inertial angle 
measurements are the only Inputs for bending mode control. 

The continuous-time resonant filter. 

2C(o.s 

L (1) 

s + 2 5 s + to!f 

where C is a damping ratio, and cOf is a radian frequency, has a 
frequency response of 1 when co = cOj, while its frequency response 
at zero and infinite frequency is 0. The sharpness of the resonant 
peak is determined by C. and the filter's tuning is fixed by fo^. A 
discrete-time realization of this filter is required for implementation 
in a digital system. 

The details of digital resonant and notch filters can 

be found in Ref. 1, where it is shown that the bilinear transformation,, 

s = (z-l)/(z+l), can be used to find the discrete-time version of 

eq. 1. This transformation eliminates the "folding" problems of 

the z -transform by mapping the entire s-plane into a horizontal 

segment with width ^ ir/ T rad/ sec (where T is the sampling 

interval), i. e., it establishes the frequency relationship, 

n = tan coT/2, and provides the corresponding digital filter as 
2 

well . The recursive digital filter corresponding to eq. 1 is 
yj « a[b(Uj - Uj_2^ ® yi-1 ‘^yi-2^ ' 

where 

a • 1 / (1 + 

b • 2COf 

c - 2 (1 - 
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9.13.1 Autopilot (cont'd) 


d = - (1 - 2 COf + 

Equation 2 has unity frequency response with zero phase 
lag when ic ~ and has zero frequency response for w - 0 or ir/T. 
Thus, equation 2 passes a sinusoid at frequency without 
modification, while frequencies to either side are attenuated. If this 
resonant filter is tuned to the bending mode frequency, then yj^ is a 
sampled estimate of the bending displacement, and a control command 
of the form. 


Ad. = e y., (3) 

is in-phase with the bending oscillation. By proper choice of the 
magnitude and sign of e, commanding the appropriate control effect- 
or with will oppose the oscillation. It will be recognized that 
this is a form of phase stabilization, and the lags mentioned earlier 
can not be neglected. This control law is proposed only for a case 
in which passive - or g'^in- stabilization is not sufficient. In the 
latter case, it would be sufficient to use equation 2 to form a 
subtractive notch filter^ 


u. = u. - y.. (4) 

to neglect explicit commands to the control effectors for bending 
stabilization, and to us ' u/ in the rigid-body estimation and 
control. 


The mechanism for the adaptation of the resonant 
filter is the tracking of a relative maximum in the power spectrum 
of the input signal. The variance of y^ forms an estimate of the 
power spectral density at 03 = OJr (with a bandwidth proportional 

* I 

to Z uJr). A resonant filter tuned to yields an 

estimate of t*^e power spectral density at oj = in the form of 

I I 

the variance of y^ . The greater variance ic assumed to lie closer 
to the spectral peak, indicating the direction for adaptation. Upon 
reaching the relative maximum, the frequency estimate will limit 
cycle about the proper value, a minor deficiency which could be 
rectified at the expense of adding a third resonant filter. A 
recursive estimate for the variance of y. can be found from 
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^ 1. 1 2 T , 
y, * y,., + k(y, -y,.,) 


( 5 ) 


where k is a constant between 0 and 1. and a similar equation 
can be applied to y^. 

It may be necessat'y to stabilize more than one 
elastic mode, and each mode may be present in more than one 
axis. In the worst case, it will be necessary to transform the 
measurements into normal cooruiuates to perform the bending 
mods estimation on the transformed coordinates, and to control 
each mode with more than one control effector. Then vector 
analogs of eq. 3 . 3. and 5 must be formed; 


£i 

= - "1-2* +£^Zi-i 

(6) 



(7) 



(8) 


In all cases. ^ ^ and £ are diagonal matrices 
with elements dependent on the damping and frequency of the 
normal modes. £.and ^ are diagonal only if u. and 
have identical dimensions and if the components of u are normal 
(in the bending modes). If this is not the case. 


i = iiS2 

is diagcmal and has clemems 2^^ . and ^ transforms u 

to normal coordinates. ^ scales ^ for proper control of tilie 
normal modes and transforms the result to control effector 
coordinates. Equation 8 is actually n scalar equations, as the 
components of ^ are indeoendent. 

For the flow charts of Section 5. it is assumed that 
there are 3 elastic modes. 3 attitude measuremente, and 3 
control effectors. 
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4 . 2 Rigid-Bod^ Pe rt^bation State_ Estim^ion 

For control during the transition, it will be necessary 
to have an estimate of the error between the desired state, as 
determined by the guidance program, and the actual state, as 
inferred jy the state measurements. The error, or perturbation 
state, can be estimated using a linear recursive filter whose 
general form is 


* iii-l 

(10) 


Ui) 


Equation 10 propagates the perturbation state 
according to the homogeneous equations of a physical model 
usin^, the state transition matrix, The effects of model 
errors and forcing ternr, s are added in eq. 11. Equations 10 and 
11 provide a completely general formulation for linear estimation; 
it is the means of determining Kj, and which defines the 
filter as optimal, near- cplimai, or ^ hoc . The Lunar Module 
rate estimator is an example of an ^ hoc estimator which can 

3 

be expressed as a special case of eq. 1C and 11 . 

For the transition phase an optimal estimator has 
been proposed, and the Section 5 flew charts make additional 
assumptions. The transition estimator is partitioned according 
to aircraft convention, i. e. , into a longitudinal set and a lateral- 
directional set. The estimator for each set incorporates all the 
variables which i.xe r.ccessary to solve the linearized equations 
normally used for the study of stability and control. Keeping in 
mind that the autopilot is primarily responsible for attitude 
control, it may prove possible to eliminate from the autopilot 
estimator some or all of the variables normally associated with 
"guidance", thus reducing the dimensions of each of the two 
estimators; however, this should not be done without assurances 
that the integrity of guidance -control interaction is maintained. 

It is assumed that the measurement variables have been transformed 
to the dssired body-axis coordinates before being read by the 
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autopilot; hence, no additional transformations are necessary, and 
the M matrix of eq. 11 is an identity matrix with the dimension 
of the measurement, z (which nted not be of the same dimension as 
x). 


The measurement vector of the longitudinal estimator 

is defined as 


_U - 


"U “ 


w 


W 


H 

- 

H 

» 

9 


e 


J , 

Measurement 

_e 

Guidance 


( 12 ) 


thus it includes axial velocity, normal velocity, and height errors 
as well as the conventional pitch and pitch rate errors. The 
control vector. 


& 


r 6 1 


r 6 1 

e 



6 



DLC 


'^DLC 



*SB 

_ ^ ®RCS _ 

Control 

0 


Guidarre 


(13) 


includes direct lift control and speed brake setting as well as 
elevator setting and pitch thruster rate, with the assumption that 
the guidance provides estimates of what these quantities should 
be (with the exception of reaction control rates). No assum;r*ion 
is made with regard to the source of the "control" quantities; 
they may be derived from autopilot and guidance commands or 
from control effector measurements. The latter is clearly 
desirable if accurate measurements are available; however, 
there is precedent for i.sing the former in the Apollo digital 
autopilots. 
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The estimated perturbation state, x, is 


X 


FAU 1 


AW 

AH 

A9 

A4 


Estimate 


( 14 ) 


The gains ^ and gg will be described in a later paragraph. The 
angle of attack estimate is derived from the total linear velocity 
estimates (summing perturbation estimates with nominal values) 
using tho conventional formula 


-1 

/-w\ 


tan 

Estimate 


tan ^ 

/-W 

/ Guidance 

W 

A Estimate 

\ ^Guidance ^ 

Estimate 


( 15 ) 


Similarly, the lateral-directional estimator includes 
the measurement vector. 


w 


" V " 


“ V ■ 



♦ 

r 

- 

r 

<P 


<0 

_ P _ 

Measurement 

_P 


Guidance 


the control vector. 


rfi 1 


fC 1 

a 

6 

6 

* 

r 


r 

^^RCS 


^‘^RCS 

_ ^"rcs _ 

Control 

_ ^RCS _ 


( 16 ) 


( 17 ) 
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Autopilot (cont'd) 


and the estimate vector 


formula: 


" AV ■ 


V = 


Ar 

A<0 


(18) 


L Ap 

The sideslip angle 


J Estimate 

estimate again uses a standard 


B - sin 


Total/ Esti 


= sin 


-1 


Estimate 
(V 


(19) 


+ V ) 
Guid A Esr 




The flow charts of Section 5 indicate that 
calculations of eq. 10 and 11 occur at different points in the 
autopilot sequence. In particular, eg. 11 occurs soon after 
the input is read, while propagation of the state (eq. 10) occurs 
after the control command has been issued. Tliis is done to 
minimize transportation lag between the computer input and 
output. Since eq. 10 is solved for the next autopilot cycle, the 
indices are incremented by one, and the equation appears in 
the flow chart as 


A'.though the longitudir 1 and lateral-directional 
estimators are each 5-dimensiom.., it is likely that many of 
the off-diagonal elements of ^ will be negligible, allowing a 
reduction in the number of additions and multiplications 
necessary to solve eq. 10 and 11. 
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Estimator gains for the transition autopilot 
prototype are computed in real-time and are the linear 

4 

optimal gains associated with Kalman and presented by 
"5 

Liebelt . The gain calculation is based upon the estimation 
of state covariances, £, which in this case also depends upon 
a ££iori knowledge of the disturbance input and measurement 
noise covariances. The covariance propagation, covariance 
estimate, and gain computation equations are 


Si ■ 

<1- Ma> 4. 

(20) 


= + ll+l. 

(21) 

*=i:=2^+l 

liii mI 'll ^ Ma III 

(22) 


Here the indices have been incremented by one, as the solutions 
apply to the next autopilot cycle. 

It will be noted that the gain matrix has been 
partitioned. Kj is the gain matrix which operates on the 
measurement residuals; ^2 operates on the control inputs, 
which are represented as estimated residuals in the gain 

g 

computations . The measurement tranformation matrix, M, 
is augmented as necessary to incorporate the control inputs 
(hence the subscript. A.). The disturbance input covaricnce 
matrix, appears in eq. 21; in general, this term prevents 
the gain matrices from vanishing under the influence of the 
inverse of the measurement noise covariance, P, in eq. 22 . 

The above notation is used for the longitudina. 
estimator; in the lateral- directional estimator, ^ S, and 
Yare substituted for £, and The state transition 
matrices are different in each case, being subscripted by x 
and V accordingly. 

Although it may not be apparent from the equations, 
this formulation of the optimal gain is well-suited to incorporating 
data at differing rates. As an example, let us examine the 
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longitudinal estimator, assuming that ^ is not measured 
and that only reaction control torques are present. Further 
assume that 9 is measured on each autopilot cycle, that U and 
W are measured every tenth cycle, and that H is measured 
every fiftieth cycle. The diagonal matrix is potentially fifth 
order; however, it is only fifth order when all measurements 
are taken simultaneously, i. e.,on an autopilot cycle when U, 


W, H, and 9 are measured simultaneously and when 
in non-zerc. When there is no control torque, M^ias a 
maximum dimension of four, but this occurs, at most, every 
fiftieth cycle. It is easy to see that measurements can be 
staggered so that the dimension of M is never greater than 
2 (assuming that the pitch attitutde measi rement, 9, is 
incorporated on every cycle). This simplifies the estimator 
greatly, particularly as a result of the reduced dimension of 
the matrix inversion in eq. 22. The price paid for the reduced 
dimensionality' on any single pass is the increased logic and 
storage required to account for the 5 measurement combinations 
(6 plus U or W or H or A9p^^g); however, tiie computation-time 
saving is considerable. This estimation scheme is demonstrated 
in Ref. 6. 


4. 3 Rigid-Body Aerodynamic Control 

It is assumed for the transition autopilot prototyrpe 
that the shuttle vehicle is equipped with three independent aero- 
dynamic controls — aileron, rudder, and elevator. In fact, each 
control surface generates secondary torques, and for the delta- 
wing configuration, the elevator and ailerons are likely to be 
combined in "elevens 

The postulated control law includes nominal aileron, 

elevator, and rudder commands obtained from the guidance 

program (where the commands may be stored, computed, or 

arbitrarily set to zero) plus perturbations derived from linear 

feed back control laws. The control gains may simply be stored 

functions of the time-tc-go, t ; however, the nonlinearity of 
7 

control effect, the wide range of angle of attack, dynamic 

8 9 

pressure, and Mach number. and the variation in *nomlnal" 
conditions from one mission to another may impose a requirement 
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for additional flexibility. This will be especially true if the 
guidance program is unable to provide nominal control settings. 

It is proposed, therefore, that the control gains 

be comprised of two parts; 1) a set of feedback gains stored as 

a function of t and 2) a factor to correct control surface effect- 

go 

iveness as a function of the estimated flight condition and the 
actual control setting. The gains dependent on t^^ can be obtained 
prior to the flight by minimizing quadratic costs on the state and 
control in the linearized model of perturbations along a nominal 
flight path. The control effectiveness factors are stored tables 
or functions which are dependent on the deviation of flight parameters 
from their nominal values. The control surface commands then 
take the form 


A A 

“ e. 

i 

= C (ff.M. 6^)fJ X. 
m e — 1 — i 

(23) 



(24) 

At 

“ r. 
i 


(25) 

A 

®i 

(A„ +A6^). 

® Guidance ® '■ 

(26) 

A 

^i 

(A^ +At). 

“Guidance “ ^ 

(27) 

®r. = 
1 

(A + A6 )• 

'^Guidance '* ^ 

(28) 


Hero, C^. C,, and C are not the conventional 
m -t n 

non-dimensional aerodynamic coefficients; they are the control 
effectiveness factors. The vectors , ^2 and f^ 
gains which are stored as functions of t 

go 

For the delta wing configuration, control moments 

about the three body axes will be supplied by the rudder^ and a 

pair of elevons. The specific control moments due to rudder 

deflection (6 ), left elevon defection (4 ), and right elevon 

,, ®L 

Some configurations have dual vertical tails and, therefore, two 
rudders. 
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deflection (i ), will be 
®R 


L = L(8 . 6 , 8 ) 

®R ^ 


M = M(6 . 6^ ) 

®L ®R 


N = N(8 , 4 , 6 ) 

®L ®R 


(29) 

(30) 

(31) 


Expanding to first order. 


L - L(6 .6 . 4 . 

r Guidance 


(32) 

®L ®R ^ 


M = M(6 .6 )_ . , 

®L ®R Guidance 


(33) 


+ A8 +Mj ^Control 

®L ^ ®R ” 


N = ‘r'caidance 


+ (N, A. +N A. +N Ae/ 
'l l 'r « ■■ 


(34) 
Control 


Then the perturbation specific control moments due 
to feedback control can be expressed as 


■ ‘'•e ‘'•c 

®L ®R ‘ 


M, M, 0 
^ ^ 0 

®L ®R 


^ Nj N, 

®L ®P • 


•> 





A6r 

m4 

H> «■ 


T • 

ill 


fT 

h ^ 


ill. 


(35a) 


J i 


9.13-29 



9.13.1 Autopilot (cont'd) 


or 






Aft. 


R 


and 


" "^6 ^ft. 


R 


M. M. 

o 0 

®L ‘ 


R 


®L ®R 


-1 


0 I2'' 


Ll 0 


T 

0 I3 


6 = 6 +A6 )• 

®L ®T e. I 

^Guidance ^ 


% • <‘e_ 

^Guidance ” 


6 = (A +A6 )• 

T* ' r 1 

I Guidance 


X 


V 


i 


'35b) 


(36) 

t?7) 

(38) 


The partials . are equal only for identical 


®L ®R 

mean deflections, as control effects are nonlinea", particularly at 
high angle of attack. This distinction is particularly important as 
a result of the coupling of longitudinal and lateral-directional 
control actions indicated in eq. 35b. The separation of control 
gains into time -dependent and flight condition-v'ependent compo- 
nents is retained in eq. 35b; however, it is clearly possible to 
evaluate the aerodynamic control partials for the nominal flight 
condition and con ^-ol settings, making them functions of time as 
well. 


Above a certain angle of attack, the rudder effect- 
iveness of most orbiter configurations will be negligible; thus, the 
3rd column of the partial matrix in eq. 35 is effectively zero. 
While the elevens (or ailerons) could be used to control yaw, there 
are only 2 control s',:rfaces for 3 axes, and the control of each 
axis is no longer independent of the other axes. This coupling 
could be beneficial ("proverse") or detrimental ("adverse"^ 'n 
general, it will be necessary to supplement yaw control at high 
angles of attack using the reaction control system. 
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4. 4 Rigid Body Reaction Jet Control 

The »*eaction control system used for attitude control 
during orbital flight and high-angle entry will serve as a backup 
control system during tne transition (as well as a primary system 
for yaw during the early transition phase). Consequently, it is 
sufficient to provide a rate damping and attitude limiting control 
law with wide deadbands fcr pitch arc roll (or sideslip), assuming 
that the aerodynamic control surfaces will exert precise control 
within the dead zones. A more precise parabolic-curve switch- 
ing logic is retained for high angle of attack yaw control. An 
important feature of the reaction control logic is that it not 
oppose the aerodynamic control inadvertently. Unless the lags 
in one system or the other arc excessive, opposition will be 
prevented by using the same attitude estimates for each set of 
control laws. 


The control law for the pitch axis commands the 
appropriate control jets when jAS] > 0Qg or when jA0| > 

At high angles of attack, sideslip angle and yaw attitude are very 
nearly normal-mode components of the lateral-directional 
mrtion?^ Thus it is appropriate to choose these as control axes, 
firing roll and yaw jets simultaneously when sideslip errors 
exceed their deadbands and firing yaw jets alone when yaw limits 
are exceeded. Sideslip rate can be approximated by a combination 
of roll and yaw rate errors, allowing the sideslip rate dampuig 
to command jets on when | A0 cos cr Atp sin a| > The 

sideslip limiting logic consists of commanding jets on when 
I jl| > ^^^sbolic switch curves similar to the Lunar Module's 

TJETLAW^ or to those presented in Ref. 12 would be used for 
yaw control. 


Below some angle of attack (to be determined) the 
sideslip control law degenerates into a redundant yaw contrc»Uer, 
and tne rudder becomes effective as well. At this point, the 
sideslip control law should be replaced by a body-axis roll law, 
firing on |A^| > or | Since the rudder becomes 

the primary yaw control effector, simple ^ and deadband logic 
then replaces the parabolic curves. 
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The selection and timing logic for the control 

3 

thrusters is similar to the Lunar Module logic and is not 

repeated here. The choice of selection logic is based on the 

number and positioning of control jets; should a highly redundant 

configuration be chosen, a linear programming approach could 

13 

be considered as an alternative to the LM design. Timing jet 
firings to the nearest millisecond, as in the LM-Alone autop.Aot 
(rather than to the resolution of the sampling interval, as in the 
LM's CSM-Oocked autopilot), is essential 1) when the reaction 
control system is the primary system and 2) when the mode of 
motion most closely associated with a control axis is unstable 
(this can occur for both longitudinal and lateral-directiciial modes). 


4. 5 Inertial and Aerodynamic Parameter Update 

Inertial and aerodynamic parameter updates are 
separated in the Section 5 flow charts because it will not always 
be necessary to update both at the same time. Although the 
inertial update is included in the transition autopilot prototype, 
this subrout ne need be done only when a main engine is firing and 
substantial changes in the vehicle mass are occuring (reaction 
control propellant usage can be ignored). In this case, the mass 
estimate is decreased as a function of the average thrust during 
the last autopilot campling interval, 

Macc. - Macc. , (Thrust) (39) 

L I* Jl L 

while moments of inertia are functions of the mass, 


XX. 

= (Massj) 

(40) 

yy^ 

= I (Mass.) 

(41) 

yy 1 

zz. 

= I (Mass. ) , 
zz l 

(42) 


i 


and reaction jet specific control moments ("jet accelerations") 
are functions of the moments of inertia. Both the moments of 
inertia and the specific control moments are reqt'ired, as these 
data are necessary for aerodynamic control as well as reaction 
control. The need to vary the relationships between mass and 


14, 15 
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inertia as functions of cargo bay loading must be determined. It 

is not anticipated that products of inertia will be estimated; 

however, cargo bay loading could have a significant effect on the 

"nose-up" product of inertia. J . 

xz 

The aerodynamic parameter s'loroutine Will update 
dynamic pres jure, q. Mach number, M. aerodynamic control 
gains, and the estimator's state transition matt ices. The 
dynamic pressure can be derived from the difference between 
measured stagnation and static pressure; it also can be derived 
from the estimated total velocity and the air density, p(H). 

A 'Wst estimate" can be formed by combining the two using 
parametric optimization. Assuming static pr<?ssure and air 
temperature are measured (or stored as a function of altitude). 


p(H) = p(H)RT(H). (43) 

where R is the universal gas constant, and 

q - kj [1/2 p(H)V|gJ +*^ 2 ^^ 


= kjqj+k2q2 


(44) 


The weightily constants, k^ and k 2 , are determined 
by the known variances in q^^ and q 2 : 


. 2 ,, 2 . 2 . 

k, - 

l, 2 2 . 2 . 

h - "aj%. > 


1 - k. 


(45) 

(46) 


Thf Mach number is determined from the total 
velocity estimate and the measured air temperature: 


M = ^gt/[kRT(H)]^/f 


(47) 


where k is the ratio of specific heats for air. 

An alternative formulation for q^ eliminates eq. 43 
and uses eq. 47. It is 
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Qj = I p(H)M^ (48) 

The control effectiveness correction parameters 
of eq. 23-25 or the inverse matrix of partial derivatives in 
eq. 35b are next computed as functions of a. M, and control 
deflection. The optimal gain vectors ^2* ^^^3 computed 
as functions of t . The product of control effectiveness correc- 
tions and the gain vectors is taken m the update subroutine in 
order to minimize computatiors in the more frequently computed 
aerodynamic control subroutine. Finally, the estimator's state 
transition matrices are updated as functions of a. Q> M, Mass 
(or Inertia), and t , using a dual approach similar to that for 
control gain computation. As in the case of the control gains, 
the purpose behind the dual approach is to maximize reliability 
and flexibility by incorporating both expected matrix values with 
measured values in the state transition matrix update. 

4. 6 Failure Detection 

This subroutine has not been examined in detail. 

The failure detection algorithms will monitor estimator residuals 
for abnormalities which would indicate failed sensors or control 
effectors. They will be tailored to the most-probable failure 
modes of the monitored instruments. 
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5. Detailed Flow Diagrams 

ThiB section contains detailed flow diagrams for the 
Unified Digital Autopilot Program '.nth enecific reference to the 
transition flight control mode. 
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Transition Read Subroutine 


(Read inputs according to Fast. 
Medium, and Slow sampling rate 
indices ) 

. Inputs: 

. Time 

. Nominal Control 

. Commanded State ( from Guidance ) 
. "Time-to-Go". t^^ 

. External Control Settings 
. Failure Discretes 
. Manual Commands 
. Measured State 

. Return to Executive 



FLOW DIAGRAM 
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I Effectors 
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2* "2 
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Figure 2d DETAILED FLOW DIAGRAM 
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Figure 2e DETAILED FLOW DIAGRAM 
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Output 6 
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Control Effectors 
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Transition Reaction Jet Control 


Pitch 


k|^’|>«pb or if I Ad command 
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Figure 2h DETAILED FLOW DIAGRAM 
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Figure 2i DETAIL'D FLOW DIAGRAM 





9*1.3. 1 Autopilot (cont*d) 
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Figure 2j DETAILED FLOW DIAGRAM 
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p(H) = Air density 
T( H ) = Air temperature 
AP = Stagnation pressure 
Static pressure 

j^Li- [j- La] " <*go* 


L«x- J ' *"■ “• V 


L =# 



Figure 2m DETAILED FLOW DIAGRAM 


Mass) 
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6. Supplementary Information 

The functional flow diagram of the Unified Digital 
Autopilot Program is largely independent of the mission phase and 
vehicle configuration; however, detailed flow diagrams depend 
on mission phase. The estimation and control algorithms are 
further dependent on vehicle configuration, while the parameters 
for control may vary with vehicle loading as well as flight condi- 
tions. Algorithms and control parameters will be further 
developed as mission profiles and vehicle characteristics are 
defined. 


In addition to standardizing functional flow, the 
Unified Digital Autopilot Program should allow sharing of sub- 
routines by autopilot logic develt^ed for various mission phases. 
For example, the transition estimator could be useful during 
cruising flight and landing approach. It also can provide 
guidance system monitoring during atmospheric boost and on-line 
estimation during orbiter powered ascent. Failure detection 
subroutines will be dedicated to specific measurement and control 
subsystems; thus the subi*outii.^ 3 will be employed as the 
subsystem? are required and may be used during more than one 
mission phase. In order for the autopilot program to be truly 
"unified", it will be necessary to consider autopilot requirements 
for all mission phases in a parallel fashiem. 

^nbire work will follow two specific paths. The 
unified autopilot concept will be refined through simulation and 
further analysiSi The autopilot requirements for specific mission 
phases will be pursued in accordance with the assigned tasks of 
NAS9- 10268, with particular emphasis on the requirements of 
high cross-range, deltr-wing configurations. 
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9.14 CRUISE AND FERIg CRUISE 

The 8rulse ^kase of the shuttle mission applies to the period 
of atmospheric flight between transition and approach. The Terry Cruise 
pnase refers to powered flight between airports. These phases include low 
alti tilde earth-relative navigation, guidance and aerosurface control. The 
Software functions required in this mission phase are the following: 

1. Powerdd flight navigation augmented by external deta to maintain 
earth-relative position. 

2. Powered flight giidance resulting in thrust level control (for 
Fc-rry only) and autopilot (attitdde) commands required to 
achieve desired intersection with the terminal approach path 
for landing. 

3. Autopilot '‘onputations to result in proper aerodynamic surface 
control and trim settings to maintain desired heading and flight 
path. 

9.14.1 Navigation TBD 

9.14.2 Guidance TBD 

9.14.3 Ciad.se Autopilot 


« 
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9.14.3.1 Mamifll Modftfl 


Tfi 6 e control system for the delta wing orbiter (as implemented by the Control 
Requirements Branch, GDC, m their Cruise and Landing Simulator) provides the 
following modes of operation: 

fl) Direct Manual (CM) 

(2) Rate Command (RC) 

O) Rate Command/ 'Attitude Hold (RCAH) 

(4) Automatic Guidance (AUTO) 

(5) Yaw Dan 5 )er (YD) 

Figures 1 through 3 present the complete control systems Implemented for pitc!!, 
roll and yaw, respectively. 

The systems presented Include certedn modifications over the previously imple- 
mented control systems (for the straight wing orbiter.) For the pitch and roll 
systems, the RC and RCAH control loops are no longer separate and e logic eleiLent 
has been Included to preclude switching to attitude hold with Icurge vehicle rates. 
In the pitch system, the attitlde hold function is not engaged unless the selected 
mode is RCAH, the hand controller is in detent (QIN), and the vehicle pitch rate 
(q) minus the ti^ coordinator piten rate (%<) Is within + .0175 rad/sec (QLIN) , 
Sij^ar conditions are also required to engage the roll attitdide hold function 
exoppt for the turn coordinate rate that is not requii?ed in the roll channel. 

The first parts of Figures 1 and 3 are presented to show the turn coordinator 
for pitch in the RCAH mode and for yaw in the RCAH mode with the yaw danQ)er dis- 
engaged. The turn coordinator as iii^)lemented provides improved arbiter response 
during banked manuevers. 

In Figure 3, a second feedback loop (in addition to the yaw damper) is shewn 
in the }caw control system. It has been included for studies to evaluate the ef- 
fectiveness of lateral acceleration feedback during cruise and landing. For base- 
line purposes, however, the loop gain (R 3 ) is zero and the loop may be disregarded. 

Table 1 presents the gains implemented for all three control systems. Table 2 
shows the logic expressions for controlling the numerous switching fimctlond used 
in the control systems. (The switch designations represent logical conditions, 
not individual switches^ and the same switch is employed in various locations). 
Integrators controlled by the logical switches are considered in reset, or tem- 
porary reset (tasO.Ol sec.), when the switch is closed. Attitude and trim commands 
are input to the integrator Initial condition tenuinals. The integrator gain is 
assumed to be plus one. 
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9.14»3.2 Digital Flight Control System 


INTRODUCTION 

The principal objective of the DFCS, horizontal aerodynamic 
flight phase, is to provide rotational control and speed brake con- 
tol of the Orbiter during those phases of the mission where the 
Orbiter flight resembles conventional aircraft. These mission 
phases are the post entry aerodynamic flight, cruise, and landing, 
and also the horizontal takeoff (without Booster). The aerodynamic 
DFCS contains an auto mode, direct manual, stability augmentation 
system (SAS) manual, and a rate command attitude hold RCAH mode. 

The horizontal aerodynamic flight phase of the DFCS is in a 
state of growth itself, even without regard to the incorporation of 
other mission phases. Modifications presently planned include 
additional tasks, such as anti-skid, brake control, nose-wheel t- 

steering, landing gear extend, throttle, a DFCS self-contained term- 
inal guidance system, and flexible body and fuel slosh control. 

This version of the DFCS provides the embrionic structure for 
the ultimate unified DFCS (see MIT Reports below) which wili. be used 
for both rotational and translational control of the Space Shuttle 
Orbiter and Orbiter/Booster ccanbined in all phases of flight, frcan 
launch ascent throu^ orbit to entry and touchdown. The DFCS 
will evolve into the unified DFCS by sequential incorporation of the 
presently separate DFCS's for the other mission phases. Thus, the 
present DFCS contains routines labeled such as: 

Transition Mode 

Entry Mode 

Orbit TVC Mode 

Orbit RCS Mode 

Insertion TVC 

Booster TVC 

and such routines, while presently empty, will contain those functions 
in the unified version. 

The DFCS provides a versatile autopilot structure while, main- 
taining simplified communications with other progrms, with sensors, 
and with control effectors the use of an exscutive-azid-subroutine 
format. (However, the present version of the DFCS ten5)orarily has de- 
moted these sub-routines, for convenience, to mere routines.) The 
DFCS reads all external variables (commands, sensors, and supplied 
data) at a single time point, copying them into dedicated storage, and 
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controls its major support subroutines (presently routines) to 
be synchronous with the autopilot cycle. As a result, the autopilot 
program is largely independent of other programs in the guidance ccm- 
puter and is equally insensitive to the characteristics of the pro- 
cessor configuration (dedicated guidance computer versus shared 
multi-processor ) . 

The sequence of autopilot functional computationB is arranged 
to minimize transportation lag, the time Interval between receiving 
a measurement and effecting a control force. While this lag may be 
largely due to equipment external to the guidance computer, the time 
required for control computation can be significant. As a consequence, 
the filter computations are separated into two sections, one being per- 
formed between the "read" and "write," that is. Just prior to the con- 
trol equations. The second filter section, including filter "push 
down," is performed after writing commands to the control surfaces, and 
therefore does not contribute to the transport lag. 

The DFCS is written in a basic Fortran version (IV) without using 
features peculiar to a particular computer installation in order to 
minimize conversion problems by other investigators. Flight versions 
or special test beds may require reprogramming into machine language 
to minimize computation time, at the expense, of course, of engineering 
readability. 


Documentation 

The baseline equations for the DFCS for the horizontal aero- 
dynamic flight phases of the Orbiter are presented herein as a Fortran 
listing, version 4 d. The listing contains a definition of the alpha- 
numerlcs, comment cards, and a compilation. cross reference printout. 

An auto*fIow chart 1s also presented, which was made from the exact 
listing shown. The listing and auto-flow chart ccxistitute version 4 d 
of the DFCS. Other explanatory diagrams herein, the six-page block 
diagrams of the control modes (version 4 b), and the DFCS interface 
diagram (version 4 b), are intended only to provide understanding for 
the more official listing (version 4 d). 

MIT Reports 

The DFCS was developed by the Draper Laboratory of MTT for MSC. 
The MIT project leader is Dr. Robert F. Stengel. Additional infor- 
mation describing the DFCS can be found in: 

a. Space Shuttle GN&C Equation Document No. 8-71 (Draper Labora- 
tory), "unified Digital Autopilot With Specific Reference to the 
Transition Rhase," by R. F. Stengel, March 1971. (Sec, 9.13.1 of this 
Report.) 
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b. SSV M.emorandum No. 71-23C-4 (Draper Laboratory), "Cruise 
Phase Autopilots for the Straight Wing Orbiter Vehicle," by A, Penchuk 
and R. W. Schlundt, May 1971. 

The first of these reports defines the intended approach for com- 
bining all the DFCS from the several mission phases into a unified 
DFCS, and is written in the context of modern control theory. The 
second report (memorandum) describes the servo design procedure and 
results whereby the analog/root-locus selection of servo loops feed- 
back signals, gains, and compensation were selected. The conversion 
from analog servo-loops to digital servo-loops using the difference 
equation approach is also detailed. The control equations for the 
delta wing orbiter in the present DFCS are based on this report, not- 
withstanding the nomenclature "straight wing" in the report +’itle. 

The CRulse And Landing Aerodynamic Simulator (CRALS) at MSC 
presently contains the DFCS (version 4 b) with minor modifications 
as required for the machine dependent interface signals. Copies of 
this variation can be obtained by request, both in list form, and as 
card decks. 


DFCS mterfacSs 

Figure 1 (which is largely self-explanatory) shows the interface 
signals between the DFCS and other units. The signals actively used 
in the present version are shown as solid lines. The dotted lines 
represent signals contained in the read and write statements of 
version 4d but which are not currently used due to zero gains loaded 
in their signal flow paths. 

Investigators who intend to luae the DFCS need not provide the 
"dotted line" group of interface signals, but if these interface sig- 
nals are not provided, they should also be omitted from the read 
statements to prevent program aborts. 

The guidance system is shown in Figure 1 as two separate systems; 
one for commands to the DFCS, and the other for computed information 
utilized by the DFCS. The computed information source need not be the 
guidance computer, but could be other sources such as the air data 
computer, or navigation system, or sensors. Of the latter set, speed, 
angle of attack, and altitude enable the DFCS to update control gains 
and filter constants. Speed and "range" (range to go until touchdown) 
enable the DFCS to control the speed brake during landing approach. 

One recurring interface problem traditionally is the arbitrary 
definitions for control surface position. Within the DFCS and at the 
interfaces, surface positions of elevator, aileron, and rudder are 
defined such as to cause positive pitch, roll, and yaw moments. This 
is opposite to* some common3y used sign conventions for elevator and 
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aileron. Eleven position polarity is positive for trailing edge 
upward, A right rudder pedal push will cause a positive rudder by 
DFCS definition, which is the same as pilot language for right rodder. 

The elevators gntl rudder follow the left hand rule, for thumb of left 
hand in direction of plus y or plus z direction. The purpose of this 
selection of signs was to simplify the verification of proper signs since 
positive gains result for the airplane trane**er function, A more coD5)alling 
reason for the sim selection was that in twitching from RCAH to Direct 
Mode, stick back (using standard signs) calls for plus pitch rate in RCAH 
but negative elevator position in Direct Mode. Using standard signs, 
numerous sign changes wor0.d confusingly appear and disappear with mode 
switching. The speed brake angle magnitude is defined in the DFCS as 
one-half of the included angle between the two surfaces. 

Interface signals are of two types, discretes and proportional 
variables. Figure i identifies the discretes and the values each can 
have. All others are proportional variables. The maximum and mini- 
mum permissible values of the interface signals have not been defined. 

The "restart" discrete is considered to be generated external to 
the DFCS. After re-initializing, the DFCS will reset the "restart" 
signal to zero, wnich informs all other programs that the DFCS infor- 
mation is trustworthy once more. 

The DFCS computes TNEXT, and informs the main computer when the 
next IFCS cycle is to be called. 

Control Modes 

The DFCS (horizontal aerodynamic flight phase) provides control 
at this stage of development for the elevons, rudder, speed brake, and 
landing gear extend. The major modes are (l) Kanual, and (2) Automatic. 

The Manual mode cont.ains three submodes: (l) Direct, (2) Stability 
augmentation system (SAS) - a form of rate command, and (3) Hate flom- 
mand Attitude Hold (rCAH). The SAS mode has been retained because of 
historical pilot familiarization with it and for purposes of evalu- 
ating comparative flight handling qualities. 

The manual modes receive commands frem the pilot’s hand controller 
and rudder pedal. In the Direct manual mode, the hand controller sig- 
nals are considered elevator and aileron position commands, lii SAS and 
RCAH manual modes, these same pilot's hand controller signals are con- 
sidered as attitude rate cemmands, pitch rate and roll rate. The rud- 
der pedals command rudder position in direct mode, and sideslip in RCAH. 

The Automatic mode receives Commands only from the guidance system. 
Provision is Included for many other guidance command signals which may 
be utilized in later versions if required by guidance. Version 4 d 
" reads" all the guidance commands, but the DFCS Internal gains are set 
at zero for those signals not currently used. 
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A yaw attitide or rate command to the DFCS la not used (at present) 
Isecause it would contradict the sideslip feedback for coordinated turns. 
Later versions may use yaw commands to avoid crabbing doring gusts at 
landing. 

A control option, which may be termed a mode, is "automatic" speed 
brake control and the elevator altitiiide trim control available only in 
the RCAH mode. Later versions will provide the auto speed brake option 
in all Tnan'UflT modes. It is s^ectable from the pilot control panel using 
FLAG 5. A trajectory profile of speed versus distance to touchdown and 
altitude versus distance, stored in core memory, is used together with 
externally provided information on speed, altitude, and range-to-go, to 
compute the speed brake commands and elevator trim commEinds. 

More detedls on this mode 6u:e contained below Intthe section entitled 
"Fortran Listing" near the end of that section (see the heading "Landing 
Approach Trajectory for Manual Modes"), 

Control Equations 

Sheets 1-6 contain block diagrams, in servo-loop format^ of the 
digitized equations for the Auto mode, a^ the BCAH mode, for each of 
three channels: elevator, edleron, and rudder. Each calculation block 
is a separate Fortran caM, except where several blocks have been en- 
closed in a larger, dotted-line block on the dleigram, which indicates 
one card. 

The calculation blocks are achieved in the DFCS programming by 
the relative time sequence of cards, which permits the DFCS cycle time 
lag to occur at the point so designated. No Fortran cards appear for 
the blocks. 

The filters are also shewn in Sheets 1-6, although located in the 
DFCS under headings called "Filter-1" and "Filter-2," Instead of in the 
"control law" routine. The filters are generalized second order digital 
filters, but at present, most of the filter gains are zero which results 
in first-order filters. 

The filter network following the summing junction of command and 
feedback is an Integrator, and a first order lead. This function is 
called (in alternative language) a "proportional plus Integral" function. 

Washout filters are so labeled on the diagram and correspond to the 
analog TS/TS+1, which is high pass filter with unity gain at high fre- 
quencies. The w6U9hout filter is used now only in the yaw rate feedback 
path for rudder control. It permits the yaw rate feedback to dangp the 
Dutch Roll node, but not to oppose a steady yaw rate turn. A pitch rate 
washout filter may be used in later versions to keep the vehicle nose 
up during banking maneuvers. 


% 
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RCAH/Elevator Channel (see Sheet No. The presently active 

part of this mode/channel has: 

a. Pilot *s hand controller longitudinal command of pitch rate. 

b. Pilot’s station trim beeper switch is a discrete inte- 

grated within the DFCS. This trims the pi '•oh rate, thus offsetting 
null biases in the hand controller, rate gyro, or any other so\irce. 

Bi direct mode, the trim button provides attitude trim. 

Pitch Rate Qyro .- Hiis signal is differenced with the hand 
controller conmand and the trim coraaand to provide an error signal. 

The pitch rate gyro signal is also "gained" and f e dback to an inner 
loop to provide serve dashing. 

The "attitude hold" aspect of the RCAH mode results only from the 
fact that with a zero r<tte cenmanded, the attitude command is "held." 

No s«ritching-in of an attitude gyro occurs at stick detent to ensure 
attitude hold. This has the advantage of permitting the RCAH mode to 
function even though all attitude references have been lost. Any 
drifting of attitude due to gusts can be corrected by pilot input. 

Additional integration in the servo loops are being considered to 
in 5 >rove the attitude hold features of RCAH. 

The noxmal acceleration feedback is not currencly used and there 
is a zero gain in the signal path. IMs feedback will be used in later 
versions to provide so-called C-star control during landing where the 
pilot's controller commands a combination of pitch rate and vertical 
acceleration. 

The normal accelerometer is not used In the Automatic mode. 

The elevator cosmand produced by the RCAH longitudinal channel 
is summed or differenced with the aileron command from the RCAH lateral 
directional channel to result in individual elevon coonands. 

The conversion of elevator and aileron signal to elevon signal 
is accomplished internal to the BFCS to permit usage of non-linear 
effectiveness gains of the elevens in future versions. 

In order to avoid stair-step motl<Hi of the control surfaces at 
the DFCS cycle period, several fixes have been proposed which utilize 
either local filtering at the actuator or rate conaands without 
position commands or some combination. 

Prefiltering doi^e locally at the rat.o gyro has been proposed to 
filter body bending without Increa,' ag t DPCS frequency. 
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Diiring bank angles (i.e., driring turna), the pitch rate gyro 
picks up body pitch rate and the RCAH gives a nose down output. To 
correct this, a vashout filter in the pitch rate gyro path can be 
employed, or a moie direct cross feed from the lateral RCAH could be 
added. The DFCS present listing contains the pitch rate washout 
filter in anticipation of this purpose, but the signal is temporarily 
deadended, as shown in Sheet no. 4. 

RCAH/Aileron Channel (sheet $) and RCAH/Rudder Channel (Sheet 6) ."» 

Together, the rudder and aileron channels constitute the lateral- 
directional RCAH control. The rudder channel is considered an inner 
loop for the aileron channel. The loop closure sequence is: 

a. Yaw rate gyro feedback to rudder, with gain and washout filter. 

b. Sideslip feedbadt to rudder with gain. The "filter" at 
present is inactive as such. 

c. Roll rate feedback with gain (aileron channel) for damping. 

d. Roll rate feedback (aileron channel) differenced with pilot *s 
hand controller to produce the RCAH roll rate error. 

e. Two feedbacks shown are presently inactive. Future versions 
may use yaw rate feedback for cross feed into aileron, and roll rate 
cross feed into rudder. 

The forward path, from roll rate command error to aileron, contains 
a gain, an integrator, and a first order lead. 

As discussed above, the pilot’s hand controller provides roll 
rate commands, and the rudder pedals provide sideslip command. 

Manual trim is provided separately fear roll rate and sideslip 
in the RCAH mode. 

Future versions may substitute lateral acceleration for the side- 
slip feedback, due to possible difficulties with sideslip sensors, 
with corresponding changes in compensation. 

Auto Mode (sheets 1, 2, and 3).- Lcbs detailed discussion of the 
Auto mode will be required since it is essentially a position loop 
closure around the RCAH mode. 

Comparisons of the Auto mode and RCAH ofode control equations can 
be made by viewing the corresponding diagrams together; i.e., Sheets 1 
and 4 for elevator. Sheets 2 and ^ for aileron, and Sheets 3 and 6 
for rudder. 
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Auto/Elevator vs. RCAH/Elevator (Sheets 1 and 4).- The Auto 
cc md Is pitch attitude from the guidance system. The forward path 
cc .uains a gain, integration, and first order lead, all of which operate 
on the pitch error signal. Damping by the pitch rate gyr' is provided. 
Inactive in the present version are the guidance ccnmands and feedback 
for angle -of -attack, and guidance command for pitch rate. 

Itonual trim is not provided in the auto mode since the function 
is accomplished internal to the guidance system. 

Auto/Aileron vs. RCAH/Aileron (Sheets 2 and ^) . - Similar remarks 
to those contained in the preceding paragraph apply to the aileron 
channel (with an appropriate change in axis ncmenclat\ire). The guidance 
system coimands for roll rate, yaw rate, and yaw attitude are currently 
inactive in the auto mode. The yaw rate gyro cross feed into the 
aileron channel is also currently inactive. 

Auto/Pudder vs. RCAH/Rudder (Sheets 3 and 6).- Presently Inactive 
for the Auto/rudder channel are the guidance commands for yaw, sideslip, 
yaw rate, and roll rate. The following cross-feeds are also inactive: 
Roll rate gyro and yaw attitude gyro feedback. 

DECS Listing* - sections of the progeam in sequence are: 


1 Entry point (for all cycles). 

103 "Pad Load"mitialization. 

If flagword ITURN = 1, this indicates an initial turn on 
Constants will be loaded from some external source by 
"read" statements. 

104 Begin standard cycle. 

External flagwords are read which indicate restart, modes, 
gain update option, and automatic speedbrake. 

Read clock and compute next DFCS cycle start time. 

2k Logic test to deteimixeif DFCS needs re-inltiallzatlon. 

If there has been indication of a restart, or a mode change, 
the DFCS must be re-inltiallzed. 
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Sequence and l/o Initialization. 

Reset the restart flags, change to the nev modes, and set 
an internal flag (ISTART) to indicate that this DFCS pass 
is an initialization cycle. At the end of the entire 
DFCS cycle, this flag will be set back no normal. 

Establish Sampling R»bes 

The DFCS has three different cycles. The control laws 
and filters are performed on each fast cycle, now set at 
.10 seconds. The speed brake commands are calculated on 
each medium cycle. The gains are updated on each slow 
cycle. A table look-up sets the three cycle times as 
functions of the DFCS mode and mission phase. 

The half-time TF2 for a fast cycle is used in the digital 
filter gains. 

The ratios of the three cycle periods are computed so that 
countdowns can be used to determine which fast pass coincides 
with a medium or slow pass. The Modulo function is used 
for the counting at the very end of the program. The 
initial count for the medlim and slow pas^s are offset 
by one cycle to avoid having a medium and slow pass occur 
simultaneously. 

For the initial pass, the gains are set equal to the 
basic gains in table KFIX, and filter constants are set to 
the basic GFIX constants. Subsequent passes will update 
the gains and constants as functions of speed, altitude, etc., 
which are stored in other tables. 

Initialize indices. 

These are the Indices for the Modulo functions which count passes 
to identify medium and slow cycles. 

Branch to read, based on mode. 

Each mode will result in reading a different set of com- 
mands, sensors, and other data. 

All "reads" are located below in a read subroutine at state- 
ment 300. After reading, the program returns to the executive 
section for the next branching test. 
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22 Branch to filter and parameter initialization. 

A logic test is made to determine if initialization should 
be bypassed. Otherwise, a branching is made on mode, 
initialization is made, and control is passed to the next 
branching. 

Whenever a branch test is made on mode, and the result is 
Manual mode, a second branch test is made on Manual mode, 
to determine which manual submode. 

Similar branching and ret\imlng to executive control is 
made for: 

56 Branch to state filter - part 1 (filter update). 

57 Branch to control law. 

59 Branch to state filters - part 2 (filter pushdown). 

60 Branch to parameter estimation. 

This branch depends on mode, and also if Flag has 
requested parameter estimation. 

The parameter estimation routines ai’e located below at 
statement 7051. These are filters for external informa- 
tion for angle-of-attadc, velocity, air density, and dynamic 
pressure. 

Afterward, table lookup is performed to obtain the iqJdate 
for the control lav gains. 

80 Branch for closeout, 

A different closeout is used for the intlal pass. 

This completes the branching section in the executive program. 

The "sub-routines” follow. In this version, as stated above, 
the subroutines are routines. 

Cooments will be made only to supplement reading of the listing. 
2000 Initialization Routines. 

87 PuUup and flare constants for manual approach. 

The strategy is to let the DFCS compute these control system 
constants to avoid human calculation errors when trajectory changes 
are made. 
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UOO Filter Routines - part 1. 

All possible filter calculations are done in Filters part 2 
following the control law calculations to avoid transport lag. 

Filters part 1 is performed immediately after "read** and 
prior to "control laws." 

The filter equations and identification of constants by 
symbol can be seen on Sheets 1 to 6 of the block diagrams. 

5000 Control Routines. 

The control laws for Auto mode and RCAH are shown in the 
block diagrams Sheets 1 to 6. 

511 Multi-Rate Speed Brake Control .- An automatic speed brake control 

mode is available (only in RCAH) to assist in landing 
approach. A reference trajectory is ctored in the DFCS 
based on altitude and velocity versus range to go to touch- 
down. Error signals of velocity and altitude are simply 
**gained** (without integration or compensation) and added to 
the pilot’s manual speed brake command. The altitude gain 
is set at zero. Additional information on this subject is 
below in the routine **Landing Approach Trajectory for 
Manual Modes** at statement 702. 

512 

^ , Elevator Trim for Landing Approach Trajectory .- DECK and DECU 
^ are the delta elevator commands due to altitude and velocity, 

respectively, from the "Landing Approach Trajectory," for 
RCAH only. The altitude and velocity errors are **gained" 
(without integration or compensation) and added to the other- 
wise total RCAH elevator command. The velocity gain is zero 
at present. The altitude error produces a low gain, slow 
loop elevator trim which tends to keep the vehicle on the 
glide slope. The loop is so slow that pilot commands to 
the RCAH are much faster. Hence, the pilot and not this 
automatic feature produce the pullup maneuver and landing 
flare. More information on this subject is below in the 
routine *'Landing Approach Trajectory for Manual Modes" at 
statement 305. 

6000 Filter Routines - Part 2.- This routine contains all the filter 
calculation that can be done before the "read" on the next 
cycle. Also included is the "push down" of the digital 
filter where the quantity at becanes the same quantity 
at for the next cycle. 

The filter equations for the RCAH and auto modes are shown 
in the block diagrams, Sheets 1-6. 
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611 Medium Cycle , filter for X- velocity. 

612 Slow Cycle , filter for altitudes. 

702 Landing Approach Trajectory for Manual Modes .- The reference 

trajectory is used only in RCAH mode, but must be maintained 
in all modes in case of switchover to RCAH. Two trajectory 
curves are stored, velocity versus range, and altitude ver- 
sus range. The velocity trajectory assists in speed brake 
control of velocity. The altitude trajectory assists in 
elevator control of altitude. 


The range-to-go is broken into four segments. Trajector'y 
information is stored at the endpoints of each segment for 
velocity, altitude, and altitude slope (slope with respect 
to range). IR is the index 1 to 4 for each segment of 
range, with 1 nearest touchdown and 4 at high altitudes. 


For altitude, the segments are; Landing flare « 1 (ending 
in touchdown with a sink rate), straight line glide slope = 2, 
approach pullup = 3, and steeper glide slope straight line = 4. 
The glide slope and pullup segments for altitude are algebraic 
quadratics (a + bx + cx^ + dx3) such that the slope and 
position of the curve coincide with the adjacent stral^t line 
segments. Coefficients for these Interpolations are computed 
in the DFCS in the Initialization routine "Generate Pullup 
and Flare Constants for Manual Approach," located near 
statement 2010. 


The velocity segments of the reference landing trajectory 
are straight lines for the three nearest touchdown. The 
fourth segment is exponential and results in a linear IAS 
(indicated airspeed) versus range. 

RGO is the range-to-go (to touchdown) and is navigation type 
data supplied to the DFCS by landing navigation aids or by the 
guidance system. RLAITO (l) is the range-to-go (reference 
trajectory) at the endpoints of each range segment. RGO 
segment nu^er IR= 1 is bounded by endpoints 1^1 and I » 2. 
Segment IR = 2 is bounded by endpoints 1=2 and 3. The 
program determines IR by comparing RGO - RLAND(I) until it 
goes negative - then IR = I - 1. 

Having found IR, the program proceeds with the interpolation 
formulas. First, however, each segment of range has its own 
measure of reference length RL, which, as the listing suggests, 
RL = RGO - RLAND(IR). 
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The results of the interpolation equations are UG and HG, a 
velocity ccamnand and altitude command. These ccmmands are 
used in the control routines (described above), to provide 
speed brake and elevator commands. 

The next version, or soon thereafter, will contain a Mach- 
trim featuas at aloitudes above 45,000 ft to provide a 
steady-state angle-of -attack along a reference t-'ajectory 
in the RCAH mode. 

8l Closeout Routine .- ISTART is the flag which indicates if this 

pass was an initialization pass for the DFCS. If this was 
a "first" pass, the executive program will begin the closeout 
routine at statement 8l and set ISTART = 0. 

For the first pass, both the count indices for medium and slow 
are 0. IS = IS + 1 offsets the slow cycle index from the 
medium cycle index so that the medium aid slow cycles will not 
occur on the same pass. 


82 Increment Indices .- Each fast pass increments the counting in- 
dices for the medium and slow cycles. 

The Modulo function is a comparison function used in count- 
downs such that if the index exceeds or equals the second 
argument, the first argument (index) is set to zero. For 
example. Modulo 5 would count 1, 2, 3, 4,0 and repeat. 

When IM = 0, the DFCS will perform the medium pass operations, 
in additic'' to the fast pass operations, and similar results 
occur for the slow pass when IS = 0. 

Autoflow Chart 

The autoflow chart is a flow diagram made automatically 
l?y computer program. Hand notations in ink have been made on the 
first page of the diagrams to explain the autoflow symbols. The 
computer mainstream of flow does not always enqphasize the same 
mainstream that the human designer had in mind, so the results 
are sometimes confusing. However, the autoflow chart is a useful 
tool. 


Pad Load and Constants 

The constants have been verified for the NR 161C. Now constants 
are being determined for the MSC 040A Orbiter. These constants have 
nee been included because they are changing, but are available upon 
request. 
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00015400 


THETA 
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00015900 
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00017200 
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00017800 




ISN 0003 


COMMON 


TNEXT ,ORM,CEM ,0AM t OSB M, RUTR IM , ETR I M, ATR IM , 

OR G , OF L G « C E R G , 0 S P G , ALP HG , T H FT AG , QG , B FT A G , P S I G , RG , ^ 

UG,U,HG,H, IM,I S,BETA,AO^AI ,THO,THI ,00,01 ,BO,BI 
FSC,PSJ,R0,RI,PH0,PHJ,P0,P1,N0,NI,U0,UI,H0,HI,GFAR 


00018300 

00018400 


00018900 

00019000 
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1SN 0004 
TSN 0005 
!SN 000ft 
ISN 0007 
TSN OOOfl 

ISN 0009 


ISN 0010 

ISN 0011 
ISN 0012 
TSN 0013 
TSN 0014 
ISN 0015 
ISN 0016 
ISN 0017 
TSN 0018 
fSN 0019 
TSN 00?0 


ISN 0021 
ISN 0022 
ISN 0023 


ISN 0024 
ISN 0025 
ISN 0026 
TSN 0027 


ISN 0028 
ISN 0C29 
ISN 0030 
ISN 0031 
ISN 003? 

ISN 0033 
ISN 0034 
ISN 0039 
TSN 0036 
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ISN 0039 
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9001 
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FORNAT U0F8.0) 

FORMAT fF8.0*l0I5> 

FORMAT f20K,8F12.4| 

RFAL K,KFIX,MEDIUM,NI ,NC,NSUF 

INTEGER FLAGl , AG? . FL AG 3, FL AG4 , FL AG5 , GE AR 
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1 CONTINUE 

C»>»FL1GHT CCNTROL SYSTEM EXECUTIVE R0UT1NE<<<<< 

C «»«*«« ■mm-mmmwm mmmwmmm 

GP TP (103,104), ITUPN 
C»»>"PAD LCAD*» INITIALIZATICN 
10? CONTINUE 

KFIX 
CFIX 

S8, FL, SLOPE, RL AND, HL AND, UL AND, RHO, SCALE 
FAST, MEDIUM, SLOW 

efit,pfit,phfit,qeit,rfit,thfit,pcomp,qcomp,eil 


READ 
READ 
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RF AC 
RF AC 
TF a 
FLAGl 
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104 CONTINUE 
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IF (FLAGl) 13,1?,?0 

13 IF (FLAG2-M0C£) 20*14*20 
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19 IF fFLA43«fl«IIIIOOI 20*30*20 

C»»>SEOUENCE C I/C INITIALIZATION 
?0 FLAGl a 0 

WRITE (6*9001) FLAGl 
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READ ISrSOGO) ORtOEL t OER tOSB* ALPH»BETA ,p* Q» R«U« H» RGO» FAIL 
GO TO 22 

C»»>TRANSITirN NOCf READ 
520 GO TO 22 

C»>>>ENTRY MrOF REAR 
330 GO TO 2 ? 

O»>>0RBIT 7VC MODE READ 
340 GO TO 22 

r.»»>ORBIT RCS MODE READ 
350 GO TO 22 

C.»>>>INSEPTirN TVC MCDE REAC 
360 GO TO ?? 

C»>»BOCSIER TVC MCDE READ 
370 GO TO 22 

C ♦♦♦♦♦♦♦♦♦♦»»»»»»»**»4 *****» **<>44*»****4***4»»4*44*4**»»4»*»»»»<i*» 

C»»>INITIALI7AT10N ROOTINES<<<<< 


00028700 

00028600 

00028900 

00029000 

00029100 

00029200 

.00029300 

00029400 

00029900 

00029600 

00029700 

00029800 

00029900 

0003(N>00 

00030100 

000302fy0 

00030300 


ISN 

0062 

2000 UG - U 

Q0030' 00 

ISH 

0083 

HG • H 

00030*00 

ISN 

0084 

NStm « 0. 

00030700 

TSN 

0085 

NOf2» » 0. 

0d03<7800 

ISN 

0086 

Nil?) = 0. 

00039900 



C»>>>GENERATF PULLUP E FLARE C06STANTS FOR PAMJAL APPROACH 

00031000 



C>»»ALTITUD£ CQNSIANTS 

.00034100 

ISN 

0087 

SLOfll • SLONE 

000/11200 

ISN 

0088 

DO 2010 t«l*2 

00031300 

ISN 

0089 

J - ?*1 

00fi91400 

ISN 

0990 

SLO(I*U = (HLAND(J*l)-HLANC( jn/(RLANP( J^1)-RLAND{ J)» 

00031500 

ISN 

0091 

FHI 1,1) = HLAND( J-1 ) 

00031600 

ISN 

0092 

FH(2,n = SLOm 

00031700 

ISN 

0093 

RL » RLANO(4»7«LANOf Joll 

ocimtoo 

ISN 

ISN 

00^‘>- 

RL2 

OiJtiiOOO' 

ISN 

0096 

HLP = SL0(I*n-FH(2,I) 

00032100 

ISN 

0097 

FHO.n = (’.♦HL-RL*HLP)/RL2 

00032200 

ISN 

0098 

2010 FH(4,n - (RL*HLP-2.*HL)/(RL2*RU 

00032300 



C»»>VE10CITV CONSTANTS 

00D32400 

ISN 

0099 

00 2012 I-lr3 

00032300 

ISN 

0100 

2012 OStOm • f01ANf0fl4ll-ULANe<IM/fRlAN0(T»!l-RtAN0f 111 

f300926O0 

ISN 

0101 

USL0(4) = AinG(ULAN0(5) /LLAN0(4) )/(RLANC(5>-RLAN0(4) ) 

00032700 



C>>>>>AUTO MODE 

0003280'" 

ISN 

0102 

210 ELBUS = 0. 

0003290O 

ISN 0103 

S8IAS •* 0. 

00033000 

ISN 

0104 

N'TSUN 9 0. 

00033100 

TSN 

0109 

Tf<S«N • 0. 

00055200 

ISN 

0106 

PSSUM = 0. 

00033500 

ISN 

0107 

ASUM = 0. 

00035400 


9 . 14 « 3.2 ^|CS (cont'd) 



»T»T«lcVl.-I«T« 








00035700 

00035800 


ISN 0130 
TSN 0131 


00036300 

00036600 


OSBC « 0. 
OHLAST * 0 


ISN 0136 
ISN 0137 


00036900 

00037000 


2136 CONTINUE 


ISN 0162 
ISN 0163 


00037500 

0003760G 


ISN 0168 
TSN 0169 


00038100 

00038200 


C>»»ORBIT TVC NQCE 
260 60 TO 60 


ISN 0152 


s 




















PAGE 009 




C>»»ORBIT RCS MOCE 

00038300 

ISN 

015i 

250 

GO TO 40 

00038400 



C»»>INSEBTION TVC NCDE 

00038500 

TSN 

0154 

260 

GO TP 40 

00038600 



C>>>>>BnnSTER TVC NCCF 

00038700 

rsN 

0155 

?70 

GO TO 40 

00038800 



C 

««]*«*««******« m**:^^^^^!************* *************** ************* 00033900 



C»»>FILTER R007INES-PART l«<« 

00039000 



C 


00039100 



C»»>MANUAL FflTFRS-PAPT 1 

00039200 



C 

AOO TRIM DISCRETES 'rC MANUAL COMMANPS 

00039300 

ISN 

0156 

400 

N0( 1) = Gm )*NI ( D+NSUP 

00039400 

ISN 

0157 


ATRIN - K(15)*ATR1M 

00039500 

ISN 

0158 


ETRIN - KI16I4ETRIN 

00039600 

ISN 

0159 


RUTRIN > K|t7r«RUTRIM 

00039700 

ISN 

0160 


PTRIW « TF2*f ATRIN^TR IMLfl II ♦ TRIMPMI 

00039800 

ISN 

0161 


QTRIM = TF?*(FTRIM*TRIML(2>> ♦ TRIMCI2) 

00039900 

ISN 

0162 


RTRIP = TF2*IRUTRIM+TRIPLm )♦ TPIMCI31 

00040000 

ISN 

0163 


OAC « KC3I«tOAM4PTRIH) 

00040100 

ISN 

0164 


0€C • K(2>«IOeN4QTRIN) 

00940200 

ISN 

0165 


ORC » Kfli«fONN«NTRINI 

0004030# 

ISN 

0166 


GO TO I50#4««f>4©02I,N»NN00 

08940400 

I SN 

0167 

4002 

Pirn = DAC-P 

00040500 

ISN 

016B 


Oim = DEC-0 

00040600 

ISN 

0169 


RUll * -R 

00040700 

ISN 

0170 


SIfll « >KTA 

00040800 

ISN 

0171 


GO TO 4104 

00040900 



C»»>AUTO FIlYfN5»l»4RT 1 

00841880 

TSN 

0172 

410 

PHjm = PHIG-PHI 

00041100 

ISN 

0173 


THI(l) = THETAG-TUFTA 

00041200 

ISN 

0174 


PSJIII s PSIG-PSI 

00041300 

ISN 

0175 


Alt 11 - Atffl$-AtPM 

00041409 

ISN NITS 



00061900 

tsw 


88041680 

ISN 

0178 


RICH = RG-R 

00041700 

ISN 

0179 


Bid) “ BETAG-BETA 

00041800 

ISN 

0180 


PHOm - CI39»4PHJm*PHSUM 

00041900 

ISN 

0181 


TMOIU > G138I*TH!(1>4THSUM 

30042000 

ISN 

0182 


PSCtll « 6(40I4PSJI1I«PSSUM 

80042100 

ISN 

0183 


AOfll ■ vlUtl IRAKI l♦ASON 

08<K??88 

ISN 

0184 

4104 

IF (IM) 412,411,412 

00042300 

ISN 

0185 

411 

UI m = UG-U 

00042400 

ISN 

0186 


UOd) - G( 5)4Uim^USUM 

00042500 

ISN 

0187 

412 

IP fISl 4l4t413«4l4 

00042600 

ISN 

0188 

413 

NIfU » HG«H 

00042700 

ISN 

0189 


MOm • Gf 6I4HII IKMSUR 

00047800 

ISN 

0190 

414 

POdl * Gm*9i{ij 

00042900 

ISN 

0191 


oom = GI2)*Q!m 

00043000 


9.14.3.2 DFCS (oont'd) 










PAGF 01 T 


ISN 0223 
ISN 0226 
ISN 0227 
ISN 02?« 
TSN 07?^ 
TSN 0230 
ISN 0231 
ISN 0232 
ISN 023 
ISN 0234 
TSN 0235 
ISN 0236 
ISN 0237 
ISN 0238 
ISN 0239 
TSN 0240 
ISN 0241 
ISN 0242 
ISN 0243 
ISN 0244 
ISN 024S 
TSN 0246 

^ ISN 0247 

^ ISN 0248 

TSN 0249 
TSN 0250 
ISN 0291 
ISN 0292 


ISN 0293 
ISN 0294 
TSN 0259 
ISN 0256 

ISN 0297 
ISN 0298 
ISN 0299 

ISN 0260 
ISN 0261 


9114 IF (IN) 513tS12*513 
912 OECU » Kan«U(?(l» 

SBU > K(12>4U0m 
91’ IF TIST 919,914,919 

514 DFCH * K(i3)*Hr'm*rLflrAS 
SBH « K(14)*HCn )>SPIAS 

515 DEC ■ CEC4CECU4CECH 
OSBC • SBU4SBH 

IF IN0DE-1> 917, 916, 917 
916 nSFC « nS9r4KT4T*nSFN 
5\y OELC » CFC-PAC 
DFPC « CFC4DAC 

WRITE (6,9500) CRC ,0E LC ,OERC ,0S BC, T NEXT 
GO TO I 91 8, 60), NODE 
518 GO TO (60, 919^919), NANMOO 
919 TRIMltl » • ATRTN 
TR1ML(2» - FTRIM 
T»'TML(3) » RUTRIH 
TRIMO(l) « PTRIF 
TKIM0(2) • QTMM 
TRIN0(3) • RTRIN 
GO TO 60 

c>»>>transition control 

520 GD TO 60 
C>>»>ENTRY CONTROL 
930 GO TO 60 
C»»>OR»IT TVC CONTROL 
540 GO TO 60 
C>>»>nPPIT PCS CONTROL 
590 GO TO 60 

C>>»>1NSERT10N TVC CONTROL 

960 GO TO 60 

C>»»BOOSTER TVC CONTROL 
970 «0 to 60 

C «*#*«*«««4**«4*4i4i**4i*«4i*4 4«i»«»«*4t****4t4r***iti*4i***4t4<4>»44i««'*****9i44* 

C»»>F11TFR ROUTINES-PART 2««< 

C ssaaas aaaaaaaa aaaeaa 

C»»>NANUAL FILTERS-RART 2 

6000 NSUM * C03>*N1 (1)46(39 )*N1(2)4G(32)9N0(U4G(34)«N0(2) 

N?T2) • Nt(l) 

NC(2) - NOm 
GO TO 61T 

C»»>AUTO FILTERS-PART 2 

610 PHSUR - 6(49)«RHJ(ll46(91)*RHJ(2)46(42)«RH0(tt46(48}4RH0(2) 

THSUM « 6(44)*THl(l)46(90)*THI(2)46(41)*THO(l)4G(47)*THO(2) 

PSSUN • 6(46)*PSJ(1 )46( 9?)*PSJ( 2»4G(43»9PSnn)4G(49)*PSrM7T 
ASUM » G(33)4Aim4G(39)*AI(2)4G(3?)*A0(ll4G(34)*A0{2) 

PHJ(?» * PHJ(1> 


00047900 

00048000 

00048100 

00040200 

00048100 

00048400 

00048500 

00048600 

00048T00 

00046800 

00048900 

00049000 

00049100 

00049200 

00049300 

00044400 

00049500 

00049600 

00049700 

00049800 

00049900 

06090000 

00050100 

00050200 

00050300 

OOOS0400 

00090900 

00090600 

00050700 

OOO5O0UC 

00050900 

00091000 

00091100 

00091200 

0005) 300 
00051400 
00051500 
00091600 
00091700 
00091800 
00051900 
00052000 
00052100 
00092200 
00092300 
00092400 
00052500 
00052600 




9.14.3.2 DFCS (cont'd) 















PAGF on 


ISN 

0301 

HG > H 

00057500 

ISN 

0302 

ELBIAS « 0. 

00037600 

ISN 

0303 

SBIAS - 0. 

00057700 

ISN 

030* 

GO TO 703 

00037800 



C>>>>>LANOING APPROACH TRAJECTORY FOR MANUAL MCnFS«<« 

00057900 



(* SSSSSCSS SSSSSSSSSS «SS S5SSSS &BS3B^ 

00058000 

ISN 

0305 

702 00 7022 1-1*5 

00053100 

ISN 

0306 

IF <RGO-RLANO< m 702**702**7022 

00058200 

ISN 

0307 

7022 CONTINUE 

00088300 

ISN 

0308 

I - 3 

00058*00 

ISN 

030P 

702* IP » I-l 

00058500 

ISN 

0310 

Rl » RGr-RLANC( IP J 

00058600 

ISN 

0311 

GO TO ( 703, 7C**7C3*70*) *IR 

00058700 



C»»>F LAKE (18-11 AND PULLUPt 18-31 

00058800 

ISN 

0312 

703 GEA8 ■ 1 

00058900 

ISN 

0313 

IP n«-2» 7032,7032,703* 

00059000 

ISN 

031* 

7032 J » 1 

00059100 

ISN 

031 5 

GC TC 7036 

00059200 

ISN 

0316 

703* J - 2 

U0059300 

ISN 

0317 

7036 8L2 » iL«Kl 

ooqit*qo 

ISN 

0318 

UG « ULAN0(I8»*USL0( 181*81 

ooqqqioo 

ISN 

0319 

MG • PHM* j»*PH(?,J»*8t^PHn,J»*RL?*PMI*,JI*8l?*8t 

00089800 

ISN 

0320 

GC TC 70*6 

00059700 



C»»>APPROACH GLIDE! IR-*I ANC PRE-FLARE GLIOFMR-?! 

00059800 

ISN 

0321 

704 IF (18-3) 70*2,7042,704* 

90059900 

ISN 

0322 

70*2 1-2 

oqqqq^qo 

ISN 

0323 

UG • UtAN0f2)*USL0(3)*8L 

00080100 

ISN 

03?* 

6EA8 • 1 

0O080P00 

ISN 

0325 

GO TC 70*6 

00060300 

ISN 

0326 

7044 J * •» 

00060*00 

ISN 

0327 

. UC « UtAN0(4)*EXP(USI.C(4)*RU 

00060500 

ISN 

0328 

N6 « HiAN0tf»MSL0(J)*8L 

0OM080O 

ISN 

0329 

fiA8 * q 

00080700 

fSN 

i>330 

70*6 flNTAS *’ WCfIM 

09080800 

ISN 

0331 

S8IAS - SB(IR) 

00060900 

ISN 

0332 

Gn TC 705 

D0061000 



c *4i«**«*«»»*«****«***4i*****4i*t*****4' **************** 

,**«**«,!*]»*« «««0006 1100 



C»>»PA8ANETE8 UPDATE 80UTINES<<<<< 

00061200 



C •..»••••• ««««.« mmmmmmmm 

00081300 

ISN 

0333 

T031 ANPAN • AMPAN*.05*( ALPM-ANEAN) 

00081*00 

ISN 

033* 

UHEAN * UMFAN+.05-IU-UMFANI 

00061500 

ISN 

0335 

IF (IS) 80,7052,00 

00061600 

ISN 

0336 

7052 PTEMP - RHO*EXP(-SCALE*H 

00061700 

ISN 

0337 

8H0V « P7ENP*UHEAN 

00061800 

ISN 

0338 

DYN8 • .5*8H0V*UMEAN 

00081900 

ISN 

0339 

60 TO 7033 

00082000 



C»»>PARAHETER FSTIMATION FOR MANUAL MODES 

00062100 

ISN 

03*0 

706 CALL TABLE(PC0MP,AMFAN,PTFMP,4) 

00062200 


2 " 


9.14.3.2 DFCo (oont'd) 
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'•T*T* 

Mi- % 






ISN 0341 
1SI| «t42 
tSU 0143 
TSN 0%44 
ISN 0145 
ISN 0146 
ISN 0347 
fill os4i 

ISN 0151 
ISN 0152 


i’^P 

GM4J » -PTFMP4FILI2) 

CALL TABLEIPFIT,AMEAN,PTEPP,3) 


00042300 
0Mii400 
09042900 
0004 2400 

00062700 

00062800 


# 9^92 »1 1 




Gill » PTEMP 
GO TO 712 


00063100 

00061400 


.yia NODE 


ISN 0362 
ISN 0363 


PTFNP4PTEHP1 
GI39»4(FIL( ll-l.l 


00064500 

00064600 


ISN 0368 
ISN 0169 


Kt8> « PTENP 

712 CALL TABLE! eFITtDYNP*PTENP,4) 


00065100 

00065200 


C»»>PARANETER ESTIMATION FOR ENTRY MODE 
730 GO TO BO 


00065700 

00065600 


ISN 0373 


C»»>PARAMETER ESTIMATION FOR INSERTION TVC MCOE 
760 GO TO 80 


00066300 

00066400 


ISN 0376 


00066900 

00067000 


ISN 0378 


81 I START « 0 





PAGF 015 


ISN 

037S 

ISN 

0380 

I SN 

0381 

ISN 

0387 

ISN 

0387 

ISN 

0384 

ISN 

0385 


IS « IS^l 

C»»>INC«EMENT INDICES 
82 IM « I»<*1 
TS - TS»1 
IM = M0C( IHt TTM) 

IS = MrcnSf ITS) 

C»»>RFTUPN TO GNtC EXECUTIVE<<<<< 
RETURN 
END 


00067100 

O00672D0 

00067300 

00067600 

00067500 

00067600 

00067700 

00067800 

00067900 


3 " 


2 " 


9.14.3.2 ^|CS (cont'd) 



LEVEL ?0.1 fMAY 71» 


ns/760 FORTRAN H 


OATF 71 ,70?/13.2'5.56 





COMPILER CPTIOMS » MAME< MAlN*42PT«02tL lNECNTs50» S lZE»OOOOK» 

C>»»irM0A» tWTfflPCl»?Y«fl<«« N ' 

C>»>>A(I , J) sTABLF, Z*INDEPENDENT VARTARLE* R*CUTPUT 
C N*NUPBER OF TAPULATEO PAIRS 

OIMENSlCMAIZfll 

DO .V' 


OOOAtOOO 

OOOAtIfM) 

00068200 

00068300 

0006A400 


ISN 0003 

ISN OdOA 
1SH M^t 
YS«f 99 ^ 

fSN 0007 
ISN 0008 
X$M 0009 

tsu 0910 
Iff 

ow? 

ISN 0013 
ISN 0014 


00060800 

00068900 




RETURN 

END 


00069400 

00069500 






>.S't 



'^-Trrr”' 7 

•V ' ■ 

:’V. , • . 

4 




' 



Vv 

* ' ■ 






EKI A ft i ^ ' 






jfr , . V » 










8t-Vl*6 


n/ri/71 


AUrjFLUi* CHART SET 


RAIN 


PAGE 01 


CFART TITLE - IN TRODLC TOR V CCRRENTS 


>>>>>SSV UNIFIED DIGITAL FLIGHT CCNTPCL SYST tP-WEBS ION AD 10-2<7-71 
F.F. STENGEL PIT C SDL 

4 «« 4 *:««***«** 4 «»««**»*«****«****#*«* 4 '**«*'«*»« *»«*••** * 41 * «**«»**»•* 


AM I I 
ALPH 
ALPHG 
AMEAN 
A0( I » 

A SUM 
ATRIM 
BETA 
BE TAG 
eFITII.J I 
Bill) 

SOI I ) 

esuM 

DAC 

CACINT 

CAM 

DEC 

CECH 

CEC INT 

OECU 

CEL 


NCMENCIATLFE 

«s3xss.%ssrssx 

ANGLE OF ATTACK (FILTEB INPUT I 

" " " IMEASUREMENT) 

" « " IGDIOANCE) 

MEAN ANGLE CF ATTACK ECR PARAM. EST . 

" " « I filter OUTPUT) 

• " » ITEMP) 

•AILEFCN* TRIM DISCRETE 
SIDESLIP ANGLE IMEASUREMENT) 

" • IGUICANCE) 

TABLE OF SIC=SLIP GAINS 

« " (FILTER INPUT) 

'• " IFILTER OUTPUT) 

'• " (TEMP) 

•ailercn* angle cc)*manc output 

KCAH integral in 'AILERON* ANGLE COMMAND 

•AILERON' ANGLE CCPMAND (MANUAL) 

'ELEVATCH* ANGLE CC)»PANO OUTPUT 
'ELEVATJK* CCMMAND DUE TO ALTITUDE FEEDBACK 
RCAH integral in 'ELEVA-'Cr' ANGLE CCPPAND 
•ELEVATCR' COMMANC DUE TO VELOCITY FEEDBACK 
LEFT FLEVON ANGLE I MF A SURE MEN T I 


9.14.3.2 DFCS (cont'd) 





CL_ 


DELC 

CELG 

CEM 

OER 

CERC 

DERG 

OR 

CRC 

ORG 

CRM 

OSB 

OSBC 

CSBG 

OSBN 

CVNP 

ELBIAS 

EL(I} 

ETRIM 

FAIL(l) 

FASTII) 

FH( ItJ I 

FIUI ) 

FLAGl 

FLAG2 


H 

M 

•• (CCMMANC) 

tl 

It 

•• (GUIDANCci 

•ELEVATOR* ANGLE COMMAND IMANUALI 

RIGHT 

ELEVCh 

angle (MEASUREMENT I 

N 

fi 

•* (COMMAND! 

M 

M 

•• (GUICANCEI 

RUCCER 

iHGLE 

(MEASUREMENT! 

M 

m 

( COMMAND) 

M 

m 

(GUICANCE) 

RUCOER 

ANGLE 

COMMAND (MANUAL! 

SPEED 

BRAKE 

SETTING (MEASUREMENT) 

M 

M 

•• (COMMAND) 

H 

m 

•• (GUIDANCE) 

SPEED 

BRAKE 

COMMAND (MANUAL) 


DYNAMIC PRESSURE ESTIMATE 

•ELEVATOR* BIAS TO TRIP SPEED BRAKE BIAS 

TABLE CFRELEVATORLeiASESATO TRIM SPE*€0 BRAKE 

•ELEVATOR* 1RIM DISCRETE 

FAILURE DISCRETES 

FAST SAMPLING INTERVAL TABLE 

CONSTANTS FOR FLARE & PULLUP ALTITUDE REFERENCE 
FILTER CONSTANTS USEC IN PARAMETER ESTIMATION 
NEM START CISCRETE I EXTERNAL ):0-N0RNAL CYCLE 

l«NFli START 

FLIGHT MODE FLAG ( EXT ERN AL I : 1«M ANU AL MODE 


2=AU TOMA TIC MODE 



9.14-50 


n/Ql /11 

CKPT TITLE - INTROCUCTaRY CCM.-IENTS 


4UTQFLQW CHART SET - MAIN 


PAGE C? 


F LAG 3 

FLAGA 

FLAGS 

GUI 

CEAR 

GFIXI I.J) 
F 

FG 

Him 
FL.HLP 
FLANOI 1 I 
HOI 1 I 
FSUM 
I 

TM 


3=TRANSITICN MODE 
A^ENTRY MODE 
5=0R8IT TVC MODE 
6=ORBIT PCS MODE 
7=INSERTION TVC 
8=S0CSTEP TVC MODE 
PANUAL MODE FLAG I EXT E . N AL » : 1=0 CRECT 

2»RATE CCPMAND 
3»RCAH 

FARAHETEB ESTIMATION FLaF ( EXTERNAL I sO» SK IP 

1»00 

MANUAL PCUE AUTO. SPEEC BRAKE FLAGIE TERNAL): 
0«NC AUTO, SPEED BRAKE 
I>ALTC. SPEED BRAKE 
FILTER GAINS FOR CURRENT FLIGHT MODE 
LANDING GEAR DISCRETE 

FILTER GAIN TABLE FOR ALL FLIGHT MODES 
ALTITUDE ABCVE RUNMAV IME ASUREMENT ) 

H ti (GUIDANCE) 

ALTITUDE ABCVE RUNWAY (FILTER INPUT I 
TEMPORARY STORAGE IN REFERENCE TRAJECTORY 
ALTITUDE PCiNTS ON LANDING APPROACH 
ALTITUDE ABCVE RUNWAY (FILTER OUTPUT) 

" n » (TEMP) 

GENERAL PUFFCSE INDEX 
MEDIUM SAMFl ING RATE INDEX 


9.14.3.2 (oont'd) 



is-vre 


IR 

IS 

I ST ART 

ITM 

ITS 

ITURN 

J 

K( I ) 

KFIX( If J I 
FANMOO 
FEOIUHI I) 
MODE 
M( II 
NO( I) 

NSUN 

F 

PCONPI If Jl 
PFITd f Jl 
PG 

PHFiri I f Jl 

PHI 

PHIG 

PHJI I I 

PHOII I 

PHSUM 

PlilS 


REFERENCE IRAJECTORY INDEX 
SLCN SAPPLIAC RATE ISCEX 
STARTING CYCLE DISCRETE 
A OF FAST CYCLES PER MEDIUM CYCLE 
# CF FAST CYCLES PER SLOW CYCLE 
TURN-ON CYCLE DISCRETE 
REFERENCE TRAJECTCBY INDEX 
CCNTFOL GAINS FOR CURRENT FLIGHT MODE 
CONTROL GAIN TABLE FCR ALL FLIGHT MCCES 
MANUAL MCOE FLA( (INTERNAL I 
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■«»- — ■*">1 


f 

9»IAS - 0. 

1 

1 


* • 

1 

1 

2140 1 NOTE 

19 

1 


I 

1 


* • 

1 

1 • 

*«•••••• 

• 0 

I 

PMSIHi - 0. 


1 


♦ • I-/0I 

I 

1 6 

BEGIN 00 LOOP 

0 

1 


1 

1 


* ELAG5 • 

4 

1 « 

2142 I « 1. 3 

0 

1 

1MSU6 - 0. 

1 

1 


* * 


I * ' 


• • 

1 


1 

1 


* « 


1 

1 


! 

PSSUM « 0. 

{ 



• * 


1 

|< 


1 


1 

1 


• 


1 

1 

20 

1 

A5VM ■ 0. 

1 

1 


1 l*l 


I *- 


„.^0 




1 


1 


1 1 

TRIMLIII « 0. 

1 


1 


1 


1 


1 


♦ 


1 


1 


1 


1 

1 



1 

07 

1 


1 


1 

1 


♦— 

^ ^ ^ ■ 


1 


1 


1 

2142 1 

21 

1 

■SUM » 0. 

1 

1 


2132 i NOTE 14 


1 *** 


• 

1 


1 

\ 

* 

«•**«**•«« 


1 1 

TRINOIII - 0. 

1 

1 

HSW - 0. 

1 


0 

EEGIN CO LOOP * 


1 ** 

■ 


1 


1 

1 

0 

2134 I - 1, 5 * 


1 

1 


1 

rsuf « 0. 

1 

i 

0 



J 

1 


1 



1 


• 


1 

1 


1 

I 

CSii7 - 0. 

1 

1 

1 


1 

1 

• 22 


1 

1 

«SUN > 0. 

1 

1 

1 





0 0 


1 


1 

1 





* 6 

NO 

1 

asuM > 0. 

1 

1 





• END OF DO 

« — 




1 





• LOOP7 • 



1 


1 





• • 



1 


1 





0 0 





-4 





0 



• « 
• 6 * 21 * 


>»»0IIBIT RCS NODE 

06,19 >• 

2$0 I 26 


• • 

. 6 . 21 . 

. . 

... 40 


»»>INSEItT10N TVC 
MODE 

C6.19 >• 

260 I 27 


. . 

. 6 . 21 . 


40 


»>»BOOSTER TVC MODE 

06,19 >* 

270 I 28 


. 6 . 21 . 


40 


».. 


40 



9 ^ 1*6 




CKFT I ITUt ' PROCtDURES 


J»>>F ILTER 
FCtTIRES-F4«T 1«<« 


)>»>PANUAL 
MLTEBS-P#RT 1 
ACC 7J«iR CISCRETbS TC 
MANUAL CCHMANOS 

Ct.21 >* 


ACC I 01 


I ROtll - I 

I COMANIIll *■ I 

I NSUM I 

I I 

I ATRIN >= I 

I M15)*A1RIH I 

I I 

! ETRIM « I 

> lcfl6l*ETRIM I 

I I 

I ftUTRlR - I 

I R(lTI*Bt1RIM I 

* — « 

I 

I 



1 


I o 
• IV 

1 

PTRIH « 


1 

} 

TP2*(ATRIM 

♦ 

1 

1 

IRIMLUn 


1 

1 

TRlMOIll 


1 

1 

1 

1 

QTRIM > 


1 

1 

1 

TF2*(ETRIM 

4> 

1 

i 

IRIMLiai) 

♦ 

1 

1 

TRINOI 21 


1 




I 

I 

I 03 

I RTRIM » I 

I IF2*1HLIRIM ♦ I 

I TRIHLOl) » > 

I TRIRCI3I I 


»>>>ALT0 
FILTERS-PART 1 

oe , 21 >* 



410 1 

07 




1 

PHJI 1) « PHIG - 

1 

1 

PHI 

1 

1 

1 

THU 11 = TMETAG 

1 

1 

THETA 

1 

1 

1 

1 

PSJUI * PSIG - 

1 

1 

1 

PSI 

1 

1 

1 

1 

AlCi « ALPHG - 

1 

1 

1 

ALPH 

1 


1 



1 

1 

OB 

1 

pim » PG - p 

1 

f 

1 

1 

1 

CKII « QG - C 

\ 

1 

I 

1 

1 

Rim • RG - R 

1 

1 

1 

1 

Sill) - BETAG - 

1 

1 

1 

BETA 

1 


1 



1 

1 

09 




1 

PHCCl) = 

1 

1 

G(3St*PHJ(l) ♦ 

1 

1 

PHSUM 

1 

1 

1 

1 

THOU) * 

1 

1 

1 

GI3Ri4THl(ll ♦ 

1 

1 

1 

THSUM 

1 

1 

I 

1 

psom = 

1 

1 

GI40)*PSJI1) ♦ 

1 

1 

PS SUM 

1 


I 


i-ly C- jf .T SET 


PAIN 


PAGE 10 


>»»TPANSI TION 
FILTERS-PART I 

06.21 >• 

A20 I 20 

• • • 

• • 

. T.Ol. 

• • 

• •• 30 





D . 


>« 

413 I 14 




- 


1 

Him - H6 - 

H 1 


1 


1 


1 

Hom > 

1 

>»»ENTRY 

1 

Gi6)*Him * 

HSUM 1 

FILTERS-PART 1 

♦- 





1 


06,21 >• 

_ 

>1 


430 1 21 


414 1 

15 

• • • 

1 

POIll • 

1 

• • 
. 7.01. 

1 

Gf 1I*PII 11 

1 

• ¥ 

1 


! 

... 50 

1 

Qom * 

1 


I 

6121*017 11 

1 

1 


1 

1 

Rom 

1 

1 


1 

1 

6I3I4RII 1) 

1 

1 



dom >= 

G(4)*BIUI 


* * 

I 

I 

I 16 »»>ORBIT rvC 

• 4 r-ILTERS-PART 1 

I COMPUTED GC TC J 

I FOR MODE I 06,21 >• 

* * 440 I 22 

I 416 10.17 I ... 




50 



»»>0KBIT lies 
FILTERS-PAKT 1 

06.21 >• 

450 I 25 

• •• 

r . 

• 7.01. 

• • 
... 50 


»»> INSERT ION TVe 
riLTERS-RART 1 

06.21 >6 

460 I 24 


. . 

. 7.01. 

. . 


. . . 


50 



wrii'b 
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CHART TITLE - PKUCfOURES 


/ 5114 / 


>J>>>eOOSTtH TVC 
FILTERS>PART 1 

C6.il 

47C I Cl 

... 


... 5 « 


}»»CCNTRCL 

«cwrucs««( 


»»>01KECt MANUAL 
CCNTML 

o.oa — >• 

5000 I 

I ORC • Km*0RM 

I 


»»>ACT0 CCNTRCL 

OT.Ol— >* 

51C I 34 

I CAC - I 

I K(3t«PHlUl) * I 

I KI5>*PCUI « I 

I K(15I*PSU(1> * I 

i K(6 I«RnASH I 

I I 

I DEC - I 

I K(2»*TH0m * I 

I M9)*CC(ll * I 

I K(4)*A0m I 

I 

I 05 

I ORC - I 

I Km*PHom ♦ I 

I KI7)*PC(ll 4 I 

I Kll6l*PS0m « I 

I K(8)«RbASH * I 

I K( 10)600(11 I 

I 

I 

I 

»»>MtLTI-RATE SPEED 
BRAKE CONTROL 
I 

511 I 06 


« KI5I6CAM 

1 i 
1 1 

1 CCMPUTEO CO TO 1 

1 

1 

1 

1 


1 

1 


1 i 

1 FOR MODE 1 

1 

1 


1 

- K(4)«0SBN 

1 1 


1 

1 


1 


-* 1 

1 5111 11.C7 1 

1 

1 


1 

1 

1 

1 

1 

1 5114 11.09 1 

1 

1 

1 

1 


914 1 

1 

1 

1 

1 

1 

1 

( 

OECH » 

• •• 

1 

1 

1 

1 

1 

k(15I*HU(1I ♦ 

• • 

1 

1 

1 

1 

1 

EL6IAS 

.11.16. 

1 

IF OUTSIDE the RANGE 

1 

1 

1 


• • 

1 

1 

1 

1 

1 

S8H - 

... 517 

1 

11. C6 >1 

1 

1 

1 

K(14l«K3m ♦ 


11.066— X 

I 

• l«9 

* « 

« * 

(-/♦I • * 

* IM • 

« * 

* * 

* • 

* 

I (») 

I 

I 

I 

I 

I 

912 I 10 

I OECU - I 

I K(11I*U0(1) I 

I I 

I SOL ■ K(12)*U0(1) I 

I 

>i 

I 

915 • 11 

* * 


(-/♦) * 
♦ — • 


* 

• * 

* 

I (0) 



o 


}»»SAS CD KCAH 
)*AAU«l CCATRCL 

ci.qa*— >• 

»0I0 I 03 


I 

I OAC - radi ♦ I I 

I K(5I«P ♦ I I 

I K4«I«IMSH I I 

I m - «c(ii ♦ I I 

I M«I*C ♦ I I 

I 1 I 

I I i 

I G«C - Mi«SM ♦ I . I 

I HITI«P « I I 

I KllflOtCIll « ORC I i 

» * I 

I I 

I I 


i i 

5111 * 07 I 

* * I 

* * I 

* * (♦! I 

• FLAGS • ♦ 

« « 

* • 

* 0 
* 

I 

I 

I 

I 

I 

5112 I OA 

I OSBC • KtAI*OSftM I 
- 

I 


♦— 


SBIAS 


I I 

I • 

I I 

♦ >1 

515 I 


13 

I DEC « DEC ♦ I 

I OECU ♦ OECH i 

I I 

I csec « SBU ♦ SBh I 

I 

I 

I 

• lA 

« * 

• * 

(-/♦!• • 

♦ * MOOC - I * 

* * 

« • 

* • 

* 


S16 


I 


I 


IBI 


I 


I 


IS 

«•* 

I 


CSIC > OSIC ♦ 
HfA»*OBBN t 

i 

1U02*— >1 

S17 I 1* 

« 

I OELC • OEC > OAC i 
I I 

I OERC > OEC 0 OAC t 

I 

I 

17 

i WRITE TO OEV / 
/ 6 / 

/ VIA FORMAT / 
t 9500 / 

/ FRCM THE LIST ! 


I 




>• 

I NOTE IB 
»•«**••*•*« 

• LIST - ORCt OELCt * 

• DERCt CSBCf TNEXT • 

• *«*•**«*** 

I 

I 

I 19 

* — _______ • 

I COMPUTED 60 TO I 
I FOR NODE I 

I SIS 11.20 I 

I BO 7.0S I 

• 0 

\ 

I 

I 

IF OUTSIDE THE RAN6E 
I 

11.19 >1 

SIB I -20 

0 — » 

I COMPUTED 60 TO | 
I FOR NANNOO | 

♦ 

•I BO 7.0S I 

I 919 11.21 I 

I 919 11.21 I 

- 0 

I 

I 

I 

IF OUTSIDE THE RAN6E 
I 

11.20B-- >1 

919 I 21 

• 0 

I TRINLIll - ATRIN I 
I I 

I TRINLI2I - ETRIN I 
I I 

I TRINLI3> - RUTRIN | 

I I 

I TRINOm - PTRIN I 
• • 

I 

I 

I 22 

* — ♦ 

I TR1M0I2I > QTRIN | 

I I 

I TRIM0I3I « RTRIN | 

I . 

I 

I 

• • • 

• • 

• 7,03. 


60 



OL'^fS 
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CE/PT TITLE - PPCCEDOOFS 


»»>THANS1T1UN 

»>»ilOCSTER TVC 


»>»AUTO 


CCKTRCL 

CUNTKOL 


FILTERS-PART 2 


'•7.TI >• 

(7,01 >* 


07,03 >* 


S20 1 01 

S70 1 36 


610 1 

OB 

• •• 

• • • 


— 

— 4 

• • 

• • 


1 PhSUR ■ 

1 

. 7.33. 

. 7.J3. 


1 C(4SI*PHJ(1) ♦ 

1 

• • 

• • 


1 CI51)«PHJI2) ♦ 

1 

• • • 60 

,,, 60 


1 G(42>*PH0(U * 

• 




1 G« 60)*PH0f 2» 

1 

1 




1 

1 

1 

09 




i THSUR • 

1 




1 G(44)*THim ♦ 

1 




1 6(50»*TNI(2) ♦ 

1 




1 G(41I*TH0(U ♦ 

1 




1 GI471*TH0(2I 

1 

»»>ENTRY CCITTROL 

>»»EILTER 





ROLTINES-PfRT 2««< 

1 


07.01 >• 



1 


530 1 02 



1 

10 

• •• 




— * 

• • 

»»>MANUAl 


1 PSSUN - 

1 

• 7.03. 

FILTERS-PART 2 


1 G(46I*PSJ(1I * 

1 

• • 



1 G(52»*PSJ(21 ♦ 

1 

• •• 60 

C7,0A*~>* 


i G(43(6PS0(U . 

1 


6000 1 

07 

1 6(49I*PS0(2I 

1 




1 

1 


1 NSUM • 

1 

1 ASUM - 

1 


1 G(33l*Kim * 

1 

1 G( 33»*AM II 4 

1 


1 G(35I«M(2I ♦ 

1 

1 G(35I*AI(2> 4 

1 


1 G(32l*NOm * 

1 

1 C(32l*A0m 4 

1 


1 G(3AI*6C(2) 

1 

1 G(34l*AO<2l 

1 


1 

1 




»»>CRBIT TVC 
CONTRUL 

C7.31 >• 

SAC I C3 


NII2) • Nim 
NU(2I > Ncm 


PHJ(2» » PHJ< I 1 
PHO(2l - PHOm 
THH2» ■ THKU 





. 7.03. 

• • 
... 60 


>»»ORBIT RCS I 

CCKtRCL ♦- 

07.01 >• 

5S0 I 04 

• . 

• 7.03. 

• • 

• a. 60 


»>»INSiR1ICN TVC 
CONmOL 

C7.0I >• 

S60 I OS 


• 7.03. 

. . 

... 60 


I 


I TH0«2I - THO<n I 



/i3.ni 



e6-tr6 


n/:i/7i 


AUTUFLOi* CHAM SET 


MAIN 


PAGE 13 


CMRl TITLE - PROCEDURES 


12 . 13 *— >* 

I 

61 ’ • 01 

« * 

* * 

I-/4I • * 

♦ * IS * 

I • • 

I • * 

I • • 

I • 

4 I (01 


vO 

■r^ 

V*» 

NJ 


»»>TR4NSITION 
FILTERS-PART 2 


»»>BOOSTER TVC 
FILTERS-PART 2 


C7.03 >• 

620 I 10 

• • 

• 7.03. 

. . 

... 70 


07.03 >* 

670 I 15 


c+ 


7.05. 


... 70 


L 

I 

1 


L 

I 

I. 

I 



614 1 

02 


HSLP - 

1 


GIIEMMIUI ♦ 

1 


6I30MPI42I * 

1 


6(12f4N0(U ♦ 

1 


6(24I«NCI2I 

1 

1 


HK2I • Rtd) 

1 

1 

1 


NCC2I « 1-cm 

1 

1 

4- 




_ I 

615 I 03 


I CCMFU 1 E 0 60 TC I 

i FOR NCOE I 

• — — « 


I 61C 13.04 I 

I 617 13.C5 I 


I 

I 

I 

IF CL1SI0E TFE RANGE 

I 

13.03 >1 

616 I 04 

I CCHPUTEO GO TO I 
I FOR PAFMOO I 


>^»>ENTfir 
FILTERS-PART 2 

07.03 >• 

630 I 11 

... 

. . 

. 7.05. 

. . 

... 70 


»»>OReiT TVC 
FILTERS-PART 2 

07.03 >• 

640 I 12 





0 



a-^i*6 


I TO 
I 618 
! 617 


7.P5, 


7,05 I 
13.09 I 
13.C5 I 



1 

* 


1 



1 


IF 

CCTSIOE RANGE 

13 

.02*— j( 

617 1 

05 

« 

1 

PS UP » 

1 

1 

1 

1 

G(13)*fllli ♦ 
G(2SI*PI(2I A 
CI7l*FCm * 

1 

1 

1 

1 

1 

6(19|6PC(2I 

1 

1 

1 

1 

CSUP - 

I 

1 

1 

1 

1 

GdAIPCKll ♦ 
GI26l*GM2t ♦ 
G(8I*€0I1I * 

1 

1 

1 


I 6<20|6«ei2t i 
I 



1 



1 

06 

1 

RSUP - 

1 

1 

G(1SI*RM1I ♦ 

1 

• 

G(2TMBI(2I * 

1 

1 

6<9I*RC(1I ♦ 

1 

1 

CI21t«RC(2l 

1 

1 


1 

1 

PII21 - Fim 


1 


1 

1 

FCI21 - fCdl 

1 


1 



1 



1 

07 


01(21 > Cl (11 

1 

1 


CC(2I > COdI 

1 

1 

1 


RK2I • RKII 

1 

1 

•> 

RC(21 « POdI 

1 

1 


I 

I 


->• 


I OB 


BSU»* - 

G(16)*Bim * 
C( 28 l*ei( 2 ) ♦ 
6(ioi*aoiii « 

GI22)*BOI2t 
81(21 - Bim 

BCI 2 I - eom 


I 

13.C6 >1 

618 I 79 


1 Q6LAST > QMASH 

1 

1 

1 

1 RhLASr > RM6SH 

* 

1 


I 

I 

I 


• •• 

I 7 . 05 ^ 


»»>ORBIT RCS 
FILTERS-PAOT 2 

C7.03 >• 

650 I 13 


• • 

* 7, 05, 

• • 

• •• TO 


»»>lNSERTtON TVC 
FILTERS-PART 2 

07,03 >• 

660 I lA 

• •• 

• • 

. 7.05. 

. . 

... 70 


. .. 


70 






n/ri/n 


4LT0FLCW CHART SFT 




PAGE lA 


Cl-flf»T TITLt - PSOCruuRCS 


>?»>L4M;ikg afprcach 

TPAJtCTGKy FD< MANUAL 
KCES««< 




L 


C7.C7 >• 

7C2 I NCTE Tl 

A********** 

* EEGIK CC LOUP * 

* 7C22 I « It 5 • 

**«*•*•«*•« 

I 

|< ♦ 

I I 

* 02 I 

* * I 

* * I 

« I 

BGC - * J 

RLAACin • I 

* * I 

* # I 

• I 

I !♦» I 


i 


7022 ♦ T3 I 

• • I 

• * I 

• • NU I 

• EKD CF OC ♦ ♦ 

• COOP2 * 

* ' * 

* « 

« 

Irfs 


t-/c» ♦ 
« 


'O 


>>>>>APPR04CH 
GLIOE(IR>AI AND 
PkE-FLAKE GlI0E(|R«2> 


14.06* — >* 

I 

7CA * 

* * 

* 

* 

• IB - 3 

« 

* 

* * 

* 

I 1-/C I 


7042 I 13 

♦ — — — — « 

I J - 2 i 

I I 

I UC - ULAN0I2I ♦ I * 

I USL0(3I*BL I I 
I I I 

I GEAR >1 II 

I « 

I I 



>• 


7044 I 14 

* — * 


J » 3 
UG > 

ULAN0(4t 

*EX-»|USLQ(4I*RLI 

HG - HLANOC |RI ♦ 
SLOfJ)*RL 

GEAR « 0 


I 

-14,10*— > I 

7046 I 15 

I ELBIAS > ELItR) I 

I I 

I SBIAS > SBCIRI I 


I 

I 

I 


*>■ 

Co 



c+ 

a 


. 7,09, 


705 


04 


* -♦ 

I 1-5 I 



w-ri.*6 


i 


I 

7024 I 95 

t IR * I - I I 

I I 

I RL - RGC - I 

I RLANOCIRI I 

*■ ♦ 

I 

I 

I '»6 

• — — • 

I CCNPLIEO CO TO I 

I FOR IR I 

I 703 14.C7 i I 

I 704 14.12 I I 

I 7C3 14.C7 I I 

I 704 14.12 I I 

» I ♦ 

I it 

I I I 


IF OCTStOe TFE RANGE 

I 

I 

I 

»»>FLARS(IF>1I ARC 


PU.LUFIIR»3I 

I 

14.0C*— >1 
703 I 07 

I GEAR - 1 r 

4 4 

I 

1 

» 


* 08 I 

• • I 

• • I 

« « l«l I 

♦ IR - 2 • ♦ 

« * 

• » 


l>/OI 



>« 

7034 I 

I J « 2 

I 

> I 

7036 I 

I RL2 - RL*RL 

I 

I U6 » ULANCIIRI * I 

I ISLOI1RI4RL I 

I I 

I l-C « FHIltJI « I 

I FH(2tJI«RL 4 I 

I Fh<3«JI*RL2 4 I 

I FH(4> J|4R12«RL I 

4i * 

I 

I 

I 



Il4.isl 
« • 

... 7046 


>* 

7032 I 11 

• — — - — • 

i J - 1 I 


I 

I 
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aUTOFLCw CHAKT SE^ 


MAIN 


CMfl TITtf - PRfiCEDJKCS 


; 3 >>>PAKA«tTE« UPCATE 
SCCT 1 AES<<<<< 


C 7 .CS )* 

7 C 51 I '•! 


I APEAK > APEAN * I 

I .C5*1ALFH - I 

I AMEAM I 

I I 

1 CPE AN « CPFAN « | 

I ,rS*lU - CMCANI I 

* » 

I 

I 

I 

• C2 

* « 

• * 

* * «-/♦» 

• IS *-* 

* * I 

* • I 

• • I 

• .... 

I ni . 7 • 

I . 12 • 

I •• • • 

I ec 


7052 I 73 


I PTEMP « I 

I PHC«EXFf>S€ALE*l-l I 

I I 

I RHCV > I 

I ETEPFPCPEAN I 

I I 

I CVNP * I 

I .‘PPHCVPINEAN I 

, - 

I 

I 


. • 

. 7.ir. 


»»>PAkAPETEK 
ESTIMATICN FUR MANUAL 
MODES 


07, ll* — >♦ 

706 I 

« « 

I 1 ! H 

18 I TABLE H 

I. I IPCOMP,ANEAN> H 
|0 I PTEMPtAl H 

II I H 

* — — — « 

I 

I 

I cs 

• — * 

I GUI * PTEMP I 

I I 

I GI13I - - I 

I PfEMPAFILIl) I 

« 

I 

I 

I 06 

• « 

III H 

IS I TABLE H 

I. I lUCOPP.OTAP. H 
|n I PTEMP, 31 H 

II I H 

0 * 

I 

I 

I C7 

I GI2I » PTEMP I 

I I 

I GIlAl » - I 

I PTEMP*FIL12> I 

« 

I 

I 

I 18 

III H 

|8 I TABLE H 

I. I (PFIT.AMEAN. H 


»»>PA«AMFTER 
ES 11 PATICN FOR AUTC 
MUCE 


C7. IC >• 

710 I 14 
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9.15 


APPROACH AND LANDING 


Approach cammences with the vehicle In the vicinity of the 
landing site and Includes the manuevers performed in the process of 
achieving the final approach trajectozy in the final approach plane. 

Final Approach terminates at touchdown on the runw^. Landing begins 
at touchdown and continues until the craft reaches ground taxi speed. 

The GN&C software functions to be performed during this phase are: 

1. Approach navigation using measurements frcm the IMD and the 
ground transponders. Both the shuttle state and the landing 
site stahe may be estimated. 

2. Approach and lAndlng guidance resulting in commands to the 

autopilot (attitude), thrust throttle commands, braking and/or 
lift flap commands. 

3. Autopilot computations to c<mqpute control surface commands in 

order to achieve the command attitude. 

9.15.1 Navigation (TBD) 

9.15.2 Guidanme 

Suhmittals are included herein for Approach Guidance and 
Final Approach Guidance. Landing Guidance requirements have not yet been 
defined. 
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9.15.2,1 Approach Guidance 


SUMMARY 

A technique la presented for guiding an aircraft In aubaonlo flight 
(with or without power) from any bearing and heading with respect to a 
runway to a reference trajectory In the final approach plane (FAP), 
that contains the runway. This technique produces three phases: an 

Initial turn toward the FAP| wlngs-level flight to the FAF^ and a final 
turn Into the FAP, Each phase Is flown with a nearly constant angle- 
of-attack and Indicated airspeed. The guidance selects the most effl-< 
dent set out of four possible sets of turn direction combinations; 
l.e.y left Initial turn with right final turn, etc. The energy manage- 
ment parameter Is the remge of the FAP entry point from the runway. 

The Initial condition requirements or the end condition constraints 
imposed on the supersonic flight phase are that the vehicle be within 
the footprint of capability (of the order of 20-30 n ml radius from the 
landing site), and be at an energy state which will remain subsonic 
during the Initial terminal area turn which la made at maximum l/D and 
45 degrees bank. (The heading angle at this point need not be con- 
strained . ) 

This guidance technique has been successfully flown on a six- 
degree s-of-freedom hybrid simulation for a straight-wing orblter. 
Results from an all-dlgltal program for both a straight and delta-wing 
orblter are presented. 


INTRODUCTION 

The Shuttle Orblter Is presumed to have re-entered the atmosphere' 
and received a navigation update follo%fing the communications blackout 
period. A post-blackout" guidance technique which Is the subject of a 
subsequent study will have steered the vehicle to a subsonic state 
(figure 1) which Is considered the beginning of the terminal area 
mission phase. 

Since the Orblter, In Its operational configuration, will not have 
a cruise capability. It Is necessary to provide a terminal area guid- 
ance technique which makes most effective use of the available energy. 
Also requlr^ In this mission phase Is the neoeeelty to operate In 
Instrument weather conditions. Because of this requirement, the guid- 
ance must be designed for ease of pilot monitoring and/or control when 
visual pilotage cues are not available. 

The guidance technique presented herein will guide either the 
powered or unpovered vehicle from a high aubaonlo energy atate In the 
most efficient manner (within the conatralnta Impoaed) from any poaltlon 
and heading with reference to a runway that la within the vehicle foot- 
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print of capability, to a reference trajectory In final approach plane. 
From that state, the Orblter can complete an unpovered landing. 

The design Includes logic which makes It compatible with good air- 
craft Instrument flight procedure and results In maximizing the periods 
of vings-level flight at nearly constant airspeed and angle-of -attack. 

The objective of tnis internal note is to present the derivation of 
guidance equations. To demonstrate this guidance, some results of an 
all-digital program are presented for both a straight and delta-wing 
orblter. 


DISCUSSION 

To guide an aircraft from any state relative to a runway to a 
linear reference trajectory in the FAP containing the runway, a two- 
tum maneuver is determined as shown on figure 2a. A constant and 
equal radius is used for both the initial turn toward the FAP and the 
final turn into the FAP. There are four possible turn direction com- 
binations of two-tum maneuvers (figure 2b). At any guidance computa- 
tion cycle, the guidance selects from these four sets that which will 
place the vehicle in the FAP with the greatest energy. The FAP entry 
point is then adjusted until this energy is the amount required to 
Intercept the FAP reference trajectory once the vehicle is in FAP. 

The two main parts of the guidance are determining the reference 
trajectory to the FAP, which is continually computed, and the attitude 
and attitude-rate connands to guide to this reference trajectory. 
Attitude commands might be processed more often than the reference tra- 
jectory determination equations. Simulations to date have used equal 
computation time intervals of two seconds. 

Determine Reference Trajectory 

Determine ground tracks . - With the initial and the final turn 
maneuvers of the same radius, it is a fairly simple problem of geometry 
to solve for the four possible ground tracks of figure 2b. The vector 
equations to be solved are shown in appendix A. The answers required 
are the angles of the initial and of the final turns and the distances 
between turns for each of the four possible turn combinations. 

If either of the opposite turn combinations (i.e., left initial, 
right final turn) has overlapping turn circles, then that particular 
combination is eliminated as a candidate for the reference trajectory. 
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Determipe altitude profilaa . - In planning the flight to the FAF| 
some margin of flight path controllability ahbuld be maintained for 
control againat errorsf i.e., wind, aero asBumptione, eto. This 
requires a nominal angle-of -attack ( a) lees than that for L/Una-r. To 
maintain thia control margin would require the nominal a to be con- 
stant. The guidance plans each segment for a specified constant a, 
and assumes the vehicle will achieve a steady-state flight path ( >S8^* 
Examples of numbers used with the straight-wing arbiter aret 



Initial 

Turn 

Uings-level 

between 

Turns 

Final 

Turn 

Final Approach with 
Landing Gear down £ 
Soeedbrakes • 50 % 

“N 

7.5 

4 

6 

4.5 

B 

8 

7.2 

a 

10 


Justification of this assumption and equations that determine magni- 
tudes of 7^^ for data load into guidance computer are presented in 
( appendix B. Equation four, which specifies the velocity required for 

7 to be constant, and equation seven, which predicts the trajectory 
with any velocity magnitude, are utilised in the guidance. For the 
initial turn, the ground plane distance X = (Turn Radius * Initial 
Turn Angle). The predicted trajectory would be 



The incremented sh5ft in predicted altitude (H) accounts for the 
general case wherein the existing vehicle energy state is not exactly 
that required to fly at the new angle-of -attack for the next segment. 

The guidance predicts an appropriate shift in altitude for an assumed 
constant 7 segment. The actual trajectory, of course, approaches the 
predicted trajectory smoothly and reaches the same end conditions, very 
closely. 

Equations four and seven are reprocessed with the predicted end con- 
ditions of each phase to predict the entire altitude profilet 
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Present Altitude 



The minliBUffl entxy point is a point in FAP which is a specified 
arbitrary time) i.e., 60 sec, before touchdown on runway. 

Select most efficient maneuvers . - Out of the four possible sets of 
turn combinations, the set that produces the highest projected altitude, 
at a given location of FAP entry point, is selected as the mosb effi- 
cient. The projected altitude error of this set is then used to compute 
a new FAP entry point. 

If the FAP entry point does not require adjustment, then this 
selected set will determine the turn direction commands. 

Ad.1ust FAP entry point . - If the moot efficient projected altitude 
error is within a deadband range (+1CX3 ft has been used}, the entry 
point will not be adjusted. Otherwise, the adjustment is made according 
to the guidance equations derived in appendix C. The lo^p control logic 
shown on the guidance flow diagram of Hgure 3 will prevent infinite 
looping in the guidance. 

Convergence failure . - The guidance could fail to compute a satis- 
factory reference trajectory as shown in the logic flow of figure 3, 
for two reasons. If the vehicle were outside the footprint of capa- 
bility with nominal flight paths, a landing might still be possible by 
flying the most efficient trajectory with L/P ma-r rather than nominal 
values. Another type of failure can occur if the vehicle were too high 
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to fly directly to the mlntmum entry point without diving (> > 15°-20°), 
and too low to fly any other two-tum maneuvers, A three-tum solution, 
as sketched on figure 2d, is achieved by oommanding an initial turn 
opposite In direction to that of the most efficient trajectory to the 
minimum entry point, until a two-tum solution is achieved. 

Reference trajectory . - The reference trajectory is basically the 
most efficient altitude profile previously defined. More specifically 
though, the ^idance output at this point is merely the altitude (Hq) 
and rate (ElcI commands for the present instant of time. During the 
initial turn, 

He = H + SHIFT 

He = -V • Sin >33 


l.e., the energy difference between present 
and desired conditions expressed as 
an altitude shift and where is 

the Vq of equation (4), appendix B. 

The initial turn is made at max L/D (although this is not a constraint) 
to conserve altitude during turn and, therefore, there is no control 
margin to remove residual altitude error (that error not removed by 
entry point adjustment, DH), but after the initial turn, 

H^ = H + SHIFT - DH 

He = -V • Sin >sa. 

The lateral part of the reference trajectory is the ground track 
that goes with the most efficient altitude profile. The lateral command 
is by definition the actual position of the vehicle, because the ground 
track is defined from the present vehicle position. For the initial 
turn, it is necessary only bo fly a constant radius turn in the com- 
manded direction in order to maintain the same projected ground track, 
and for the flight between turns to fly in the commanded direction 
toward the final turn. There is an exception, though, for the final 
turn, during which time the entry point is no longer computed; i.e., 
remains fixed. The lateral position is then commanded to be a fixed 
radius from the center of the turn. 

Initial tumt The lateral axis definition of the reference 
trajectory Is chosen as the perpendicular to the true airspeed (VA). 


where 


SHIFT = 


2 2 
V - V 

c 


2g 
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Y, = Y = 0 



In a following section, an acceleration command, Yq will "be com- 
puted Independent of Y and Y to hold a constant radius turn. 

Flight between turns: The lateral axis definition of the 

reference trajectory is chosen as the perpendicular to the selected 
most efficient flight direction vector (TO) as computed In appendix A 
(l.e., FD = or AD, etc.). 

Y^ = Y = 0 

io = 0 

Y = Vk • g U) 

where LV Is the local vertical and U( ) means the unit of the 
vector within the parenthesis. 

Final turn: The final turn can be flovm In the same semi- 

ppen-loop manner as the Initial turn. But for more precise guidance in 
the presence of errors, tlie lateral axis direction la chosen along the 
vector between the present vehicle position (^) and the present loci- 
tlon of the center of the air turn circle v TRNCNT ) of radius R. During 
the final turn, TRNCNT Is a fixed vector only if the aprlorl wind Is 
aero (discussed In later section). 

Y„ = 0 

Y = R - I T RMCNT - RP 

Y, = 0 

Y = VA • U f TRMCNT - 

Phase control . - Definition of the reference trajectory and vehicle 
posltlnr relative to It has just been shown to be phase dependent. 
Initial turn (IT), between turns (FD, flight direction )> or final turn 
(ft). An estimate of phase time Is necessary for phase control. 

From the geometry (defined In appendix A) of i.he most efficient 
ground track selected by guidance, the ground plane distances for each 
phase can be computed. In the process of predicting (using appendix B) 
the trajectory, airspeeds (horiaontal component) at each phase change 
point have been computed. 
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The guidance trajectory determination asaumea Inatantaneoua attitude 
responae at phase change polntai l.e,, 0 = 35-— ^0, To compensate for 
attitude response time delay, the smrieuvera are led. The phase (PH) 
control logic with examples for lead times ore: 

At first processing of Quid, PH = IT 

At subsequent processing, 

If PK = n and TGO(IT) < 3, then PH = FD 

If PH = FD and TGO(IT) > 5, then PH = IT 

If PH = FD and TG0(FD'> < 3, then PH = FT 

If PH = FT and TGO(FT) < 3, then EXIT Term Area Guid| 

i.e., PH = Final Approach Guld. 

Wind considerations . - The guidance will compensate for wind errors 
in two simultaneous manners. With error compensation networks (to be 
discussed), the angle-of-attack will be adjusted in order to hold the 
reference trajectory. Also, the FAP entry point will readjust during 
IT or FD phases if the vehicle is blown off the predicted course. 

The performance of the guidance in the presence of wind is greatly 
enhanced by providing the guidance with knowledge of the wind (apriori 
wind magnitude and direction versus altitude). The reference trajec- 
tory can be adjusted for wind in order to maintain nominal angle-of- 
attack, in which case, a- margin is allotted for wind uncertainty of 
apriori profile rather than on full magnitude of win*'.. 

In addition to apriori wind profile, a routine is desirable which 
will compute the average wind over two specified variable levels of 
altitude. 

To incorporate this wind in the guidance, appendix A requires modi- 
fication of final turn circle (S) as shown on figure 2c. To keep the 
turn bank angles nominal, the final turn circle center is first rotated 
so that ground speed at end of final turn is in FAP, and second, trans- 
lated against the wind ly the sum cf the average wind of each of three 
maneuvers times the time of each maneuver. 

Guide to the Reference Trajectory 

Each time the reference trajectory equations of the ^previous section 
are processed, vehicle state commands (h^, h^, Y^, and Y^) are generated 
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a n d prosent vehicle state components (h, h, Y and Y) are defined. The 
task Is now to guide to the commanded state. The guidance equations of 
this section are derived in detail and presented in appendix D. A 
functional description of the guidance output commands (two components 
of Euler attitude and three body rates), and guidance-autopilot inter- 
face follows. 

Attitude co"»«Ands . - Pitch (B^) and roll (0^) attitude commands of 

the standard aircraft Euler attitude sequence; ^ 0, «re generated. 
An azimuth command, is not generated. The control requirement on 
^ is that sideslip be zero relative to actual, not navigated airspeed. 

To control the vehicle state, a second-order control law is u ed, 

«o = - 9) + “ 9) + 93s 

which corresponds to a system of the form shown: 


Qfs 



where Q = is either h or Y 

= natural frequency of desired response (of the order of 
10 times slower than vehicle attitude response) 

f = damping ratio of desired response 
" dQc 

Qgg = = steady state acceleration required to maintain the 

two error signals at zero. For example, this term 
for the lateral axis during a turn would be the 
centripetal acceleration. 

The acceleration commands and Yg are converted respectively to 
angle-of -attack Oq and beuik angle 0^^ commands. Pitch 0^ and roll 0c 
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commands are computed from these tvc angles and the navigated airspeed | 
i.e,, VA = ^(Nav gnd speed) + VW (apriorl wind). 


Attitude rate co™>*n<^" - - The state vector acceleration commands 

(Qq) are by definition variable. The vehicle can be connanded to a 

linearly varying acceleration state by qpoputlng body attitude rates 

■which will produce the commanded jerk (O^), where an expression for 

jerk crmmar.d is ob'talned by differentiating Q , 

0 


Qj = Wn^(9„ - Q) + 2tWn(Qg, - q’ ) + Q'.a 


Body rate commands, Pq, and r^, as functions of jerk oomaand are 
obtained by differentiating (as done In appendix D) and 

Error compenaatlon . - Acceleration errors can exist from errors In 
the aero coefficient assumptions, navigation, wind, attitude response 
time delays, etc. Any of these errors that exist In a form of a bias 
can be detected by comparing the present vehicle state with the state 
that should exist had the acceleration and jerk coonands of the prev- 
ious computation been flown, A bias error term Is obtained by weighting 
and accumulating this measured error, compensated for expected acceler- 
ation error due to attitude time delcy. This bias is then fed back in'to 
the acceleration consnand so that the actual acceleration will converge 
on the command acceleration Ir th«- presence of any bias errors. 

Guidance-Autopilot interface . - The guidance can operate at a rela- 
tl-'’’ely large computation cycle time such as 2 sec or probably even 
larger. There are autopilot Interface equations that must be processed 
more often as shown on figure 4. If the sequence of attitude of the 
Inertial measurement unit (IMU) is not s'tandard aircraft relative 'to 
landing site, then the three angles must be converted to the required 
6 and 0. Attitude error Is then converted to body axes error. Coordi- 
nated turn rate commands must also be added to the guidance rate com- 
mands. The autopilot Is shown on the figure simply as a second-order 
system merely *10 show that the guidance commands on the autopilot aire 
attitude ramps in body axes pitch and roll. Body axis yaw oust control 
to a rate and also for sideslip i3 = 0. 0 could be replaced with 

either an Integral of accelerometer measurement or merely a signal pro- 
portional 'to bank angle. 


Selec /ion of Turn Radius 

With respect to altitude lost per unit of turn angle, a 45^ bank 
turn la most efficient. With respect to comfort though, a 30® bank 
turn is more desirable. A constant radius turn tends to satisfy both. 
For example, a 9200 ft radius for the stralght-wlng orblter requires 
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45® bank for an initial turn at 20,000 ft altitude with max L/D, and a 
30® bank for a final turn If at 4000 ft altitude with L/D < 

This same situation exists for the delta-wing orbiter with a 12,000 ft 
radius of turn. These numbers were used for the results that follow in 
another section. At Initial turn altitudes greater them 20,000 ft, the 
bank coanand was limited to 45°. This results In a minor adjustment of 
entry point In FAP because the actual Is greater then guided turn 
radius. The actual converges to the guided trajectory at either the 
completion of initial turn or altitude less than 20,000 ft, whichever 
occurs first. 

Turn radius could bo selected so tliat It would bo flyablo with 45° 
hunk at the wwHwnim altitude for terminal area guidance; l.e., about 
40,000 ft for delta-wing. But this two to throe times larger radius 
would give an Inefficient final turn which would require of the order 
of 15° bank. 


DEMONSTRATION RUNS 
Description of Simulation 

The results presented were made on a slx-degrees-of -freedom 1108 
computer program. This program contains aerodynamic forces as a func- 
tion of angle-of -attack for a Mach number of ,25. In place of aero- 
dynamic moments, an autopilot Is assumed which drives the vehicle to 
the command attitude and rate with a second-order response, where the 
natural frequency Is 1 rad/sec with ,707 damping. The results of this 
program, though, are representative of results obtained on a more 
detailed real-time hybrid simulation. 

Guidance reference trajectory computations and attitude commands are 
processed every 2 sec. The autopilot and environment equations are pro- 
cessed every 1/16 sec. The natural frequency of the guidance commands 
for driving to the reference trajectory Is selected to be 1/10 that of 
attitude response, or .1 rad/sec with .707 damping. 

Phase (IT, FT), or FT) control is monitored in the fast time loop so 
that phase change will occur within 1/1 6 sec of the time computed In 
the slow time loop. The g«ildance equations are processed Immediately 
at phase change. 

This report presents only the guldeuice to FAP. The end condition 
energy state should be such that the reference trajectory In FAP will 
be Intersected after entry Into FAP, The simulation contains a complete 
guidance to touchdown, and the results show part of the FAP guidance in 
order to demonstrate the intersection with the reference teajectory. 
Simulation data applicable to all runs made Is shown on table I. 
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Results 

Results of runs are presented to demonstrate i 

For stralght-ving orbiter . - 

(1) Ground tracks from various initial position relative to runway* 

(2) Effects of wind errors on trajectory and control parameters, 

(3) Effects of very high winds with guidance knowledge of such, on 
trajectory and control parameters. 

(4) Altitude convergence to reference trajectory. 

For delta-wing orbiter . - 

(1) A nominal ground .track from a maximum energy state initial con- 
dition for terminal area guidance. 

(2) Same as (l) except with a guidance modification for a maximum 
entry point. 

Initial txasition . - Run 1 on figure 5, which consists of two right- 
turn maneuvers to an entry point at 50,000 ft range from runway, is the 
nominal run that is used in the next sections. Run 2 demonstrates a 
left turn-right turn maneuver. For run 3, which is already in the FAP, 
the altitude high error at MEP would be greater than 100 ft if the 
vehicle were to proceed directly to MEP with nominal angle-of -attack. 

But also, the energy state is too low for any other two-tum solution 
to MEP; i.e., a 360^ initial turn. The guidance sets the flag lOPOST 
to 1 when it detects this condition, and the vehicle is commanded to 
make an opposite turn; i.e., initial turn direction opposite in direc- 
tion to the most efficient path to MEP. After approximately 110® of 
turn, a two-tum (right turn-left turn) solution is achieved, and the 
flag is then set to zero. For the above case, a possibility under study 
exists to modulate speed brakes and thus increase the flight path angle. 
This technique could also provide a satisfactory two-trm solution and 
minimize the necessity for a non-standard procedure. 

The initial energy state ^ jt run 4 is too low for nominal angles- 
of -attack flight to MEP (lOPOST = -1). The vehicle is then coimnanded 
along the best path, but at max L/D. At a later point, the vehicle 
intersects the reference trajectory, at which point nominal angles-of- 
attack are flown (lOPOST = 0) . 

Viiya errors . - Run 1 (no wind) is replctted on figure 6 glong with 
+40 ft/sec wind errors (zero wind assumed). This wind is 45° relative 
to the runway. Entry point varied between 46,000 and 52,000 ft. Alti- 
tude is shown to converge on the reference trajectory in FAP. 
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Control parameter variation la shown on figure 6b for the nominal 
€uid one wind error case. Angle-of -attack holds very near the refei^nce 
values for each phase (see table I for comparison). 

Anriori wind . - With guidance computer knowledge of the high winds, 
a successful descent was made as shown on figure 7. Note the altitude 
climbed during the initial turn, the reason being that this run started 
as run 1 with seune ground speed, and in effect, it instantaneously 
encountered the wind which then increased the Initial airspeed by the 
wind velocity. Note also the change of final turn direction from run 1 . 
This change of direction was net one which occurred after Initiation, 
but. was computed as the best during the first pass through the guidance. 
The turns no longer appear circular because the plot la relative to the 
ground. Relative to the moving air mass, the turns should remain as 
circles. 

One thing the guidance does with its knowledge of wind is adjust 
the reference trajectory slope in FAP in order to maintain the nominal <>< 
for that phase. This slope adjustment is shown on figure 7a. The 
altitude la shown to converge on the reference trajectory in FAP. 

Control parameter variation is shown on figure 7b. Angle-of-attack 
converges very near the reference values for each phase (compare with 
figure 6b). 

Altitude convergence . - The guidance computes altitude command as a 
function of the potential energy equivalent of the kinetic energy dif- 
ference between present vehicle state and that required to hold the 
reference flight path with nominal angle-of-attack. At phase changes 
then, there can be a discontinuity in He as shown on figure 8a. for 
run 1. Note, though, for the nearly constant a (see figure 6b), H con- 
verges to He (figure 8a). The discontinuity at other phase change 
points is small because of the choice of targeting conditions, but this 
discontinuity is not required to be small. 

The same type plot is shown on figure 8b for the apriori wind run. 
Heie there are discontinuities at each phase change, but, still, the 
trajectory is shown to converge smoothly to the reference. The reason 
for these discontinuities on figure 8b and not on figure 8a is that 
the nominal final turn angle-of-attack is replaced with the higher 
initial turn a. The guidance automatically makes this targeting change 
whenever it estimates the altitude of the final turn to be greater than 
15,000 ft (13,000 ft for delta-wlngK This keeps bank angle required 
from exceeding the maximum in order to hold turn radius. 
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Noainal delta-vine orbiter. - Fop a oonatpalnt of Mach < ,9 during 
the Initial turn with a = 7.5''f there Is a marlmia energy state from 
which the initial turn oan be started. The runs on figure 9 are appror 
imately from this muT-iimim energy state. Note that even though the 

12.000 ft turn radius cannot be flown at h = 41 >000 ft, a smooth con- 
vergence to a solution is achieved. The entry point with the nominal 
guidance is 30,000 ft from the landing site. 

There should not be any fear that the entry point is too far from 
the runway, because there is considerable margin in controllability 
against errors once in the FAP. For example, at the 80,000 ft point 
for the nominal trajectory, the landing gear gees down and speedbrakes 
at one-half maximum drag. And even if this were not enough control 
margin, then the reference trajectory slope to MSP could be made 
itoeper — runs were made with the same slope to and from MEP. But, to 
demonstrate that entry point can be limited, the dashed run on figure 9 
activated the opposite turn logic whenever entry point greater than 

60.000 ft. 


CONCUJSIONS 

A guidance technique is derived which will guide an unpowered air- 
craft in a nearly optimum manner from any bearing and heading with 
respect to a runway to a reference trajectory in the final approach 
plane from \diich a landing can be made. The entry point into this 
plane is computed so that the aircraft can fly at nearly constant, 
specified angles-of -attack. Constant turn radii are utilised. 

Simulations of both straight and delta-wing orbiters denonstrate the 
ability of the technique to guide in the presence of large known winds 
and of wind errors. 

Subsequent study will be concerned with procedural variations within 
the guidance which might be desirable for use under different weather 
and wind environments. 
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Altitude 400,000 ft 



Time, minutes 


Figure 1 Guidance phases. 
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Air speed 




(b) Select best maneuvers out of four possible sets. 



(c) Displace arget circle by function of known wing. 



(d) Three turn maneuver capability. 

Figure 2. • Sketches demonstrating ground tracks considered by guidance. 
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Figure 3. -Flow diagram of guidance to final approach plane. 
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Figure 4 Guidance-autopilot interface. 
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Figure 5.- Terminal area guidance of straight wing vehicle from various 

positions relative to runway 
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Elapsed time, sec 

Figure 6b.- Terminal area guidance of straight wing vehicle, nominal and with wind error. 

Trajectory control parameters. 


9.15-22 




Final 

appr 




Figure 7a.- Terminal area guidance of straight wing vehicle wivh apriori wind 

equal actual wind. - Trajectory. 
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Figure 7b. ^ Terminal a a guidance of straight wing vehicle with apriori wind equal 
actual wind.- Trajectory control parameters. 
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Figure 8a.- Terminal area guidance of straight 
wing vehicle/ nominal - altitude and command. 
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Figure 8b. - Tenaln'\l r rea guidance of atralght-wing vehicle 
with aprlorl wind equal actual wind - altitude and command 
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Figure 9.- Terminal area guidance of delta wing vehicle. 
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Ground Track Vector Equations 
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a = ASIK ( 


2R 


') 


|D - A| 

AD = (Is - Al- 2 • R • SIN a) • U (D - A) + (2 * R • COS a) 
. U ((D - A) X R)) 

If |D - A|‘^2R, se'« flei ‘:o eliminate this turn comb. 



BC = (|C - B|- 2 ‘ R • SIN b) * U (C - B) + (2 • R • COS b) 

• S (k X (f: - B) ) 

If |C - B |^2R| oX^Bii.nA^G ^uth cosib* 

Initial Turn Angles - Solve for + angle (O to 360°) between V and 

vector between turns. 

Final Turn Angles - Solve for + angle (0 to 360°) between vector 

between turns and 1. 
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Trajectory Prediction of Aircraft Flying Constant 
Angle-of-Attack and Turn Radius 


Steady state flight path an^de . - The equations presented hsre are 
not (except equations (4 and 7)) a direct part of the guidance, but are 
equations to be externally processed for data load In the flight plan 
portion of the onboard guldemce. The guidance makes simple use of the 
fact that the steady state flight path angle Is nearly constant for 
constant angle-of -attack . The approximation 


(y = xTScos"^^ satisfactory. 

To hold a constant turn radius (R), the centripetal acceleration, 


V - Ccs7 

but lift. 



L = -^ ClS pV^ 


• • » 


(1) Sin 0 = 


2M Cos^ > 


RClSp 

and by differentiating (1) 

^ ^ . Tan 0 

Standard atmosphere air density, 

1 

(2) P = .CK32378 e ^ " 6.875 • 10 


(3) 


by differentiating (2) 

-i _ 6.875 • * V • Sin T 


.235(1 - 6.875 • 10 


-6 


-6 


h) 


h) 


If there Is a steady value of > , then 

-y_Q_g* C 08 >-j;j cos 0 

- 


SIGN CONVENTION 
y (+) for Descent 
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or 


L ri 

g COO ^ ^ 00* P 


fros this *rd the lift equation velocity lo conelrained. 


(4) V = 


- | 2 • W • g~ 

Vc, • s 


C 08 > 


P • COO 


T 


The rate of change of velocity by definltionj 


V = g * oln 7 - g 


or 


g 3in 7 = g + V 


Also, for a steady 7 , 

e < 

7=0 


or 


d(L • coo 0) = 0 = L Coo 0 - L Sin 0 • 0 

* • • 

i = L- -^{1+2i|) = LTanfJ'?l = -L' Ta^n 0 ’ y 
(5) V ^ • 7 (1 + T»n^ 0) 


The flight path that producea 7 = 0 lo then^ 

_ 1 . V 


~ + V 

ff Sln7 M 

Tan y = v t r = 7 


g Coo 7 


g Coo 0 g Coo 0 


Sin 7 


g Cog ' Sin 7 


by rearranging and collecting Tan 7 terns, this equation becomes 

Ivi Iterative procesoing ''f equations (1 to 6) provides the state of 
vel^ ..\ty and flight path that muot ''xlot for the first and second deriv- 
ative of flight path to be zero at any altitude, radius of turn, and 
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angle-of -attack. Two questions to be answered empirically are: 

a. Is ther '3 correlation between the predicted *y of equation (6) 
and actual 7 ? 

b. Is there a steady state value for 7 ? 

Results of equations (1 - 6) for three conditions of or and R are 
shown as a function of altitude on figure B-1 . The aero-coeff iclent 
used are those of a straight-wing orblter. The wlngs-level case 
(Infinite radius) corresponds to within the resolution of the plot with 
an actual flight. A steady state value of 7 is not achieved, but for 
practical guidance purposes ^ the variation Is sufficiently small so 
that a constant 7 can be assumed. For turgs under 20,000 ft altitude, 
If a and R are selected so that bcmk < 45°, the 7 holds fairly 
steady at 8°. The simplified estimate of 


^ Cos 0 

Is shown to be greatly In error. 

Turns starting at higher altitudes and with turn radii so that max- 
imum banks are 46.5 degrees are shov^i on figure B-2, along with actual 
results for c.^e case. More variation is seen here for 7 , but still, 
for practical guidance purposes, a steady state 7 can be assumed. For 
example, at 50,000 ft altitude, a 180® turn would lose 8000 ft altitude 
with a resultant 7 from 4.4 to 5.75. The guidance can estimate alti- 
tude at end of maneuver and 7 (5.75) for that h can be assumed. The 
angle -of -attack required to hold this assumed 7 for o^= 7.5 will be 
< 7.5. This is a conservative approach because the kinetic energy at 
end of maneuver must be > aim condition of guidance. 

Tra.1 ecto^ prediction . - At any altitude with a velocity as given 
by equation (4), the future flight path can be predicted and It Is suf- 
ficiently constant so that a linear trajectory (7ga) can be assumed. 
The remaining problem, though, la to predict the trajectory from a state 
with V ^ equation (4). A solution fs to shift the linear trajectory by 
the potential energy equivalent of the kinetic energy difference of 
vehicles present state and that of equation (4)| 

^ (V^ - V^^) = WAh 
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= V or equation (4) 
W = Mg 



2g 


the pradicted referonce trajectory for constant a 
K = h . - I . T« > 

‘C “ 2g * 




flight 


Soma results of this assumption are shovn on figure B-3. A of 

'^.16 was selected from figure B-1 for an altitude of 10,000 ft. The 
cnree actual trajectories of figure B-3 were started at h = 20,000 ft, 
y ~ 6,8, and V = 579.3 (from equation 4)> 632.4 and 989.5 coiresponding 
to 0, 1000, and 10,000 ft equivalent potential energy difference. The 
ft case steadies out on the reference trajectory, but the shift 
for the 10,000 ft case is only about 50 percent of pr^icted. 

Energy is dissipated during flight because of drag force (D). 

AE = J" D dx 

The trajectcry prediction assumes no additional energy dissipation over 
this nominal dissipation, but actually, there is additional drag with a 
greater velocity, and therefore, 

AE = J (D + AD) dx 

It is possible to predict this additional energy loss to better predict 
the trajectory but an iteration solution would be required. Also, this 
better solution indicates that the 1000 ft case is only 25 ft in error, 
and 1000 ft is about the maximum value encountered with the guidance 
system presented in this report, and therefore, the simpler prediction 
(7) is valid. 
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Figure B-1 Predicted flight path for constant angle of attack flight. 
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Compute New FAP Entry Point 


Gl^sn: 1. A vector DBA from exit point of Initial turn to entry 

point of final turn. 

2. Exit point of final turn (R6X = FAP entry point). 

3. TAN of flight path angles; to FAP and K 2 In FAP. 

4 . Projected vehicle high altitude error (DH) at minimum 
FAP entry point (R5X). 



poiN'T 


Assumptions: 1. A new DDA^ solution for a new RSXff will have tho 

same DDAY component. 

2. The sum of both the Initial and final turn angle 
remain constant with variation of R6X, and there- 
fore, altitude loss of turns Is approximately 
Independent of R6X. 

NOTE; These are better afsumptlons for turns with same direction 

(l.e., left li^tial and final turns) than for opposite turns 
— thir should merely mean that convergence on proper R6X will 
be faster If vehicle Is In position for alike turns. 

Because of the assumptlors, the probltj becomes one of determining 
’,ntersectlon(s) of two straight lines — trajectory to FAP and In FAP. 
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RaiMT R •»-* 

s&Ay 

FflP 

f?^.Xvs R^X Wpy 

Q, *■ y 

= R5X - R6X + DDAX 

The constant altitude difference between point A and the desired 
altitude ® R5 is 

AH = >/dDAX^ + DDAY^ + K2(Q., - DDAX) + DH 

The objective of the new R6 Xjj Is V have a zero DH, therefore, 

K., VdDAXjj^ + DDAY^ + (ij.j - DDAIn)K2 = K., V DDAX^ + DDAY^ 
t K 2 (Q., - DDAX) + DH 
rearranging equation 

(1) K 2 VdDAX^^ + DDAY^ = K 2 • DDAXj, + 
where 

(2) 
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By squaring equation (1) 

K1^ (DDAXjj^ + DDAT^) = (K 2 * DDAXjj)^ + 2Q^ • K 2 • DDAX + 

which is of quadrrliic form 

A * DDAXn^ + B • DDAI + C = 0 
where 

A = 

B = - 2 Q 2 K 2 
C = (K<| • DDAY)^ - 

the solution of w hich is 

-B ±\/b^ - 4AC 


DDAXn = 


2A 


this equation presents a problem though, when A = 0, which is 
likely for a no-wind condition; therefore, a better form is 

(-B + \Ib^ - 4AC) 

(-B + - 4AC) 



DDAIjj = 


-B +\B^ - 4AC 
2A 


Before using equation (3), some constraints must be applied . 
The roots of (3) must be real, 

0 

I I ■■ii I II i,md 

1 - 4AC/B^ > 0 
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There are at most two real roots of equation (3)» and there can be 
eight sets of two real roots depending on the signs of the coefficients 
(Ay By & C). Four of these sets where B 1s positive, although real for 
equation (3) are not real for the physical problem; l.e. , equation (1). 
These can be avoided with constraint 

B < 0 


or 



Two other sets of roots where A Is negative contain tv' real roots 
for equation (3), but only one Is valid for equation (1), The valid 
one Is obtained by using the positive sign In equation (3). 

The remaining two sets are valid for both equations, but the more 
negative root (FAP entry tt greatest distance from R5X) Is desired, and 
It also r squires positive sign In equation (3), and therefore. If no 
constraint violated. 


DDAX^ = 


-2C/I 


1 Vi - 4AC/B^ 


R6Xk = R5X - Qi DDAXjj 
.tf (R6 In > R5I) R6X = R5X 


If any constraint was violated, then FAF entry point set to minimum. 


R6X = R5X 
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Guidance to Reference Trajectory 


• e 

The reference trajectory (H^, Yq) and actual vehicle etate 

(H, H, I, Y) have been defined for vertical and lateral directions. 

The guidance equations are derived that produce guidance output coa- 
■ends of of standard aircraft sequence 0 and body rate 

comands Pq, qQ, r^. A compilation of just the guidance equations is 
suraarised at the end of this appendix. 

General Physical Eouatione Utilised 

Airspeed 

V = V (Nav relative to ground) - W (xpriori wind) 

V = /!/ 

Dynamic pressurs 
q = 


Lift coefficient 

(1 ) “ 

Uft 

(2) L = Cl S q 

Air3p««! ftcs.l.ratlon 
^ = (V - v„)/at 

wh.r. 1. ccBputatloa Gyoli> tla. 2 Mc) 

and Vq is airspeed at previous cc«putation time. 

Flight path 

7 = Tan'** (VZ/ Vvj^ + Vy^) Note* >+ DOWN 
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Drag 

(3) D = g sin 7 - V 

M 


Total vertical acceleration 

(4) % ” M ^v Coa>+ ^ Sln7 - g 

where 

0^ = bank angle with reference to V. 
froB (3) and (4) 

e • 

(5) L g w _ ^ V Sin 7 + g Coo 7 

M ~ Cob 

a component of vertical acceleration defined as 




(6) H = Hx + V Sin 7 

angle of attack from (1, 2, 5; and 6) 

vf " X 

(7) a = il-a (• -H + g Cos 7 ) - -r— 

' ' pCLA • S • q • Coa 0^ P^La 


by differentiating (7) 


(8) a = 


M 


• • • 

(^.*y 


00 

H 


pCU . S . q 7 Cob 0^ \ Coa 7'*' ^Coa7 ** 

- (chy * * + T - ^y>) 


also rate of change of angle-of -attack 


(9) a = 7 Coa 0^ + q 

where q la body pitch rate. 


>0 I 4^ ' 


g Cos 7 ) 
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rate of change of flight path 

• g Cob 7 - ^ Cos 0 ^ 

~ V 

fin) y = — =S— 

V Cos 7 

lateral acoelaration 

(11) i' = ^Sln0^ 


Bank angle froa (5> 6, and 8) 

(.21 0 =Tan*V'-B ^ ^) 

\( ;j- ^ g Cob > / 


Cob 7 


by differentiating (12) 

(13) 0 g - -- — (f-Tan0 (=^+7TBn> 

(csfr - 8 Cos-y))^ 

Vertical ailB guidance . - To guide the vehicle altitude to the coa- 
aand reference trajectory with apeclfled reanonse charaoterlatlcB of 
natural frequency (VN) and daaplng ratio ( T )> a aecond-order accelera- 
tion ooBwand la generated 


’^T = - H) + 2 rw„(Hg - H) - V Sin 7 


The last tera of this equation la the ateady atate acceleration 
required to hold the error algnala of the first two t^raa at sero. But, 
If this acceleration la applied to equation (5)f the V tera will cancel, 
and therefore, just the following coaponaQt of the ooBoand need be 
generated . 
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(14) 


Ho = Wn^(Ho - H) + 2 f W^(H^ - H) 


(15) 


Eq[uatlon (7) will give the resulting for this accel command, 

M 


“c 


PCLA * S • i • Cos ^Cos 7 + 8 Cos >) - 


®c . Cio 


PCLA 


where 0^^, %*lch is discussed later, must be generated first for 
use in this equation. 

Given 7 , «c» ^ c' local vertical component of a unit 

vector .^long the comnanaed direction of the fuselage is, 

(16) UXBLV = Sino^ • Cos 0^0 * Cos7 - Coso^ • Sin 7 


The pitch component of this vector for the standard aircraft 
sequence of ^ 0 0 is, 

(17) -e^ = Sin"** (UIBLV) 

The acceleration command (14) is a variable, the slope of which 
at a given time. 


He 



(He - H) + 2fW^ - H) 


% = -V Sin 7 


dt 


ss 


(as defined for reference trajectory) 


assume. 


H = Hgj = Hq + V Sin 7 


• • » 


He = V (He - H) - 2fWn (He) 


• i 
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Substituting this into equation (8) and neglecting the last two 
terms of (8), which is an Imperlcally determined good approximation, 


(19) 




M 


PCLA 


^ 

q . Cos 'Cost 


From (9 and 10), the pitch rate command, 


( 20 ) 


H. 


Cos 


0 . 


9c = 


VC 


'.OS y 


+ 


9 



If these Cvommends were perfectly applied, then the total vertical 
acceleration would be. 


.. .» . •»' 

% = Hq - V Sin 7 • t 


Assume, though, that there are bias type errors present resul"^ 'ng 
from errors in aero coefficients, navigation, wind, etc., 

the actual acceleration, 


- ftj + H, 


This error can be computed by comparing the altitude rate from one 
computation cycle to the next with the rate that should exist, 

At 



( 21 ) 


0 • 

H. 


^ * 2 / 

= (H - Ho - • At + V Sin 7 • At - -jHco 
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The objective la to remove blaa errora, but not tranalent error auch 
ea caused attitude reaponae time delaja. For a atep change of attitude 
coanand, the Inetantaneoua error would be. 


( 22 ) 


'u i - PCLA * 
- 


M 


(«co - 


Aa that ayatem reaponda to the command, the error will reduce, ao 
that the net expected error over /It see will be some fraction (i.e., 
,4) of . The error to be compensated should then be reduced by the 
expected error, 


(23) 


= H, - .4 • 


By accumulating this error through a gain, I.e., .8, 


(24) HbiaS = Hbias + .8 • 

and feeding back a corrected comaand Into (13) only 

(25) Ho = h’o (of (U)) - Hbias 

the bias will b) removed and the actual acceleration will converge 
on the desired level of (14). 


lateral Axle Guidance . - To guide the vehicle lateral position tc the 
command reference trajectory with specified response characteristics of 
natural frequency (Wj^) and damping ratio {t ), a. second order accelera- 
tion command la generated. 


Vo = W„2 (Yo - I) + 2 (Yo 


e • • 

- + I'm 


The last term of this equation la the steady state acceleration 
required to hold the error signals of the first two terms at zero. This 
last term erlsts only for turn maneuvers and la the centripetal force, 
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9 . 1 ? . 2 . 1 Approach Guidance ( cont ’ d ) 


(26) 


y = (V ‘ Goa (radtua of turn) 

8S 


Tha rafarenot traj actors aaotion of tha guidanca haa aalacted tha 
y^xla diraotlon ao that tha poaitlon and rata oonnaanda and are 
aaro, 

(27) I’o = -Wn^ - 2 f W„i + 1,, 


The number of tarma of thla aquation that are utilized ia phaae 
dependant. The full aquation la uaad for tha final turn, but only tha 
la at tarn for initial turn and only aacond term for flight batvaan 
tuma. Tha acceleration and jerk computation phaae dependency ia shovn 
in detail at tha and of thia appendix, but juat tha final turn phaae 
will be diacuaaad hare. 

a 

Tha bank angle cornmnd ia obtained by substituting command values 
into aquation (12), 


( 28 ) 

Given 7 and 0 , tha local vertical component of a unit vector 

along tha commanded direction of tha wings is, 

(29) UYELV = -Goa 7 • Sin 0^^ 

Tha roll component of thia vector for the standard aircraft sequence 
of 4^e 0 is, 

(30) 0^ = Sin"** (-UIBLV/Gos e) 

w C 

a 

The rate of change of I which ia rafarancad to a rotating axis system, 



= Tan 


-1 


Y. 


1 


( 

Coa7 


t g G087) 
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The rate of change of the acceleration coumand 

¥•9 Oe , 90 *9 017 • nnfl ^ M* V 

(31) Yo = - 2 > W„ do - log) + " r '~~ + HoT“») 

Substituting this Into equation (13) end neglecting the negligible 
last term. 


(32) 


Ka = d„ - Tan 0^^ • H„/Cos1') Co8^0,„/(;J^ + g Cos 7 ) 


vc' 'Coa'y 


Projecting this on body axes, 


(33) Pc = 0^c ' “c 
r^ - 0YC * ®c 


As done with the vertical axis, a bias acceleration error can be 
measured and compensated with feedback through the acceleration command. 


(34) 


V • • • • • 

y, = (y - Yc - Y 


“• +2 

At - 


(35) 


YBIAS = YBIAS + i Y* 


where the gain of ^ Is empirically determined for stable removal of 
errors 

(36) Yc = Yc (of (27)) - YBIAS 


This compensated value of acceleration commcmd Is used only In 
equation (28). 

The guidance to the reference trajectory equations and logic are now 
summarized for each axis. 


\ 
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VertlQAl Axle Guldmoe Eouiitiong 


Hq and Bq fron reference trajectory section of Guidance 

J. 

e = .002378 -235 1“ d ' 6.875 • 10-6 . H) 


- „ Flight Path Sign Conv . 

(+) Down 

V = (V - Vo)At, where d^t is comp cycle time; i.e., 2 sec 


Vo « V 

\ = W„2(H^ - H) + 2 f U„(Ho - H) 

He = M„2 (Ho - ii) - 2fW„(H^) 

= (H - - Hc^At + V • Sin 7 • At - iflcoAt^)/^^ 


• • 



•• •«! , ^ 

^C» ^00 ”” ^c* ^ ® 


h ' = ^cu • s • q . (0 . 

H 

«o = e 

• I I « • • / 

» .4 • H« 

If < 0), then = 0. 

If (Not at a phase change, and oq < a^mm^) , 


then, HfilAS - HBIAS .8 



He = - HBIAS 


Note: Compute 0^^ from lat. 

guidance first for use 
here* 
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M (Hq/Cos ^ + g Cob y ) Clq 
S 5 PcLA * ^vo " PCU 


If (Flight director phase between turns and nax L/t) discrete from 
riferance trajectory section of guidance), then oq = “cmax* 

UXBCZ = -Sin Oq • Coo 0^^ Coo y + Cosofe • Sin 7 

= Sin*’*' (-UXBC2) 

^co • 0V C + 

% V CosT' 




J Limit to a 
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Utarftl Axia Quidmot Eouatlone 

« • 

^ 0 * ^ 0 * ^ r«f«r«no« tr«j«otorr seowion of guidaaoe 

Umlt Y to ±2000 ft Flight. Path Sign Copy. 

(+) DOWN 

• • i»i 

IpS » ^S8 • 0 

If (Icitiai or final turn), than * (V • CoaY)^/R, 

^aa * (V CoaY + Hqo TmY ) 

Ic * (T^ - I) + 2 rWo (io - Y) 

**• tt *«• 

= -2 fU„ . + I„ 

• • • • ^ t 

Yc = Ic + ^aa 


If (final turn or in PAP, and / I /<2000 ft) 

2 


» • •« < 


than, To = Yo - . Y 


V. = (i - io - Io • At - ^^)/At 


In * I 


a « • 4 



To * To 

YBUS = YBIAS + i k\ 

If (initial tun or flight diraotor phaaa) YBIAS - 0 
Yc * Yo - YBIAS 


If (Initial or final tun and a laft tun ooimand froa rafaranca trajao- 
tory aactio.n of guid) 
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»» .< *“ 
thwi, Yg = -Tq and 

If (initial turn and an oppoalta turn oouaand for thrae-tum aanauver), 

*« • • **• 
thw, = -!<, »ai = -Yg 

Llait to 

ulaCZ = Sla * Css y 

0Q = Sin“^ (DYBCZ/Coa 0^) Notat Compute 6 q first, then use here. 
Kc- (Yq - Tan • H^j/Coa y) Goa^ 7 + g Coa> ) 

Po = Kc * ®c* **0 = * Sin 
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9.15.2.1 Approach Guidance (cont'd) 


TABLE i. - SIMULATION INPUT DATA FOR DEMONSTRATION RUNS 



VEHICLE 



STRAIGHT WING ORBITER 

DELTA WING ORBITER 

Ihltlal 

Condi- 

tions 

H 

20000 ft 

41000 ft 

vg 

562 ft/aec 

703 ft/aec 

Y 

7 deg 

6 deg 

Quid. 

Know- 

ledge 

of 

Trim 

Aero. 

Term. 

Area 

.1954 H 
.1954 ^ 


-.026746 + 1.958083 • 

Final 

Apgr. 


-.014995 t 1.936128 • 

S 

1175 ft^ 

6650 ft^ 

M 

4041 slugs 

8330 slugs 

Guid. 

Target- 

ing 

mm 

n 

IT 

FD 

FT 

^ 

Appr. 

1 

Ln 

i u 

r i 

— 

Approach 

®N 

7 . 5 “ 


6° 

4.5° 

m 

5 ° 

6° 

6.3° 

Yss 

8° 

m 

8° 

10° 

7.46~ 

7.83° 

8.22° 

14 . 5 ° 

Speed 

Brakes 

0 

0 

0 

30° 

0 

0 

0 

20° 

Landg. 

Gear 

Up 

Up 

Up 

Down 

Up 

Up 

Up 

Down 

Turn 

Radius 

9200 ft 

12000 ft 


X 

-21000 ft 

-22000 ft 

h 

2600 ft 

4200 ft 

other 
Quid . 
Items 

Min H 

c 


g 

Max<y 

c 

8° 

8° 

Rsx 

46 ^ 

46° 

R 6 

which 

(Ux 

for 

turr 

i 20000 ft 

i 20000 ft 
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9.15.2.2 Final Approach Guidance 


1. Introduction 

These equations are submitted as candidates to fulfill the 
unpowered Final Approach Guidance requirements for the space-shuttle 
Orbiter, They include Autoland lateral and longitudinal guidance 
equations. The scheme is all-inertial j navigation aids are used 
only to update the navigated vehicle state. Pitch rate and speed- 
brake commands are computed and issued to control in-plane approach. 
Lateral position error and its integral plus heading-angle error 
are used to form the vehicle roll command. (There is no decrab or 
wings level manuever; the assumption is made that the gear is de- 
signed to accommodate the stress for crabbed landings in design 
winds.) 

2. Functional Diagram 

A general description of Autoland Guidance is conteiined in 
Ref, 1. Figure 1 is a functional diagram. Fig. 2 is a block 
diagram. ^For general iniormatlon, the autopilots being used in 
simulation runs are included in Fig. 2.) 

Inputs to the Guidance module (OG-7) are from the Final Approach 
and Guidance Navigation module (0N-5)j the inputs are the navigated 
sta\'.e in the Earth-fixed landing coordinate system. From this are 
calculated the range to toi’chdown target, altitude, inertial ve- 
locity magnitude, inertial flight-path angle, lateral position and 
heading angle. Outputs are p|.tcb rate command, speed-brake position 
command and vehicle roll command to the autopilot. The guidance roll 
command drives a roll-rate aileron-autcpilot inner loop with roll 
attitude outer loop. Roll rete command is interconnected to a rate 
command rudder autopilot with turn coordination and normal ac- 
celeration inputs. The acceleration and heading-angle signals are 
instrumented in holding the orbirer to the final approach plane in 
crosswinds. 
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9.15.2.2 Final Approach Guidance (cont’d) 
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Figure 2b Block Diagram, Lateral 


9.15-62 








9.15.2.2 Final Approach Guidance (cont’d) 


3. Coordinate Syatem 

The autoland guidance uses vehicle position and velocity re- 
lative to a runway coordinate system, as shown in Figure 3. Figure 3 
also indicates longitudinal sign convention for the equations. The 
"altitude of the IMU" at touchdown is represented in the equations 
as c.g. altitude. 

4. Souations and Flow 

4.1 The longitudinal guidance equations are presented in 
Figure 4. 

The guidance routine (as currently implemented) is entered 
four times per second, (This frequency is to be verified by planned 
guidance/autopilot interaction and wind gust response studies.) 
Stability of the guidance equations and minimization of the vehicle 
rotational motions set the minimum rate. There are no switching 
requirements in the scheme. 

The vehicle’s state is in a landing coordinate system. On 
the initial pass, the target sink rate, touchdown ground speed, 
and shallow flight path angle are used to compute reference tra- 
jectory constants. The reference trajectory subroutine is entered 
with iihe targeted range at piillup ccmipletion, and flare and touch- 
down constants, to find the reference altitude at puUup completion. 
At TARGET, a steep flight path angle is coonputed which will carry 
the vehicle from its current altitude through puHup onto the 

shallow slope. This capability may also be used to replan the 

/. 

reference trajectory during the steep phase to minimize transients 
from large navigation updates or to replan when forced by emergency 
conditions. A linear reference velocity is defined to carry the 
vehicle from its current velocity to its targeted value at the 
completion of the pullup. 
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9.15.2.2 Final Approach Guidance (cont*d) 


On a normal call, after the vehicle state is determined, if 
phase is steep the reference trajectory is defined through pullup. 
Vnien the pullup-complete range has been reached, pnase is set to 
shallow and a new reference velocity profile is defined to decele- 
rate the vehicle to its targeted touchdown velocity. The two 
velocity slopes approximate the e3q)ected vehicle deleleration 
along the steep and shallow glide paths with gear down and half 
speed brakes. The switch to shallow is not critical since the two 
portions of the reference trajectory overlap smoothly. 

At LOOPS, velocity dependent gains are computed, A pitch 
rate command is formed add limited. A speed brake command is 
formed such that zero error gets half brakes and both commends 
are issued to the autopilot. 

The reference trajectory subroutine is entered with one of 
two sets of targets; pullup over and touchdown. Calling arguments 
are as follows: 

V - Current vehicle velocity, used for |T . 

R - Range to go to end of current maneuver. PuUi^, 
R-R|j^and touchdown, R, 

p - Constant which sets maneuver curvature, normal 
acceleration level. ■ 

m - Negative of difference between flight path angles 
entering and leaving the maneuver. Steep minus 
shallow and shallow minus touchdown. 

t, 

mb" I^B.gative of flight path angle leaving the maneuver. 
Shallow glide and touchdown glide slope. 

hb- Altitude at which maneuver is complete and vehicle 
flying -mb . Altitude at pullup complete and eg 
height at touchdown. 
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9.15.2,2 Final Approach GuJ.dance (cont'd) 


hj. - Reference altitude 

Vr “ Reference flight path angle 

# 

- Reference flight path angle rate 

ERF is an approximation to the error function devised in 
Ref. 2. Its use is illustrated in Fig. 4. 

Linear approximations of the system transfer functions are 
presented in Fig, 4* 

4.2 The lateral guidance equations are presented in Fig. 5. 

On the initial pass, the roll gain, crossrange integral gain, and 
the heading gain are stored. On a normal pass the crossrange gain, 
Ky, is calculated as a function of velocity. When altitude becomes 
less than 50 ft, the roll command gain is decreased from 1 to .5 
over a 2-second period. The roll command is the sum of a cross- 
range, integral of crossrange, and velocity heading angle term. 

It is limited and issued to -he autopilot. 

5. Constants/Variables Summary 

Figure 6 summariaes variables and constants. Precision and 
quantization requirements have not been finalized. Figure 7 
describes the principal variables. Figure 8 lists typical input 
values. 

6. FORTRAN Coding 

The equations have been coded in FORTRAN for simulation cn 
the SSFS (Space Shuttle Functional Simulator), 

7. Simulation Status 

The scheme is operational on the SSFS. 
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9.15.2.2 Final Approach Guidance (cont‘d) 
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Figure 6c Constants/Variable Summary 
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^.15.2.2 Final Approach Guidance (cont’d) 


h 

Y 

V 
a 


Y 


L 



Vehicle altitude, ft. 

Vehicle lateral displacement 
from runway center line, ft. 

Vehicle velocity, fps . 

Velocity heading angle 

V 

tan a = - ^ , rad. 

Lateral displacement ini’egral 
threshold, ft. 

Roll command, rad. 

Roll command limi , rad. 



■"yi* 


Guidance gains, rad/ft, 1/sec, 
rad/rad. 


Figure 7c Variable Description 
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9.15.2.2 Final Approach Guidance (cont'd) 
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9.15.2,2 Final Approach Giildance (cont'd) 


8, Testing Statua 

A touchdown dispersion study has been performed. Some results 
are shown in Figure 9. For nests beginning at 1740 ft. altitude 
and 15,000 ft. range, crossrange dispersion in +10 ft. and crocs- 
range rate dispersion is +2fps. Vehicle yaw at touchdown is d^- 
tersiined essentially by the crosswind. A 20 -knot wind results 
in a 7.2 -degree yaw. 

9. Derj-vation 

Lateral guidance is basically an all-inertial version of Ref, 5. 

10. Assumptions 

NR 161C delta wing aero and mass properties (Ref. 4). 

Unpowersd landing 
Error Models (per Ref. 3). 

11. Further Activliy 

The integral-of-crossrange term in the lateral guidance limits 
steady state errors. Unfortunately, it also retards response to 

dynamic errors. Simulation runs indicato that in^roved performance 

■' c. 

may result from freezing(claiig)ing) the integral at some point before 
touchdown. 


9.15-^! 



15-82 


I'igure 9. 'I'ouclidown Dispersions 
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MSC FORM 376 (Kiv Aur, 65) 
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9.15.2.2 Final Approach cont ’ d ) 











9.15.2.2 Final Approach Guidance (cont'd) 
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